This notebook is for performing downstream analysis of CYCLOPS. It includes running cosinor regression and assessing differential rhythmicity (two ways for both), calling the enrichR API via python, using fgsea in R, running PSEA, outputing files for making KEGG map diagrams, etc.

Load Packages

This notebook makes use of python via anaconda. Please make sure you install anaconda to proceed:

conda_list()

Select your conda environment of choice, in my case, item 2 from list above:

conda_list()[[2]][2] %>% 
  use_condaenv(required = TRUE)
print(paste("Using cond env:",conda_list()[[1]][2] ))
[1] "Using cond env: base"

Setup

Select the CYCLOPS ordering path and path to data file and path to DEseq2 differential expression results:

#path to cyclops ordering
path_to_cyclops_ordering = "/Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit/"

#path to normalized counts
path_to_tmm_file = "/Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/tmms/scROSMAP/cogdx_controls/EdgeR_filt_normed/ExcSubtypes3and5_FiltByExprDefault_TMM.csv"

# optional path to file containing differential expression results
deseq_de_filename = "/Users/henryhollis/Box Sync/Henry_stuff/AD_project/scROSMAP/simple_differential_expr/Exc3_5/Exc3_5_filtByExpr_combatSeq_edgeR_DE_exactTest.csv"

path_to_rosmap_clin = "~/Box Sync/Henry_stuff/AD_project/scROSMAP/Meta_data/cleaned_rosmap_meta_cogdxConds.csv"
metabolon_filename = "~/Box Sync/Henry_stuff/AD_project/scROSMAP/Metabolon/ROSMAP_Metabolon_HD4_Brain514_assay_data.csv"
metabolon_datakey = "~/Box Sync/Henry_stuff/AD_project/scROSMAP/Metabolon/ROSMAP_Metabolon_HD4_Data_Dictionary.csv"

#BHQ corrected p cutoff to be considered cycling in analysis:
isCyclingBHQCutoff = 0.1
percentile = 0.025
useBatch = T

Cosinor Regression and Differential Rhythmicity

Run cosinor regression and differential rhythmicity analysis with custom script and with compareRhythms library. I make use of Ensembl and Entrez IDs so I need the following translation dictionaries:

# read in translation dict between gene symbols and other alias
Ensembl_dict = readxl::read_xlsx("scROSMAP_ENSEMBL_dict.xlsx")

#run custom cycling and differential rhythmicity analysis
print("***Running cosinor and differential rhythms analysis***")
[1] "***Running cosinor and differential rhythms analysis***"
run_cycling_and_dr_analysis(path_to_cyclops_ordering, path_to_tmm_file, isCyclingSigCutoff = isCyclingBHQCutoff, percentile = percentile, useBatch = useBatch)

Running is_cycling() on cond_subset: cond_0
NOTE: Using batches in regression.

[>-------------------------------------------------------------------------------------------------------------------------------------------]   1%
[=>------------------------------------------------------------------------------------------------------------------------------------------]   1%
[=>------------------------------------------------------------------------------------------------------------------------------------------]   2%
[==>-----------------------------------------------------------------------------------------------------------------------------------------]   2%
[===>----------------------------------------------------------------------------------------------------------------------------------------]   3%
[====>---------------------------------------------------------------------------------------------------------------------------------------]   3%
[====>---------------------------------------------------------------------------------------------------------------------------------------]   4%
[=====>--------------------------------------------------------------------------------------------------------------------------------------]   4%
[=====>--------------------------------------------------------------------------------------------------------------------------------------]   5%
[======>-------------------------------------------------------------------------------------------------------------------------------------]   5%
[=======>------------------------------------------------------------------------------------------------------------------------------------]   5%
[=======>------------------------------------------------------------------------------------------------------------------------------------]   6%
[========>-----------------------------------------------------------------------------------------------------------------------------------]   6%
[========>-----------------------------------------------------------------------------------------------------------------------------------]   7%
[=========>----------------------------------------------------------------------------------------------------------------------------------]   7%
[==========>---------------------------------------------------------------------------------------------------------------------------------]   8%
[===========>--------------------------------------------------------------------------------------------------------------------------------]   8%
[===========>--------------------------------------------------------------------------------------------------------------------------------]   9%
[============>-------------------------------------------------------------------------------------------------------------------------------]   9%
[============>-------------------------------------------------------------------------------------------------------------------------------]  10%
[=============>------------------------------------------------------------------------------------------------------------------------------]  10%
[==============>-----------------------------------------------------------------------------------------------------------------------------]  10%
[==============>-----------------------------------------------------------------------------------------------------------------------------]  11%
[===============>----------------------------------------------------------------------------------------------------------------------------]  11%
[===============>----------------------------------------------------------------------------------------------------------------------------]  12%
[================>---------------------------------------------------------------------------------------------------------------------------]  12%
[=================>--------------------------------------------------------------------------------------------------------------------------]  13%
[==================>-------------------------------------------------------------------------------------------------------------------------]  13%
[==================>-------------------------------------------------------------------------------------------------------------------------]  14%
[===================>------------------------------------------------------------------------------------------------------------------------]  14%
[===================>------------------------------------------------------------------------------------------------------------------------]  15%
[====================>-----------------------------------------------------------------------------------------------------------------------]  15%
[=====================>----------------------------------------------------------------------------------------------------------------------]  15%
[=====================>----------------------------------------------------------------------------------------------------------------------]  16%
[======================>---------------------------------------------------------------------------------------------------------------------]  16%
[======================>---------------------------------------------------------------------------------------------------------------------]  17%
[=======================>--------------------------------------------------------------------------------------------------------------------]  17%
[========================>-------------------------------------------------------------------------------------------------------------------]  18%
[=========================>------------------------------------------------------------------------------------------------------------------]  18%
[=========================>------------------------------------------------------------------------------------------------------------------]  19%
[==========================>-----------------------------------------------------------------------------------------------------------------]  19%
[==========================>-----------------------------------------------------------------------------------------------------------------]  20%
[===========================>----------------------------------------------------------------------------------------------------------------]  20%
[============================>---------------------------------------------------------------------------------------------------------------]  20%
[============================>---------------------------------------------------------------------------------------------------------------]  21%
[=============================>--------------------------------------------------------------------------------------------------------------]  21%
[=============================>--------------------------------------------------------------------------------------------------------------]  22%
[==============================>-------------------------------------------------------------------------------------------------------------]  22%
[===============================>------------------------------------------------------------------------------------------------------------]  23%
[================================>-----------------------------------------------------------------------------------------------------------]  23%
[================================>-----------------------------------------------------------------------------------------------------------]  24%
[=================================>----------------------------------------------------------------------------------------------------------]  24%
[=================================>----------------------------------------------------------------------------------------------------------]  25%
[==================================>---------------------------------------------------------------------------------------------------------]  25%
[===================================>--------------------------------------------------------------------------------------------------------]  25%
[===================================>--------------------------------------------------------------------------------------------------------]  26%
[====================================>-------------------------------------------------------------------------------------------------------]  26%
[====================================>-------------------------------------------------------------------------------------------------------]  27%
[=====================================>------------------------------------------------------------------------------------------------------]  27%
[======================================>-----------------------------------------------------------------------------------------------------]  28%
[=======================================>----------------------------------------------------------------------------------------------------]  28%
[=======================================>----------------------------------------------------------------------------------------------------]  29%
[========================================>---------------------------------------------------------------------------------------------------]  29%
[========================================>---------------------------------------------------------------------------------------------------]  30%
[=========================================>--------------------------------------------------------------------------------------------------]  30%
[==========================================>-------------------------------------------------------------------------------------------------]  30%
[==========================================>-------------------------------------------------------------------------------------------------]  31%
[===========================================>------------------------------------------------------------------------------------------------]  31%
[===========================================>------------------------------------------------------------------------------------------------]  32%
[============================================>-----------------------------------------------------------------------------------------------]  32%
[=============================================>----------------------------------------------------------------------------------------------]  33%
[==============================================>---------------------------------------------------------------------------------------------]  33%
[==============================================>---------------------------------------------------------------------------------------------]  34%
[===============================================>--------------------------------------------------------------------------------------------]  34%
[===============================================>--------------------------------------------------------------------------------------------]  35%
[================================================>-------------------------------------------------------------------------------------------]  35%
[=================================================>------------------------------------------------------------------------------------------]  35%
[=================================================>------------------------------------------------------------------------------------------]  36%
[==================================================>-----------------------------------------------------------------------------------------]  36%
[==================================================>-----------------------------------------------------------------------------------------]  37%
[===================================================>----------------------------------------------------------------------------------------]  37%
[====================================================>---------------------------------------------------------------------------------------]  38%
[=====================================================>--------------------------------------------------------------------------------------]  38%
[=====================================================>--------------------------------------------------------------------------------------]  39%
[======================================================>-------------------------------------------------------------------------------------]  39%
[======================================================>-------------------------------------------------------------------------------------]  40%
[=======================================================>------------------------------------------------------------------------------------]  40%
[========================================================>-----------------------------------------------------------------------------------]  40%
[========================================================>-----------------------------------------------------------------------------------]  41%
[=========================================================>----------------------------------------------------------------------------------]  41%
[=========================================================>----------------------------------------------------------------------------------]  42%
[==========================================================>---------------------------------------------------------------------------------]  42%
[===========================================================>--------------------------------------------------------------------------------]  43%
[============================================================>-------------------------------------------------------------------------------]  43%
[============================================================>-------------------------------------------------------------------------------]  44%
[=============================================================>------------------------------------------------------------------------------]  44%
[=============================================================>------------------------------------------------------------------------------]  45%
[==============================================================>-----------------------------------------------------------------------------]  45%
[===============================================================>----------------------------------------------------------------------------]  45%
[===============================================================>----------------------------------------------------------------------------]  46%
[================================================================>---------------------------------------------------------------------------]  46%
[================================================================>---------------------------------------------------------------------------]  47%
[=================================================================>--------------------------------------------------------------------------]  47%
[==================================================================>-------------------------------------------------------------------------]  48%
[===================================================================>------------------------------------------------------------------------]  48%
[===================================================================>------------------------------------------------------------------------]  49%
[====================================================================>-----------------------------------------------------------------------]  49%
[====================================================================>-----------------------------------------------------------------------]  50%
[=====================================================================>----------------------------------------------------------------------]  50%
[======================================================================>---------------------------------------------------------------------]  50%
[======================================================================>---------------------------------------------------------------------]  51%
[=======================================================================>--------------------------------------------------------------------]  51%
[=======================================================================>--------------------------------------------------------------------]  52%
[========================================================================>-------------------------------------------------------------------]  52%
[=========================================================================>------------------------------------------------------------------]  53%
[==========================================================================>-----------------------------------------------------------------]  53%
[==========================================================================>-----------------------------------------------------------------]  54%
[===========================================================================>----------------------------------------------------------------]  54%
[===========================================================================>----------------------------------------------------------------]  55%
[============================================================================>---------------------------------------------------------------]  55%
[=============================================================================>--------------------------------------------------------------]  55%
[=============================================================================>--------------------------------------------------------------]  56%
[==============================================================================>-------------------------------------------------------------]  56%
[==============================================================================>-------------------------------------------------------------]  57%
[===============================================================================>------------------------------------------------------------]  57%
[================================================================================>-----------------------------------------------------------]  58%
[=================================================================================>----------------------------------------------------------]  58%
[=================================================================================>----------------------------------------------------------]  59%
[==================================================================================>---------------------------------------------------------]  59%
[==================================================================================>---------------------------------------------------------]  60%
[===================================================================================>--------------------------------------------------------]  60%
[====================================================================================>-------------------------------------------------------]  60%
[====================================================================================>-------------------------------------------------------]  61%
[=====================================================================================>------------------------------------------------------]  61%
[=====================================================================================>------------------------------------------------------]  62%
[======================================================================================>-----------------------------------------------------]  62%
[=======================================================================================>----------------------------------------------------]  63%
[========================================================================================>---------------------------------------------------]  63%
[========================================================================================>---------------------------------------------------]  64%
[=========================================================================================>--------------------------------------------------]  64%
[=========================================================================================>--------------------------------------------------]  65%
[==========================================================================================>-------------------------------------------------]  65%
[===========================================================================================>------------------------------------------------]  65%
[===========================================================================================>------------------------------------------------]  66%
[============================================================================================>-----------------------------------------------]  66%
[============================================================================================>-----------------------------------------------]  67%
[=============================================================================================>----------------------------------------------]  67%
[==============================================================================================>---------------------------------------------]  68%
[===============================================================================================>--------------------------------------------]  68%
[===============================================================================================>--------------------------------------------]  69%
[================================================================================================>-------------------------------------------]  69%
[================================================================================================>-------------------------------------------]  70%
[=================================================================================================>------------------------------------------]  70%
[==================================================================================================>-----------------------------------------]  70%
[==================================================================================================>-----------------------------------------]  71%
[===================================================================================================>----------------------------------------]  71%
[===================================================================================================>----------------------------------------]  72%
[====================================================================================================>---------------------------------------]  72%
[=====================================================================================================>--------------------------------------]  73%
[======================================================================================================>-------------------------------------]  73%
[======================================================================================================>-------------------------------------]  74%
[=======================================================================================================>------------------------------------]  74%
[=======================================================================================================>------------------------------------]  75%
[========================================================================================================>-----------------------------------]  75%
[=========================================================================================================>----------------------------------]  75%
[=========================================================================================================>----------------------------------]  76%
[==========================================================================================================>---------------------------------]  76%
[==========================================================================================================>---------------------------------]  77%
[===========================================================================================================>--------------------------------]  77%
[============================================================================================================>-------------------------------]  78%
[=============================================================================================================>------------------------------]  78%
[=============================================================================================================>------------------------------]  79%
[==============================================================================================================>-----------------------------]  79%
[==============================================================================================================>-----------------------------]  80%
[===============================================================================================================>----------------------------]  80%
[================================================================================================================>---------------------------]  80%
[================================================================================================================>---------------------------]  81%
[=================================================================================================================>--------------------------]  81%
[=================================================================================================================>--------------------------]  82%
[==================================================================================================================>-------------------------]  82%
[===================================================================================================================>------------------------]  83%
[====================================================================================================================>-----------------------]  83%
[====================================================================================================================>-----------------------]  84%
[=====================================================================================================================>----------------------]  84%
[=====================================================================================================================>----------------------]  85%
[======================================================================================================================>---------------------]  85%
[=======================================================================================================================>--------------------]  85%
[=======================================================================================================================>--------------------]  86%
[========================================================================================================================>-------------------]  86%
[========================================================================================================================>-------------------]  87%
[=========================================================================================================================>------------------]  87%
[==========================================================================================================================>-----------------]  88%
[===========================================================================================================================>----------------]  88%
[===========================================================================================================================>----------------]  89%
[============================================================================================================================>---------------]  89%
[============================================================================================================================>---------------]  90%
[=============================================================================================================================>--------------]  90%
[==============================================================================================================================>-------------]  90%
[==============================================================================================================================>-------------]  91%
[===============================================================================================================================>------------]  91%
[===============================================================================================================================>------------]  92%
[================================================================================================================================>-----------]  92%
[=================================================================================================================================>----------]  93%
[==================================================================================================================================>---------]  93%
[==================================================================================================================================>---------]  94%
[===================================================================================================================================>--------]  94%
[===================================================================================================================================>--------]  95%
[====================================================================================================================================>-------]  95%
[=====================================================================================================================================>------]  95%
[=====================================================================================================================================>------]  96%
[======================================================================================================================================>-----]  96%
[======================================================================================================================================>-----]  97%
[=======================================================================================================================================>----]  97%
[========================================================================================================================================>---]  98%
[=========================================================================================================================================>--]  98%
[=========================================================================================================================================>--]  99%
[==========================================================================================================================================>-]  99%
[==========================================================================================================================================>-] 100%
[===========================================================================================================================================>] 100%

Running is_cycling() on cond_subset: cond_1
NOTE: Using batches in regression.

[>-------------------------------------------------------------------------------------------------------------------------------------------]   1%
[=>------------------------------------------------------------------------------------------------------------------------------------------]   1%
[=>------------------------------------------------------------------------------------------------------------------------------------------]   2%
[==>-----------------------------------------------------------------------------------------------------------------------------------------]   2%
[===>----------------------------------------------------------------------------------------------------------------------------------------]   3%
[====>---------------------------------------------------------------------------------------------------------------------------------------]   3%
[====>---------------------------------------------------------------------------------------------------------------------------------------]   4%
[=====>--------------------------------------------------------------------------------------------------------------------------------------]   4%
[=====>--------------------------------------------------------------------------------------------------------------------------------------]   5%
[======>-------------------------------------------------------------------------------------------------------------------------------------]   5%
[=======>------------------------------------------------------------------------------------------------------------------------------------]   5%
[=======>------------------------------------------------------------------------------------------------------------------------------------]   6%
[========>-----------------------------------------------------------------------------------------------------------------------------------]   6%
[========>-----------------------------------------------------------------------------------------------------------------------------------]   7%
[=========>----------------------------------------------------------------------------------------------------------------------------------]   7%
[==========>---------------------------------------------------------------------------------------------------------------------------------]   8%
[===========>--------------------------------------------------------------------------------------------------------------------------------]   8%
[===========>--------------------------------------------------------------------------------------------------------------------------------]   9%
[============>-------------------------------------------------------------------------------------------------------------------------------]   9%
[============>-------------------------------------------------------------------------------------------------------------------------------]  10%
[=============>------------------------------------------------------------------------------------------------------------------------------]  10%
[==============>-----------------------------------------------------------------------------------------------------------------------------]  10%
[==============>-----------------------------------------------------------------------------------------------------------------------------]  11%
[===============>----------------------------------------------------------------------------------------------------------------------------]  11%
[===============>----------------------------------------------------------------------------------------------------------------------------]  12%
[================>---------------------------------------------------------------------------------------------------------------------------]  12%
[=================>--------------------------------------------------------------------------------------------------------------------------]  13%
[==================>-------------------------------------------------------------------------------------------------------------------------]  13%
[==================>-------------------------------------------------------------------------------------------------------------------------]  14%
[===================>------------------------------------------------------------------------------------------------------------------------]  14%
[===================>------------------------------------------------------------------------------------------------------------------------]  15%
[====================>-----------------------------------------------------------------------------------------------------------------------]  15%
[=====================>----------------------------------------------------------------------------------------------------------------------]  15%
[=====================>----------------------------------------------------------------------------------------------------------------------]  16%
[======================>---------------------------------------------------------------------------------------------------------------------]  16%
[======================>---------------------------------------------------------------------------------------------------------------------]  17%
[=======================>--------------------------------------------------------------------------------------------------------------------]  17%
[========================>-------------------------------------------------------------------------------------------------------------------]  18%
[=========================>------------------------------------------------------------------------------------------------------------------]  18%
[=========================>------------------------------------------------------------------------------------------------------------------]  19%
[==========================>-----------------------------------------------------------------------------------------------------------------]  19%
[==========================>-----------------------------------------------------------------------------------------------------------------]  20%
[===========================>----------------------------------------------------------------------------------------------------------------]  20%
[============================>---------------------------------------------------------------------------------------------------------------]  20%
[============================>---------------------------------------------------------------------------------------------------------------]  21%
[=============================>--------------------------------------------------------------------------------------------------------------]  21%
[=============================>--------------------------------------------------------------------------------------------------------------]  22%
[==============================>-------------------------------------------------------------------------------------------------------------]  22%
[===============================>------------------------------------------------------------------------------------------------------------]  23%
[================================>-----------------------------------------------------------------------------------------------------------]  23%
[================================>-----------------------------------------------------------------------------------------------------------]  24%
[=================================>----------------------------------------------------------------------------------------------------------]  24%
[=================================>----------------------------------------------------------------------------------------------------------]  25%
[==================================>---------------------------------------------------------------------------------------------------------]  25%
[===================================>--------------------------------------------------------------------------------------------------------]  25%
[===================================>--------------------------------------------------------------------------------------------------------]  26%
[====================================>-------------------------------------------------------------------------------------------------------]  26%
[====================================>-------------------------------------------------------------------------------------------------------]  27%
[=====================================>------------------------------------------------------------------------------------------------------]  27%
[======================================>-----------------------------------------------------------------------------------------------------]  28%
[=======================================>----------------------------------------------------------------------------------------------------]  28%
[=======================================>----------------------------------------------------------------------------------------------------]  29%
[========================================>---------------------------------------------------------------------------------------------------]  29%
[========================================>---------------------------------------------------------------------------------------------------]  30%
[=========================================>--------------------------------------------------------------------------------------------------]  30%
[==========================================>-------------------------------------------------------------------------------------------------]  30%
[==========================================>-------------------------------------------------------------------------------------------------]  31%
[===========================================>------------------------------------------------------------------------------------------------]  31%
[===========================================>------------------------------------------------------------------------------------------------]  32%
[============================================>-----------------------------------------------------------------------------------------------]  32%
[=============================================>----------------------------------------------------------------------------------------------]  33%
[==============================================>---------------------------------------------------------------------------------------------]  33%
[==============================================>---------------------------------------------------------------------------------------------]  34%
[===============================================>--------------------------------------------------------------------------------------------]  34%
[===============================================>--------------------------------------------------------------------------------------------]  35%
[================================================>-------------------------------------------------------------------------------------------]  35%
[=================================================>------------------------------------------------------------------------------------------]  35%
[=================================================>------------------------------------------------------------------------------------------]  36%
[==================================================>-----------------------------------------------------------------------------------------]  36%
[==================================================>-----------------------------------------------------------------------------------------]  37%
[===================================================>----------------------------------------------------------------------------------------]  37%
[====================================================>---------------------------------------------------------------------------------------]  38%
[=====================================================>--------------------------------------------------------------------------------------]  38%
[=====================================================>--------------------------------------------------------------------------------------]  39%
[======================================================>-------------------------------------------------------------------------------------]  39%
[======================================================>-------------------------------------------------------------------------------------]  40%
[=======================================================>------------------------------------------------------------------------------------]  40%
[========================================================>-----------------------------------------------------------------------------------]  40%
[========================================================>-----------------------------------------------------------------------------------]  41%
[=========================================================>----------------------------------------------------------------------------------]  41%
[=========================================================>----------------------------------------------------------------------------------]  42%
[==========================================================>---------------------------------------------------------------------------------]  42%
[===========================================================>--------------------------------------------------------------------------------]  43%
[============================================================>-------------------------------------------------------------------------------]  43%
[============================================================>-------------------------------------------------------------------------------]  44%
[=============================================================>------------------------------------------------------------------------------]  44%
[=============================================================>------------------------------------------------------------------------------]  45%
[==============================================================>-----------------------------------------------------------------------------]  45%
[===============================================================>----------------------------------------------------------------------------]  45%
[===============================================================>----------------------------------------------------------------------------]  46%
[================================================================>---------------------------------------------------------------------------]  46%
[================================================================>---------------------------------------------------------------------------]  47%
[=================================================================>--------------------------------------------------------------------------]  47%
[==================================================================>-------------------------------------------------------------------------]  48%
[===================================================================>------------------------------------------------------------------------]  48%
[===================================================================>------------------------------------------------------------------------]  49%
[====================================================================>-----------------------------------------------------------------------]  49%
[====================================================================>-----------------------------------------------------------------------]  50%
[=====================================================================>----------------------------------------------------------------------]  50%
[======================================================================>---------------------------------------------------------------------]  50%
[======================================================================>---------------------------------------------------------------------]  51%
[=======================================================================>--------------------------------------------------------------------]  51%
[=======================================================================>--------------------------------------------------------------------]  52%
[========================================================================>-------------------------------------------------------------------]  52%
[=========================================================================>------------------------------------------------------------------]  53%
[==========================================================================>-----------------------------------------------------------------]  53%
[==========================================================================>-----------------------------------------------------------------]  54%
[===========================================================================>----------------------------------------------------------------]  54%
[===========================================================================>----------------------------------------------------------------]  55%
[============================================================================>---------------------------------------------------------------]  55%
[=============================================================================>--------------------------------------------------------------]  55%
[=============================================================================>--------------------------------------------------------------]  56%
[==============================================================================>-------------------------------------------------------------]  56%
[==============================================================================>-------------------------------------------------------------]  57%
[===============================================================================>------------------------------------------------------------]  57%
[================================================================================>-----------------------------------------------------------]  58%
[=================================================================================>----------------------------------------------------------]  58%
[=================================================================================>----------------------------------------------------------]  59%
[==================================================================================>---------------------------------------------------------]  59%
[==================================================================================>---------------------------------------------------------]  60%
[===================================================================================>--------------------------------------------------------]  60%
[====================================================================================>-------------------------------------------------------]  60%
[====================================================================================>-------------------------------------------------------]  61%
[=====================================================================================>------------------------------------------------------]  61%
[=====================================================================================>------------------------------------------------------]  62%
[======================================================================================>-----------------------------------------------------]  62%
[=======================================================================================>----------------------------------------------------]  63%
[========================================================================================>---------------------------------------------------]  63%
[========================================================================================>---------------------------------------------------]  64%
[=========================================================================================>--------------------------------------------------]  64%
[=========================================================================================>--------------------------------------------------]  65%
[==========================================================================================>-------------------------------------------------]  65%
[===========================================================================================>------------------------------------------------]  65%
[===========================================================================================>------------------------------------------------]  66%
[============================================================================================>-----------------------------------------------]  66%
[============================================================================================>-----------------------------------------------]  67%
[=============================================================================================>----------------------------------------------]  67%
[==============================================================================================>---------------------------------------------]  68%
[===============================================================================================>--------------------------------------------]  68%
[===============================================================================================>--------------------------------------------]  69%
[================================================================================================>-------------------------------------------]  69%
[================================================================================================>-------------------------------------------]  70%
[=================================================================================================>------------------------------------------]  70%
[==================================================================================================>-----------------------------------------]  70%
[==================================================================================================>-----------------------------------------]  71%
[===================================================================================================>----------------------------------------]  71%
[===================================================================================================>----------------------------------------]  72%
[====================================================================================================>---------------------------------------]  72%
[=====================================================================================================>--------------------------------------]  73%
[======================================================================================================>-------------------------------------]  73%
[======================================================================================================>-------------------------------------]  74%
[=======================================================================================================>------------------------------------]  74%
[=======================================================================================================>------------------------------------]  75%
[========================================================================================================>-----------------------------------]  75%
[=========================================================================================================>----------------------------------]  75%
[=========================================================================================================>----------------------------------]  76%
[==========================================================================================================>---------------------------------]  76%
[==========================================================================================================>---------------------------------]  77%
[===========================================================================================================>--------------------------------]  77%
[============================================================================================================>-------------------------------]  78%
[=============================================================================================================>------------------------------]  78%
[=============================================================================================================>------------------------------]  79%
[==============================================================================================================>-----------------------------]  79%
[==============================================================================================================>-----------------------------]  80%
[===============================================================================================================>----------------------------]  80%
[================================================================================================================>---------------------------]  80%
[================================================================================================================>---------------------------]  81%
[=================================================================================================================>--------------------------]  81%
[=================================================================================================================>--------------------------]  82%
[==================================================================================================================>-------------------------]  82%
[===================================================================================================================>------------------------]  83%
[====================================================================================================================>-----------------------]  83%
[====================================================================================================================>-----------------------]  84%
[=====================================================================================================================>----------------------]  84%
[=====================================================================================================================>----------------------]  85%
[======================================================================================================================>---------------------]  85%
[=======================================================================================================================>--------------------]  85%
[=======================================================================================================================>--------------------]  86%
[========================================================================================================================>-------------------]  86%
[========================================================================================================================>-------------------]  87%
[=========================================================================================================================>------------------]  87%
[==========================================================================================================================>-----------------]  88%
[===========================================================================================================================>----------------]  88%
[===========================================================================================================================>----------------]  89%
[============================================================================================================================>---------------]  89%
[============================================================================================================================>---------------]  90%
[=============================================================================================================================>--------------]  90%
[==============================================================================================================================>-------------]  90%
[==============================================================================================================================>-------------]  91%
[===============================================================================================================================>------------]  91%
[===============================================================================================================================>------------]  92%
[================================================================================================================================>-----------]  92%
[=================================================================================================================================>----------]  93%
[==================================================================================================================================>---------]  93%
[==================================================================================================================================>---------]  94%
[===================================================================================================================================>--------]  94%
[===================================================================================================================================>--------]  95%
[====================================================================================================================================>-------]  95%
[=====================================================================================================================================>------]  95%
[=====================================================================================================================================>------]  96%
[======================================================================================================================================>-----]  96%
[======================================================================================================================================>-----]  97%
[=======================================================================================================================================>----]  97%
[========================================================================================================================================>---]  98%
[=========================================================================================================================================>--]  98%
[=========================================================================================================================================>--]  99%
[==========================================================================================================================================>-]  99%
[==========================================================================================================================================>-] 100%
[===========================================================================================================================================>] 100%

Running is_cycling_method2() (analogue of compareRhythms) on all subjects
NOTE: Using batches in regression.

[--------------------------------------------------------------------------------------------------------------------------------------------]   0%
[>-------------------------------------------------------------------------------------------------------------------------------------------]   0%
[>-------------------------------------------------------------------------------------------------------------------------------------------]   1%
[=>------------------------------------------------------------------------------------------------------------------------------------------]   1%
[=>------------------------------------------------------------------------------------------------------------------------------------------]   2%
[==>-----------------------------------------------------------------------------------------------------------------------------------------]   2%
[===>----------------------------------------------------------------------------------------------------------------------------------------]   3%
[====>---------------------------------------------------------------------------------------------------------------------------------------]   3%
[====>---------------------------------------------------------------------------------------------------------------------------------------]   4%
[=====>--------------------------------------------------------------------------------------------------------------------------------------]   4%
[=====>--------------------------------------------------------------------------------------------------------------------------------------]   5%
[======>-------------------------------------------------------------------------------------------------------------------------------------]   5%
[=======>------------------------------------------------------------------------------------------------------------------------------------]   5%
[=======>------------------------------------------------------------------------------------------------------------------------------------]   6%
[========>-----------------------------------------------------------------------------------------------------------------------------------]   6%
[========>-----------------------------------------------------------------------------------------------------------------------------------]   7%
[=========>----------------------------------------------------------------------------------------------------------------------------------]   7%
[==========>---------------------------------------------------------------------------------------------------------------------------------]   8%
[===========>--------------------------------------------------------------------------------------------------------------------------------]   8%
[===========>--------------------------------------------------------------------------------------------------------------------------------]   9%
[============>-------------------------------------------------------------------------------------------------------------------------------]   9%
[============>-------------------------------------------------------------------------------------------------------------------------------]  10%
[=============>------------------------------------------------------------------------------------------------------------------------------]  10%
[==============>-----------------------------------------------------------------------------------------------------------------------------]  10%
[==============>-----------------------------------------------------------------------------------------------------------------------------]  11%
[===============>----------------------------------------------------------------------------------------------------------------------------]  11%
[===============>----------------------------------------------------------------------------------------------------------------------------]  12%
[================>---------------------------------------------------------------------------------------------------------------------------]  12%
[=================>--------------------------------------------------------------------------------------------------------------------------]  13%
[==================>-------------------------------------------------------------------------------------------------------------------------]  13%
[==================>-------------------------------------------------------------------------------------------------------------------------]  14%
[===================>------------------------------------------------------------------------------------------------------------------------]  14%
[===================>------------------------------------------------------------------------------------------------------------------------]  15%
[====================>-----------------------------------------------------------------------------------------------------------------------]  15%
[=====================>----------------------------------------------------------------------------------------------------------------------]  15%
[=====================>----------------------------------------------------------------------------------------------------------------------]  16%
[======================>---------------------------------------------------------------------------------------------------------------------]  16%
[======================>---------------------------------------------------------------------------------------------------------------------]  17%
[=======================>--------------------------------------------------------------------------------------------------------------------]  17%
[========================>-------------------------------------------------------------------------------------------------------------------]  18%
[=========================>------------------------------------------------------------------------------------------------------------------]  18%
[=========================>------------------------------------------------------------------------------------------------------------------]  19%
[==========================>-----------------------------------------------------------------------------------------------------------------]  19%
[==========================>-----------------------------------------------------------------------------------------------------------------]  20%
[===========================>----------------------------------------------------------------------------------------------------------------]  20%
[============================>---------------------------------------------------------------------------------------------------------------]  20%
[============================>---------------------------------------------------------------------------------------------------------------]  21%
[=============================>--------------------------------------------------------------------------------------------------------------]  21%
[=============================>--------------------------------------------------------------------------------------------------------------]  22%
[==============================>-------------------------------------------------------------------------------------------------------------]  22%
[===============================>------------------------------------------------------------------------------------------------------------]  23%
[================================>-----------------------------------------------------------------------------------------------------------]  23%
[================================>-----------------------------------------------------------------------------------------------------------]  24%
[=================================>----------------------------------------------------------------------------------------------------------]  24%
[=================================>----------------------------------------------------------------------------------------------------------]  25%
[==================================>---------------------------------------------------------------------------------------------------------]  25%
[===================================>--------------------------------------------------------------------------------------------------------]  25%
[===================================>--------------------------------------------------------------------------------------------------------]  26%
[====================================>-------------------------------------------------------------------------------------------------------]  26%
[====================================>-------------------------------------------------------------------------------------------------------]  27%
[=====================================>------------------------------------------------------------------------------------------------------]  27%
[======================================>-----------------------------------------------------------------------------------------------------]  28%
[=======================================>----------------------------------------------------------------------------------------------------]  28%
[=======================================>----------------------------------------------------------------------------------------------------]  29%
[========================================>---------------------------------------------------------------------------------------------------]  29%
[========================================>---------------------------------------------------------------------------------------------------]  30%
[=========================================>--------------------------------------------------------------------------------------------------]  30%
[==========================================>-------------------------------------------------------------------------------------------------]  30%
[==========================================>-------------------------------------------------------------------------------------------------]  31%
[===========================================>------------------------------------------------------------------------------------------------]  31%
[===========================================>------------------------------------------------------------------------------------------------]  32%
[============================================>-----------------------------------------------------------------------------------------------]  32%
[=============================================>----------------------------------------------------------------------------------------------]  33%
[==============================================>---------------------------------------------------------------------------------------------]  33%
[==============================================>---------------------------------------------------------------------------------------------]  34%
[===============================================>--------------------------------------------------------------------------------------------]  34%
[===============================================>--------------------------------------------------------------------------------------------]  35%
[================================================>-------------------------------------------------------------------------------------------]  35%
[=================================================>------------------------------------------------------------------------------------------]  35%
[=================================================>------------------------------------------------------------------------------------------]  36%
[==================================================>-----------------------------------------------------------------------------------------]  36%
[==================================================>-----------------------------------------------------------------------------------------]  37%
[===================================================>----------------------------------------------------------------------------------------]  37%
[====================================================>---------------------------------------------------------------------------------------]  38%
[=====================================================>--------------------------------------------------------------------------------------]  38%
[=====================================================>--------------------------------------------------------------------------------------]  39%
[======================================================>-------------------------------------------------------------------------------------]  39%
[======================================================>-------------------------------------------------------------------------------------]  40%
[=======================================================>------------------------------------------------------------------------------------]  40%
[========================================================>-----------------------------------------------------------------------------------]  40%
[========================================================>-----------------------------------------------------------------------------------]  41%
[=========================================================>----------------------------------------------------------------------------------]  41%
[=========================================================>----------------------------------------------------------------------------------]  42%
[==========================================================>---------------------------------------------------------------------------------]  42%
[===========================================================>--------------------------------------------------------------------------------]  43%
[============================================================>-------------------------------------------------------------------------------]  43%
[============================================================>-------------------------------------------------------------------------------]  44%
[=============================================================>------------------------------------------------------------------------------]  44%
[=============================================================>------------------------------------------------------------------------------]  45%
[==============================================================>-----------------------------------------------------------------------------]  45%
[===============================================================>----------------------------------------------------------------------------]  45%
[===============================================================>----------------------------------------------------------------------------]  46%
[================================================================>---------------------------------------------------------------------------]  46%
[================================================================>---------------------------------------------------------------------------]  47%
[=================================================================>--------------------------------------------------------------------------]  47%
[==================================================================>-------------------------------------------------------------------------]  48%
[===================================================================>------------------------------------------------------------------------]  48%
[===================================================================>------------------------------------------------------------------------]  49%
[====================================================================>-----------------------------------------------------------------------]  49%
[====================================================================>-----------------------------------------------------------------------]  50%
[=====================================================================>----------------------------------------------------------------------]  50%
[======================================================================>---------------------------------------------------------------------]  50%
[======================================================================>---------------------------------------------------------------------]  51%
[=======================================================================>--------------------------------------------------------------------]  51%
[=======================================================================>--------------------------------------------------------------------]  52%
[========================================================================>-------------------------------------------------------------------]  52%
[=========================================================================>------------------------------------------------------------------]  53%
[==========================================================================>-----------------------------------------------------------------]  53%
[==========================================================================>-----------------------------------------------------------------]  54%
[===========================================================================>----------------------------------------------------------------]  54%
[===========================================================================>----------------------------------------------------------------]  55%
[============================================================================>---------------------------------------------------------------]  55%
[=============================================================================>--------------------------------------------------------------]  55%
[=============================================================================>--------------------------------------------------------------]  56%
[==============================================================================>-------------------------------------------------------------]  56%
[==============================================================================>-------------------------------------------------------------]  57%
[===============================================================================>------------------------------------------------------------]  57%
[================================================================================>-----------------------------------------------------------]  58%
[=================================================================================>----------------------------------------------------------]  58%
[=================================================================================>----------------------------------------------------------]  59%
[==================================================================================>---------------------------------------------------------]  59%
[==================================================================================>---------------------------------------------------------]  60%
[===================================================================================>--------------------------------------------------------]  60%
[====================================================================================>-------------------------------------------------------]  60%
[====================================================================================>-------------------------------------------------------]  61%
[=====================================================================================>------------------------------------------------------]  61%
[=====================================================================================>------------------------------------------------------]  62%
[======================================================================================>-----------------------------------------------------]  62%
[=======================================================================================>----------------------------------------------------]  63%
[========================================================================================>---------------------------------------------------]  63%
[========================================================================================>---------------------------------------------------]  64%
[=========================================================================================>--------------------------------------------------]  64%
[=========================================================================================>--------------------------------------------------]  65%
[==========================================================================================>-------------------------------------------------]  65%
[===========================================================================================>------------------------------------------------]  65%
[===========================================================================================>------------------------------------------------]  66%
[============================================================================================>-----------------------------------------------]  66%
[============================================================================================>-----------------------------------------------]  67%
[=============================================================================================>----------------------------------------------]  67%
[==============================================================================================>---------------------------------------------]  68%
[===============================================================================================>--------------------------------------------]  68%
[===============================================================================================>--------------------------------------------]  69%
[================================================================================================>-------------------------------------------]  69%
[================================================================================================>-------------------------------------------]  70%
[=================================================================================================>------------------------------------------]  70%
[==================================================================================================>-----------------------------------------]  70%
[==================================================================================================>-----------------------------------------]  71%
[===================================================================================================>----------------------------------------]  71%
[===================================================================================================>----------------------------------------]  72%
[====================================================================================================>---------------------------------------]  72%
[=====================================================================================================>--------------------------------------]  73%
[======================================================================================================>-------------------------------------]  73%
[======================================================================================================>-------------------------------------]  74%
[=======================================================================================================>------------------------------------]  74%
[=======================================================================================================>------------------------------------]  75%
[========================================================================================================>-----------------------------------]  75%
[=========================================================================================================>----------------------------------]  75%
[=========================================================================================================>----------------------------------]  76%
[==========================================================================================================>---------------------------------]  76%
[==========================================================================================================>---------------------------------]  77%
[===========================================================================================================>--------------------------------]  77%
[============================================================================================================>-------------------------------]  78%
[=============================================================================================================>------------------------------]  78%
[=============================================================================================================>------------------------------]  79%
[==============================================================================================================>-----------------------------]  79%
[==============================================================================================================>-----------------------------]  80%
[===============================================================================================================>----------------------------]  80%
[================================================================================================================>---------------------------]  80%
[================================================================================================================>---------------------------]  81%
[=================================================================================================================>--------------------------]  81%
[=================================================================================================================>--------------------------]  82%
[==================================================================================================================>-------------------------]  82%
[===================================================================================================================>------------------------]  83%
[====================================================================================================================>-----------------------]  83%
[====================================================================================================================>-----------------------]  84%
[=====================================================================================================================>----------------------]  84%
[=====================================================================================================================>----------------------]  85%
[======================================================================================================================>---------------------]  85%
[=======================================================================================================================>--------------------]  85%
[=======================================================================================================================>--------------------]  86%
[========================================================================================================================>-------------------]  86%
[========================================================================================================================>-------------------]  87%
[=========================================================================================================================>------------------]  87%
[==========================================================================================================================>-----------------]  88%
[===========================================================================================================================>----------------]  88%
[===========================================================================================================================>----------------]  89%
[============================================================================================================================>---------------]  89%
[============================================================================================================================>---------------]  90%
[=============================================================================================================================>--------------]  90%
[==============================================================================================================================>-------------]  90%
[==============================================================================================================================>-------------]  91%
[===============================================================================================================================>------------]  91%
[===============================================================================================================================>------------]  92%
[================================================================================================================================>-----------]  92%
[=================================================================================================================================>----------]  93%
[==================================================================================================================================>---------]  93%
[==================================================================================================================================>---------]  94%
[===================================================================================================================================>--------]  94%
[===================================================================================================================================>--------]  95%
[====================================================================================================================================>-------]  95%
[=====================================================================================================================================>------]  95%
[=====================================================================================================================================>------]  96%
[======================================================================================================================================>-----]  96%
[======================================================================================================================================>-----]  97%
[=======================================================================================================================================>----]  97%
[========================================================================================================================================>---]  98%
[=========================================================================================================================================>--]  98%
[=========================================================================================================================================>--]  99%
[==========================================================================================================================================>-]  99%
[==========================================================================================================================================>-] 100%
[===========================================================================================================================================>] 100%

Running Mesor_differences()
NOTE: Using batches in regression.

[--------------------------------------------------------------------------------------------------------------------------------------------]   0%
[>-------------------------------------------------------------------------------------------------------------------------------------------]   0%
[>-------------------------------------------------------------------------------------------------------------------------------------------]   1%
[=>------------------------------------------------------------------------------------------------------------------------------------------]   1%
[=>------------------------------------------------------------------------------------------------------------------------------------------]   2%
[==>-----------------------------------------------------------------------------------------------------------------------------------------]   2%
[===>----------------------------------------------------------------------------------------------------------------------------------------]   3%
[====>---------------------------------------------------------------------------------------------------------------------------------------]   3%
[====>---------------------------------------------------------------------------------------------------------------------------------------]   4%
[=====>--------------------------------------------------------------------------------------------------------------------------------------]   4%
[=====>--------------------------------------------------------------------------------------------------------------------------------------]   5%
[======>-------------------------------------------------------------------------------------------------------------------------------------]   5%
[=======>------------------------------------------------------------------------------------------------------------------------------------]   5%
[=======>------------------------------------------------------------------------------------------------------------------------------------]   6%
[========>-----------------------------------------------------------------------------------------------------------------------------------]   6%
[========>-----------------------------------------------------------------------------------------------------------------------------------]   7%
[=========>----------------------------------------------------------------------------------------------------------------------------------]   7%
[==========>---------------------------------------------------------------------------------------------------------------------------------]   8%
[===========>--------------------------------------------------------------------------------------------------------------------------------]   8%
[===========>--------------------------------------------------------------------------------------------------------------------------------]   9%
[============>-------------------------------------------------------------------------------------------------------------------------------]   9%
[============>-------------------------------------------------------------------------------------------------------------------------------]  10%
[=============>------------------------------------------------------------------------------------------------------------------------------]  10%
[==============>-----------------------------------------------------------------------------------------------------------------------------]  10%
[==============>-----------------------------------------------------------------------------------------------------------------------------]  11%
[===============>----------------------------------------------------------------------------------------------------------------------------]  11%
[===============>----------------------------------------------------------------------------------------------------------------------------]  12%
[================>---------------------------------------------------------------------------------------------------------------------------]  12%
[=================>--------------------------------------------------------------------------------------------------------------------------]  13%
[==================>-------------------------------------------------------------------------------------------------------------------------]  13%
[==================>-------------------------------------------------------------------------------------------------------------------------]  14%
[===================>------------------------------------------------------------------------------------------------------------------------]  14%
[===================>------------------------------------------------------------------------------------------------------------------------]  15%
[====================>-----------------------------------------------------------------------------------------------------------------------]  15%
[=====================>----------------------------------------------------------------------------------------------------------------------]  15%
[=====================>----------------------------------------------------------------------------------------------------------------------]  16%
[======================>---------------------------------------------------------------------------------------------------------------------]  16%
[======================>---------------------------------------------------------------------------------------------------------------------]  17%
[=======================>--------------------------------------------------------------------------------------------------------------------]  17%
[========================>-------------------------------------------------------------------------------------------------------------------]  18%
[=========================>------------------------------------------------------------------------------------------------------------------]  18%
[=========================>------------------------------------------------------------------------------------------------------------------]  19%
[==========================>-----------------------------------------------------------------------------------------------------------------]  19%
[==========================>-----------------------------------------------------------------------------------------------------------------]  20%
[===========================>----------------------------------------------------------------------------------------------------------------]  20%
[============================>---------------------------------------------------------------------------------------------------------------]  20%
[============================>---------------------------------------------------------------------------------------------------------------]  21%
[=============================>--------------------------------------------------------------------------------------------------------------]  21%
[=============================>--------------------------------------------------------------------------------------------------------------]  22%
[==============================>-------------------------------------------------------------------------------------------------------------]  22%
[===============================>------------------------------------------------------------------------------------------------------------]  23%
[================================>-----------------------------------------------------------------------------------------------------------]  23%
[================================>-----------------------------------------------------------------------------------------------------------]  24%
[=================================>----------------------------------------------------------------------------------------------------------]  24%
[=================================>----------------------------------------------------------------------------------------------------------]  25%
[==================================>---------------------------------------------------------------------------------------------------------]  25%
[===================================>--------------------------------------------------------------------------------------------------------]  25%
[===================================>--------------------------------------------------------------------------------------------------------]  26%
[====================================>-------------------------------------------------------------------------------------------------------]  26%
[====================================>-------------------------------------------------------------------------------------------------------]  27%
[=====================================>------------------------------------------------------------------------------------------------------]  27%
[======================================>-----------------------------------------------------------------------------------------------------]  28%
[=======================================>----------------------------------------------------------------------------------------------------]  28%
[=======================================>----------------------------------------------------------------------------------------------------]  29%
[========================================>---------------------------------------------------------------------------------------------------]  29%
[========================================>---------------------------------------------------------------------------------------------------]  30%
[=========================================>--------------------------------------------------------------------------------------------------]  30%
[==========================================>-------------------------------------------------------------------------------------------------]  30%
[==========================================>-------------------------------------------------------------------------------------------------]  31%
[===========================================>------------------------------------------------------------------------------------------------]  31%
[===========================================>------------------------------------------------------------------------------------------------]  32%
[============================================>-----------------------------------------------------------------------------------------------]  32%
[=============================================>----------------------------------------------------------------------------------------------]  33%
[==============================================>---------------------------------------------------------------------------------------------]  33%
[==============================================>---------------------------------------------------------------------------------------------]  34%
[===============================================>--------------------------------------------------------------------------------------------]  34%
[===============================================>--------------------------------------------------------------------------------------------]  35%
[================================================>-------------------------------------------------------------------------------------------]  35%
[=================================================>------------------------------------------------------------------------------------------]  35%
[=================================================>------------------------------------------------------------------------------------------]  36%
[==================================================>-----------------------------------------------------------------------------------------]  36%
[==================================================>-----------------------------------------------------------------------------------------]  37%
[===================================================>----------------------------------------------------------------------------------------]  37%
[====================================================>---------------------------------------------------------------------------------------]  38%
[=====================================================>--------------------------------------------------------------------------------------]  38%
[=====================================================>--------------------------------------------------------------------------------------]  39%
[======================================================>-------------------------------------------------------------------------------------]  39%
[======================================================>-------------------------------------------------------------------------------------]  40%
[=======================================================>------------------------------------------------------------------------------------]  40%
[========================================================>-----------------------------------------------------------------------------------]  40%
[========================================================>-----------------------------------------------------------------------------------]  41%
[=========================================================>----------------------------------------------------------------------------------]  41%
[=========================================================>----------------------------------------------------------------------------------]  42%
[==========================================================>---------------------------------------------------------------------------------]  42%
[===========================================================>--------------------------------------------------------------------------------]  43%
[============================================================>-------------------------------------------------------------------------------]  43%
[============================================================>-------------------------------------------------------------------------------]  44%
[=============================================================>------------------------------------------------------------------------------]  44%
[=============================================================>------------------------------------------------------------------------------]  45%
[==============================================================>-----------------------------------------------------------------------------]  45%
[===============================================================>----------------------------------------------------------------------------]  45%
[===============================================================>----------------------------------------------------------------------------]  46%
[================================================================>---------------------------------------------------------------------------]  46%
[================================================================>---------------------------------------------------------------------------]  47%
[=================================================================>--------------------------------------------------------------------------]  47%
[==================================================================>-------------------------------------------------------------------------]  48%
[===================================================================>------------------------------------------------------------------------]  48%
[===================================================================>------------------------------------------------------------------------]  49%
[====================================================================>-----------------------------------------------------------------------]  49%
[====================================================================>-----------------------------------------------------------------------]  50%
[=====================================================================>----------------------------------------------------------------------]  50%
[======================================================================>---------------------------------------------------------------------]  50%
[======================================================================>---------------------------------------------------------------------]  51%
[=======================================================================>--------------------------------------------------------------------]  51%
[=======================================================================>--------------------------------------------------------------------]  52%
[========================================================================>-------------------------------------------------------------------]  52%
[=========================================================================>------------------------------------------------------------------]  53%
[==========================================================================>-----------------------------------------------------------------]  53%
[==========================================================================>-----------------------------------------------------------------]  54%
[===========================================================================>----------------------------------------------------------------]  54%
[===========================================================================>----------------------------------------------------------------]  55%
[============================================================================>---------------------------------------------------------------]  55%
[=============================================================================>--------------------------------------------------------------]  55%
[=============================================================================>--------------------------------------------------------------]  56%
[==============================================================================>-------------------------------------------------------------]  56%
[==============================================================================>-------------------------------------------------------------]  57%
[===============================================================================>------------------------------------------------------------]  57%
[================================================================================>-----------------------------------------------------------]  58%
[=================================================================================>----------------------------------------------------------]  58%
[=================================================================================>----------------------------------------------------------]  59%
[==================================================================================>---------------------------------------------------------]  59%
[==================================================================================>---------------------------------------------------------]  60%
[===================================================================================>--------------------------------------------------------]  60%
[====================================================================================>-------------------------------------------------------]  60%
[====================================================================================>-------------------------------------------------------]  61%
[=====================================================================================>------------------------------------------------------]  61%
[=====================================================================================>------------------------------------------------------]  62%
[======================================================================================>-----------------------------------------------------]  62%
[=======================================================================================>----------------------------------------------------]  63%
[========================================================================================>---------------------------------------------------]  63%
[========================================================================================>---------------------------------------------------]  64%
[=========================================================================================>--------------------------------------------------]  64%
[=========================================================================================>--------------------------------------------------]  65%
[==========================================================================================>-------------------------------------------------]  65%
[===========================================================================================>------------------------------------------------]  65%
[===========================================================================================>------------------------------------------------]  66%
[============================================================================================>-----------------------------------------------]  66%
[============================================================================================>-----------------------------------------------]  67%
[=============================================================================================>----------------------------------------------]  67%
[==============================================================================================>---------------------------------------------]  68%
[===============================================================================================>--------------------------------------------]  68%
[===============================================================================================>--------------------------------------------]  69%
[================================================================================================>-------------------------------------------]  69%
[================================================================================================>-------------------------------------------]  70%
[=================================================================================================>------------------------------------------]  70%
[==================================================================================================>-----------------------------------------]  70%
[==================================================================================================>-----------------------------------------]  71%
[===================================================================================================>----------------------------------------]  71%
[===================================================================================================>----------------------------------------]  72%
[====================================================================================================>---------------------------------------]  72%
[=====================================================================================================>--------------------------------------]  73%
[======================================================================================================>-------------------------------------]  73%
[======================================================================================================>-------------------------------------]  74%
[=======================================================================================================>------------------------------------]  74%
[=======================================================================================================>------------------------------------]  75%
[========================================================================================================>-----------------------------------]  75%
[=========================================================================================================>----------------------------------]  75%
[=========================================================================================================>----------------------------------]  76%
[==========================================================================================================>---------------------------------]  76%
[==========================================================================================================>---------------------------------]  77%
[===========================================================================================================>--------------------------------]  77%
[============================================================================================================>-------------------------------]  78%
[=============================================================================================================>------------------------------]  78%
[=============================================================================================================>------------------------------]  79%
[==============================================================================================================>-----------------------------]  79%
[==============================================================================================================>-----------------------------]  80%
[===============================================================================================================>----------------------------]  80%
[================================================================================================================>---------------------------]  80%
[================================================================================================================>---------------------------]  81%
[=================================================================================================================>--------------------------]  81%
[=================================================================================================================>--------------------------]  82%
[==================================================================================================================>-------------------------]  82%
[===================================================================================================================>------------------------]  83%
[====================================================================================================================>-----------------------]  83%
[====================================================================================================================>-----------------------]  84%
[=====================================================================================================================>----------------------]  84%
[=====================================================================================================================>----------------------]  85%
[======================================================================================================================>---------------------]  85%
[=======================================================================================================================>--------------------]  85%
[=======================================================================================================================>--------------------]  86%
[========================================================================================================================>-------------------]  86%
[========================================================================================================================>-------------------]  87%
[=========================================================================================================================>------------------]  87%
[==========================================================================================================================>-----------------]  88%
[===========================================================================================================================>----------------]  88%
[===========================================================================================================================>----------------]  89%
[============================================================================================================================>---------------]  89%
[============================================================================================================================>---------------]  90%
[=============================================================================================================================>--------------]  90%
[==============================================================================================================================>-------------]  90%
[==============================================================================================================================>-------------]  91%
[===============================================================================================================================>------------]  91%
[===============================================================================================================================>------------]  92%
[================================================================================================================================>-----------]  92%
[=================================================================================================================================>----------]  93%
[==================================================================================================================================>---------]  93%
[==================================================================================================================================>---------]  94%
[===================================================================================================================================>--------]  94%
[===================================================================================================================================>--------]  95%
[====================================================================================================================================>-------]  95%
[=====================================================================================================================================>------]  95%
[=====================================================================================================================================>------]  96%
[======================================================================================================================================>-----]  96%
[======================================================================================================================================>-----]  97%
[=======================================================================================================================================>----]  97%
[========================================================================================================================================>---]  98%
[=========================================================================================================================================>--]  98%
[=========================================================================================================================================>--]  99%
[==========================================================================================================================================>-]  99%
[==========================================================================================================================================>-] 100%
[===========================================================================================================================================>] 100%
[============================================================================================================================================] 100%
                                                                                                                                                   

Running diff_rhyth_AD_severity()
NOTE: Using batches in regression.
Warning: The `x` argument of `as_tibble.matrix()` must have unique column names if `.name_repair` is omitted as of tibble 2.0.0.
Using compatibility `.name_repair`.

Running diff_rhyth() on seedlist of size: 2287
NOTE: Using batches in regression.

[==>-----------------------------------------------------------------------------------------------------------------------------------------]   2%
[===>----------------------------------------------------------------------------------------------------------------------------------------]   3%
[====>---------------------------------------------------------------------------------------------------------------------------------------]   3%
[====>---------------------------------------------------------------------------------------------------------------------------------------]   4%
[=====>--------------------------------------------------------------------------------------------------------------------------------------]   4%
[=====>--------------------------------------------------------------------------------------------------------------------------------------]   5%
[======>-------------------------------------------------------------------------------------------------------------------------------------]   5%
[=======>------------------------------------------------------------------------------------------------------------------------------------]   5%
[=======>------------------------------------------------------------------------------------------------------------------------------------]   6%
[========>-----------------------------------------------------------------------------------------------------------------------------------]   6%
[========>-----------------------------------------------------------------------------------------------------------------------------------]   7%
[=========>----------------------------------------------------------------------------------------------------------------------------------]   7%
[==========>---------------------------------------------------------------------------------------------------------------------------------]   8%
[===========>--------------------------------------------------------------------------------------------------------------------------------]   8%
[===========>--------------------------------------------------------------------------------------------------------------------------------]   9%
[============>-------------------------------------------------------------------------------------------------------------------------------]   9%
[============>-------------------------------------------------------------------------------------------------------------------------------]  10%
[=============>------------------------------------------------------------------------------------------------------------------------------]  10%
[==============>-----------------------------------------------------------------------------------------------------------------------------]  10%
[==============>-----------------------------------------------------------------------------------------------------------------------------]  11%
[===============>----------------------------------------------------------------------------------------------------------------------------]  11%
[===============>----------------------------------------------------------------------------------------------------------------------------]  12%
[================>---------------------------------------------------------------------------------------------------------------------------]  12%
[=================>--------------------------------------------------------------------------------------------------------------------------]  13%
[==================>-------------------------------------------------------------------------------------------------------------------------]  13%
[==================>-------------------------------------------------------------------------------------------------------------------------]  14%
[===================>------------------------------------------------------------------------------------------------------------------------]  14%
[===================>------------------------------------------------------------------------------------------------------------------------]  15%
[====================>-----------------------------------------------------------------------------------------------------------------------]  15%
[=====================>----------------------------------------------------------------------------------------------------------------------]  15%
[=====================>----------------------------------------------------------------------------------------------------------------------]  16%
[======================>---------------------------------------------------------------------------------------------------------------------]  16%
[======================>---------------------------------------------------------------------------------------------------------------------]  17%
[=======================>--------------------------------------------------------------------------------------------------------------------]  17%
[========================>-------------------------------------------------------------------------------------------------------------------]  18%
[=========================>------------------------------------------------------------------------------------------------------------------]  18%
[=========================>------------------------------------------------------------------------------------------------------------------]  19%
[==========================>-----------------------------------------------------------------------------------------------------------------]  19%
[==========================>-----------------------------------------------------------------------------------------------------------------]  20%
[===========================>----------------------------------------------------------------------------------------------------------------]  20%
[============================>---------------------------------------------------------------------------------------------------------------]  20%
[============================>---------------------------------------------------------------------------------------------------------------]  21%
[=============================>--------------------------------------------------------------------------------------------------------------]  21%
[=============================>--------------------------------------------------------------------------------------------------------------]  22%
[==============================>-------------------------------------------------------------------------------------------------------------]  22%
[===============================>------------------------------------------------------------------------------------------------------------]  23%
[================================>-----------------------------------------------------------------------------------------------------------]  23%
[================================>-----------------------------------------------------------------------------------------------------------]  24%
[=================================>----------------------------------------------------------------------------------------------------------]  24%
[=================================>----------------------------------------------------------------------------------------------------------]  25%
[==================================>---------------------------------------------------------------------------------------------------------]  25%
[===================================>--------------------------------------------------------------------------------------------------------]  25%
[===================================>--------------------------------------------------------------------------------------------------------]  26%
[====================================>-------------------------------------------------------------------------------------------------------]  26%
[====================================>-------------------------------------------------------------------------------------------------------]  27%
[=====================================>------------------------------------------------------------------------------------------------------]  27%
[======================================>-----------------------------------------------------------------------------------------------------]  28%
[=======================================>----------------------------------------------------------------------------------------------------]  28%
[=======================================>----------------------------------------------------------------------------------------------------]  29%
[========================================>---------------------------------------------------------------------------------------------------]  29%
[========================================>---------------------------------------------------------------------------------------------------]  30%
[=========================================>--------------------------------------------------------------------------------------------------]  30%
[==========================================>-------------------------------------------------------------------------------------------------]  30%
[==========================================>-------------------------------------------------------------------------------------------------]  31%
[===========================================>------------------------------------------------------------------------------------------------]  31%
[===========================================>------------------------------------------------------------------------------------------------]  32%
[============================================>-----------------------------------------------------------------------------------------------]  32%
[=============================================>----------------------------------------------------------------------------------------------]  33%
[==============================================>---------------------------------------------------------------------------------------------]  33%
[==============================================>---------------------------------------------------------------------------------------------]  34%
[===============================================>--------------------------------------------------------------------------------------------]  34%
[===============================================>--------------------------------------------------------------------------------------------]  35%
[================================================>-------------------------------------------------------------------------------------------]  35%
[=================================================>------------------------------------------------------------------------------------------]  35%
[=================================================>------------------------------------------------------------------------------------------]  36%
[==================================================>-----------------------------------------------------------------------------------------]  36%
[==================================================>-----------------------------------------------------------------------------------------]  37%
[===================================================>----------------------------------------------------------------------------------------]  37%
[====================================================>---------------------------------------------------------------------------------------]  38%
[=====================================================>--------------------------------------------------------------------------------------]  38%
[=====================================================>--------------------------------------------------------------------------------------]  39%
[======================================================>-------------------------------------------------------------------------------------]  39%
[======================================================>-------------------------------------------------------------------------------------]  40%
[=======================================================>------------------------------------------------------------------------------------]  40%
[========================================================>-----------------------------------------------------------------------------------]  40%
[========================================================>-----------------------------------------------------------------------------------]  41%
[=========================================================>----------------------------------------------------------------------------------]  41%
[=========================================================>----------------------------------------------------------------------------------]  42%
[==========================================================>---------------------------------------------------------------------------------]  42%
[===========================================================>--------------------------------------------------------------------------------]  43%
[============================================================>-------------------------------------------------------------------------------]  43%
[============================================================>-------------------------------------------------------------------------------]  44%
[=============================================================>------------------------------------------------------------------------------]  44%
[=============================================================>------------------------------------------------------------------------------]  45%
[==============================================================>-----------------------------------------------------------------------------]  45%
[===============================================================>----------------------------------------------------------------------------]  45%
[===============================================================>----------------------------------------------------------------------------]  46%
[================================================================>---------------------------------------------------------------------------]  46%
[================================================================>---------------------------------------------------------------------------]  47%
[=================================================================>--------------------------------------------------------------------------]  47%
[==================================================================>-------------------------------------------------------------------------]  48%
[===================================================================>------------------------------------------------------------------------]  48%
[===================================================================>------------------------------------------------------------------------]  49%
[====================================================================>-----------------------------------------------------------------------]  49%
[====================================================================>-----------------------------------------------------------------------]  50%
[=====================================================================>----------------------------------------------------------------------]  50%
[======================================================================>---------------------------------------------------------------------]  50%
[======================================================================>---------------------------------------------------------------------]  51%
[=======================================================================>--------------------------------------------------------------------]  51%
[=======================================================================>--------------------------------------------------------------------]  52%
[========================================================================>-------------------------------------------------------------------]  52%
[=========================================================================>------------------------------------------------------------------]  53%
[==========================================================================>-----------------------------------------------------------------]  53%
[==========================================================================>-----------------------------------------------------------------]  54%
[===========================================================================>----------------------------------------------------------------]  54%
[===========================================================================>----------------------------------------------------------------]  55%
[============================================================================>---------------------------------------------------------------]  55%
[=============================================================================>--------------------------------------------------------------]  55%
[=============================================================================>--------------------------------------------------------------]  56%
[==============================================================================>-------------------------------------------------------------]  56%
[==============================================================================>-------------------------------------------------------------]  57%
[===============================================================================>------------------------------------------------------------]  57%
[================================================================================>-----------------------------------------------------------]  58%
[=================================================================================>----------------------------------------------------------]  58%
[=================================================================================>----------------------------------------------------------]  59%
[==================================================================================>---------------------------------------------------------]  59%
[==================================================================================>---------------------------------------------------------]  60%
[===================================================================================>--------------------------------------------------------]  60%
[====================================================================================>-------------------------------------------------------]  60%
[====================================================================================>-------------------------------------------------------]  61%
[=====================================================================================>------------------------------------------------------]  61%
[=====================================================================================>------------------------------------------------------]  62%
[======================================================================================>-----------------------------------------------------]  62%
[=======================================================================================>----------------------------------------------------]  63%
[========================================================================================>---------------------------------------------------]  63%
[========================================================================================>---------------------------------------------------]  64%
[=========================================================================================>--------------------------------------------------]  64%
[=========================================================================================>--------------------------------------------------]  65%
[==========================================================================================>-------------------------------------------------]  65%
[===========================================================================================>------------------------------------------------]  65%
[===========================================================================================>------------------------------------------------]  66%
[============================================================================================>-----------------------------------------------]  66%
[============================================================================================>-----------------------------------------------]  67%
[=============================================================================================>----------------------------------------------]  67%
[==============================================================================================>---------------------------------------------]  68%
[===============================================================================================>--------------------------------------------]  68%
[===============================================================================================>--------------------------------------------]  69%
[================================================================================================>-------------------------------------------]  69%
[================================================================================================>-------------------------------------------]  70%
[=================================================================================================>------------------------------------------]  70%
[==================================================================================================>-----------------------------------------]  70%
[==================================================================================================>-----------------------------------------]  71%
[===================================================================================================>----------------------------------------]  71%
[===================================================================================================>----------------------------------------]  72%
[====================================================================================================>---------------------------------------]  72%
[=====================================================================================================>--------------------------------------]  73%
[======================================================================================================>-------------------------------------]  73%
[======================================================================================================>-------------------------------------]  74%
[=======================================================================================================>------------------------------------]  74%
[=======================================================================================================>------------------------------------]  75%
[========================================================================================================>-----------------------------------]  75%
[=========================================================================================================>----------------------------------]  75%
[=========================================================================================================>----------------------------------]  76%
[==========================================================================================================>---------------------------------]  76%
[==========================================================================================================>---------------------------------]  77%
[===========================================================================================================>--------------------------------]  77%
[============================================================================================================>-------------------------------]  78%
[=============================================================================================================>------------------------------]  78%
[=============================================================================================================>------------------------------]  79%
[==============================================================================================================>-----------------------------]  79%
[==============================================================================================================>-----------------------------]  80%
[===============================================================================================================>----------------------------]  80%
[================================================================================================================>---------------------------]  80%
[================================================================================================================>---------------------------]  81%
[=================================================================================================================>--------------------------]  81%
[=================================================================================================================>--------------------------]  82%
[==================================================================================================================>-------------------------]  82%
[===================================================================================================================>------------------------]  83%
[====================================================================================================================>-----------------------]  83%
[====================================================================================================================>-----------------------]  84%
[=====================================================================================================================>----------------------]  84%
[=====================================================================================================================>----------------------]  85%
[======================================================================================================================>---------------------]  85%
[=======================================================================================================================>--------------------]  85%
[=======================================================================================================================>--------------------]  86%
[========================================================================================================================>-------------------]  86%
[========================================================================================================================>-------------------]  87%
[=========================================================================================================================>------------------]  87%
[==========================================================================================================================>-----------------]  88%
[===========================================================================================================================>----------------]  88%
[===========================================================================================================================>----------------]  89%
[============================================================================================================================>---------------]  89%
[============================================================================================================================>---------------]  90%
[=============================================================================================================================>--------------]  90%
[==============================================================================================================================>-------------]  90%
[==============================================================================================================================>-------------]  91%
[===============================================================================================================================>------------]  91%
[===============================================================================================================================>------------]  92%
[================================================================================================================================>-----------]  92%
[=================================================================================================================================>----------]  93%
[==================================================================================================================================>---------]  93%
[==================================================================================================================================>---------]  94%
[===================================================================================================================================>--------]  94%
[===================================================================================================================================>--------]  95%
[====================================================================================================================================>-------]  95%
[=====================================================================================================================================>------]  95%
[=====================================================================================================================================>------]  96%
[======================================================================================================================================>-----]  96%
[======================================================================================================================================>-----]  97%
[=======================================================================================================================================>----]  97%
[========================================================================================================================================>---]  98%
[=========================================================================================================================================>--]  98%
[=========================================================================================================================================>--]  99%
[==========================================================================================================================================>-]  99%
[==========================================================================================================================================>-] 100%
[===========================================================================================================================================>] 100%
[============================================================================================================================================] 100%
                                                                                                                                                   

Running diff_rhyth() on seedlist of size: 831
NOTE: Using batches in regression.

[=======>------------------------------------------------------------------------------------------------------------------------------------]   6%
[========>-----------------------------------------------------------------------------------------------------------------------------------]   6%
[========>-----------------------------------------------------------------------------------------------------------------------------------]   7%
[=========>----------------------------------------------------------------------------------------------------------------------------------]   7%
[==========>---------------------------------------------------------------------------------------------------------------------------------]   8%
[===========>--------------------------------------------------------------------------------------------------------------------------------]   8%
[===========>--------------------------------------------------------------------------------------------------------------------------------]   9%
[============>-------------------------------------------------------------------------------------------------------------------------------]   9%
[============>-------------------------------------------------------------------------------------------------------------------------------]  10%
[=============>------------------------------------------------------------------------------------------------------------------------------]  10%
[==============>-----------------------------------------------------------------------------------------------------------------------------]  10%
[==============>-----------------------------------------------------------------------------------------------------------------------------]  11%
[===============>----------------------------------------------------------------------------------------------------------------------------]  11%
[===============>----------------------------------------------------------------------------------------------------------------------------]  12%
[================>---------------------------------------------------------------------------------------------------------------------------]  12%
[=================>--------------------------------------------------------------------------------------------------------------------------]  13%
[==================>-------------------------------------------------------------------------------------------------------------------------]  13%
[==================>-------------------------------------------------------------------------------------------------------------------------]  14%
[===================>------------------------------------------------------------------------------------------------------------------------]  14%
[===================>------------------------------------------------------------------------------------------------------------------------]  15%
[====================>-----------------------------------------------------------------------------------------------------------------------]  15%
[=====================>----------------------------------------------------------------------------------------------------------------------]  15%
[=====================>----------------------------------------------------------------------------------------------------------------------]  16%
[======================>---------------------------------------------------------------------------------------------------------------------]  16%
[======================>---------------------------------------------------------------------------------------------------------------------]  17%
[=======================>--------------------------------------------------------------------------------------------------------------------]  17%
[========================>-------------------------------------------------------------------------------------------------------------------]  18%
[=========================>------------------------------------------------------------------------------------------------------------------]  18%
[=========================>------------------------------------------------------------------------------------------------------------------]  19%
[==========================>-----------------------------------------------------------------------------------------------------------------]  19%
[==========================>-----------------------------------------------------------------------------------------------------------------]  20%
[===========================>----------------------------------------------------------------------------------------------------------------]  20%
[============================>---------------------------------------------------------------------------------------------------------------]  20%
[============================>---------------------------------------------------------------------------------------------------------------]  21%
[=============================>--------------------------------------------------------------------------------------------------------------]  21%
[=============================>--------------------------------------------------------------------------------------------------------------]  22%
[==============================>-------------------------------------------------------------------------------------------------------------]  22%
[===============================>------------------------------------------------------------------------------------------------------------]  23%
[================================>-----------------------------------------------------------------------------------------------------------]  23%
[================================>-----------------------------------------------------------------------------------------------------------]  24%
[=================================>----------------------------------------------------------------------------------------------------------]  24%
[=================================>----------------------------------------------------------------------------------------------------------]  25%
[==================================>---------------------------------------------------------------------------------------------------------]  25%
[===================================>--------------------------------------------------------------------------------------------------------]  25%
[===================================>--------------------------------------------------------------------------------------------------------]  26%
[====================================>-------------------------------------------------------------------------------------------------------]  26%
[====================================>-------------------------------------------------------------------------------------------------------]  27%
[=====================================>------------------------------------------------------------------------------------------------------]  27%
[======================================>-----------------------------------------------------------------------------------------------------]  28%
[=======================================>----------------------------------------------------------------------------------------------------]  28%
[=======================================>----------------------------------------------------------------------------------------------------]  29%
[========================================>---------------------------------------------------------------------------------------------------]  29%
[========================================>---------------------------------------------------------------------------------------------------]  30%
[=========================================>--------------------------------------------------------------------------------------------------]  30%
[==========================================>-------------------------------------------------------------------------------------------------]  30%
[==========================================>-------------------------------------------------------------------------------------------------]  31%
[===========================================>------------------------------------------------------------------------------------------------]  31%
[===========================================>------------------------------------------------------------------------------------------------]  32%
[============================================>-----------------------------------------------------------------------------------------------]  32%
[=============================================>----------------------------------------------------------------------------------------------]  33%
[==============================================>---------------------------------------------------------------------------------------------]  33%
[==============================================>---------------------------------------------------------------------------------------------]  34%
[===============================================>--------------------------------------------------------------------------------------------]  34%
[===============================================>--------------------------------------------------------------------------------------------]  35%
[================================================>-------------------------------------------------------------------------------------------]  35%
[=================================================>------------------------------------------------------------------------------------------]  35%
[=================================================>------------------------------------------------------------------------------------------]  36%
[==================================================>-----------------------------------------------------------------------------------------]  36%
[==================================================>-----------------------------------------------------------------------------------------]  37%
[===================================================>----------------------------------------------------------------------------------------]  37%
[====================================================>---------------------------------------------------------------------------------------]  38%
[=====================================================>--------------------------------------------------------------------------------------]  38%
[=====================================================>--------------------------------------------------------------------------------------]  39%
[======================================================>-------------------------------------------------------------------------------------]  39%
[======================================================>-------------------------------------------------------------------------------------]  40%
[=======================================================>------------------------------------------------------------------------------------]  40%
[========================================================>-----------------------------------------------------------------------------------]  40%
[========================================================>-----------------------------------------------------------------------------------]  41%
[=========================================================>----------------------------------------------------------------------------------]  41%
[=========================================================>----------------------------------------------------------------------------------]  42%
[==========================================================>---------------------------------------------------------------------------------]  42%
[===========================================================>--------------------------------------------------------------------------------]  43%
[============================================================>-------------------------------------------------------------------------------]  43%
[============================================================>-------------------------------------------------------------------------------]  44%
[=============================================================>------------------------------------------------------------------------------]  44%
[=============================================================>------------------------------------------------------------------------------]  45%
[==============================================================>-----------------------------------------------------------------------------]  45%
[===============================================================>----------------------------------------------------------------------------]  45%
[===============================================================>----------------------------------------------------------------------------]  46%
[================================================================>---------------------------------------------------------------------------]  46%
[================================================================>---------------------------------------------------------------------------]  47%
[=================================================================>--------------------------------------------------------------------------]  47%
[==================================================================>-------------------------------------------------------------------------]  48%
[===================================================================>------------------------------------------------------------------------]  48%
[===================================================================>------------------------------------------------------------------------]  49%
[====================================================================>-----------------------------------------------------------------------]  49%
[====================================================================>-----------------------------------------------------------------------]  50%
[=====================================================================>----------------------------------------------------------------------]  50%
[======================================================================>---------------------------------------------------------------------]  50%
[======================================================================>---------------------------------------------------------------------]  51%
[=======================================================================>--------------------------------------------------------------------]  51%
[=======================================================================>--------------------------------------------------------------------]  52%
[========================================================================>-------------------------------------------------------------------]  52%
[=========================================================================>------------------------------------------------------------------]  53%
[==========================================================================>-----------------------------------------------------------------]  53%
[==========================================================================>-----------------------------------------------------------------]  54%
[===========================================================================>----------------------------------------------------------------]  54%
[===========================================================================>----------------------------------------------------------------]  55%
[============================================================================>---------------------------------------------------------------]  55%
[=============================================================================>--------------------------------------------------------------]  55%
[=============================================================================>--------------------------------------------------------------]  56%
[==============================================================================>-------------------------------------------------------------]  56%
[==============================================================================>-------------------------------------------------------------]  57%
[===============================================================================>------------------------------------------------------------]  57%
[================================================================================>-----------------------------------------------------------]  58%
[=================================================================================>----------------------------------------------------------]  58%
[=================================================================================>----------------------------------------------------------]  59%
[==================================================================================>---------------------------------------------------------]  59%
[==================================================================================>---------------------------------------------------------]  60%
[===================================================================================>--------------------------------------------------------]  60%
[====================================================================================>-------------------------------------------------------]  60%
[====================================================================================>-------------------------------------------------------]  61%
[=====================================================================================>------------------------------------------------------]  61%
[=====================================================================================>------------------------------------------------------]  62%
[======================================================================================>-----------------------------------------------------]  62%
[=======================================================================================>----------------------------------------------------]  63%
[========================================================================================>---------------------------------------------------]  63%
[========================================================================================>---------------------------------------------------]  64%
[=========================================================================================>--------------------------------------------------]  64%
[=========================================================================================>--------------------------------------------------]  65%
[==========================================================================================>-------------------------------------------------]  65%
[===========================================================================================>------------------------------------------------]  65%
[===========================================================================================>------------------------------------------------]  66%
[============================================================================================>-----------------------------------------------]  66%
[============================================================================================>-----------------------------------------------]  67%
[=============================================================================================>----------------------------------------------]  67%
[==============================================================================================>---------------------------------------------]  68%
[===============================================================================================>--------------------------------------------]  68%
[===============================================================================================>--------------------------------------------]  69%
[================================================================================================>-------------------------------------------]  69%
[================================================================================================>-------------------------------------------]  70%
[=================================================================================================>------------------------------------------]  70%
[==================================================================================================>-----------------------------------------]  70%
[==================================================================================================>-----------------------------------------]  71%
[===================================================================================================>----------------------------------------]  71%
[===================================================================================================>----------------------------------------]  72%
[====================================================================================================>---------------------------------------]  72%
[=====================================================================================================>--------------------------------------]  73%
[======================================================================================================>-------------------------------------]  73%
[======================================================================================================>-------------------------------------]  74%
[=======================================================================================================>------------------------------------]  74%
[=======================================================================================================>------------------------------------]  75%
[========================================================================================================>-----------------------------------]  75%
[=========================================================================================================>----------------------------------]  75%
[=========================================================================================================>----------------------------------]  76%
[==========================================================================================================>---------------------------------]  76%
[==========================================================================================================>---------------------------------]  77%
[===========================================================================================================>--------------------------------]  77%
[============================================================================================================>-------------------------------]  78%
[=============================================================================================================>------------------------------]  78%
[=============================================================================================================>------------------------------]  79%
[==============================================================================================================>-----------------------------]  79%
[==============================================================================================================>-----------------------------]  80%
[===============================================================================================================>----------------------------]  80%
[================================================================================================================>---------------------------]  80%
[================================================================================================================>---------------------------]  81%
[=================================================================================================================>--------------------------]  81%
[=================================================================================================================>--------------------------]  82%
[==================================================================================================================>-------------------------]  82%
[===================================================================================================================>------------------------]  83%
[====================================================================================================================>-----------------------]  83%
[====================================================================================================================>-----------------------]  84%
[=====================================================================================================================>----------------------]  84%
[=====================================================================================================================>----------------------]  85%
[======================================================================================================================>---------------------]  85%
[=======================================================================================================================>--------------------]  85%
[=======================================================================================================================>--------------------]  86%
[========================================================================================================================>-------------------]  86%
[========================================================================================================================>-------------------]  87%
[=========================================================================================================================>------------------]  87%
[==========================================================================================================================>-----------------]  88%
[===========================================================================================================================>----------------]  88%
[===========================================================================================================================>----------------]  89%
[============================================================================================================================>---------------]  89%
[============================================================================================================================>---------------]  90%
[=============================================================================================================================>--------------]  90%
[==============================================================================================================================>-------------]  90%
[==============================================================================================================================>-------------]  91%
[===============================================================================================================================>------------]  91%
[===============================================================================================================================>------------]  92%
[================================================================================================================================>-----------]  92%
[=================================================================================================================================>----------]  93%
[==================================================================================================================================>---------]  93%
[==================================================================================================================================>---------]  94%
[===================================================================================================================================>--------]  94%
[===================================================================================================================================>--------]  95%
[====================================================================================================================================>-------]  95%
[=====================================================================================================================================>------]  95%
[=====================================================================================================================================>------]  96%
[======================================================================================================================================>-----]  96%
[======================================================================================================================================>-----]  97%
[=======================================================================================================================================>----]  97%
[========================================================================================================================================>---]  98%
[=========================================================================================================================================>--]  98%
[=========================================================================================================================================>--]  99%
[==========================================================================================================================================>-]  99%
[==========================================================================================================================================>-] 100%
[===========================================================================================================================================>] 100%
[============================================================================================================================================] 100%
                                                                                                                                                   

Running diff_rhyth() on seedlist of size: 10771
NOTE: Using batches in regression.

[>-------------------------------------------------------------------------------------------------------------------------------------------]   0%
[>-------------------------------------------------------------------------------------------------------------------------------------------]   1%
[=>------------------------------------------------------------------------------------------------------------------------------------------]   1%
[=>------------------------------------------------------------------------------------------------------------------------------------------]   2%
[==>-----------------------------------------------------------------------------------------------------------------------------------------]   2%
[===>----------------------------------------------------------------------------------------------------------------------------------------]   3%
[====>---------------------------------------------------------------------------------------------------------------------------------------]   3%
[====>---------------------------------------------------------------------------------------------------------------------------------------]   4%
[=====>--------------------------------------------------------------------------------------------------------------------------------------]   4%
[=====>--------------------------------------------------------------------------------------------------------------------------------------]   5%
[======>-------------------------------------------------------------------------------------------------------------------------------------]   5%
[=======>------------------------------------------------------------------------------------------------------------------------------------]   5%
[=======>------------------------------------------------------------------------------------------------------------------------------------]   6%
[========>-----------------------------------------------------------------------------------------------------------------------------------]   6%
[========>-----------------------------------------------------------------------------------------------------------------------------------]   7%
[=========>----------------------------------------------------------------------------------------------------------------------------------]   7%
[==========>---------------------------------------------------------------------------------------------------------------------------------]   8%
[===========>--------------------------------------------------------------------------------------------------------------------------------]   8%
[===========>--------------------------------------------------------------------------------------------------------------------------------]   9%
[============>-------------------------------------------------------------------------------------------------------------------------------]   9%
[============>-------------------------------------------------------------------------------------------------------------------------------]  10%
[=============>------------------------------------------------------------------------------------------------------------------------------]  10%
[==============>-----------------------------------------------------------------------------------------------------------------------------]  10%
[==============>-----------------------------------------------------------------------------------------------------------------------------]  11%
[===============>----------------------------------------------------------------------------------------------------------------------------]  11%
[===============>----------------------------------------------------------------------------------------------------------------------------]  12%
[================>---------------------------------------------------------------------------------------------------------------------------]  12%
[=================>--------------------------------------------------------------------------------------------------------------------------]  13%
[==================>-------------------------------------------------------------------------------------------------------------------------]  13%
[==================>-------------------------------------------------------------------------------------------------------------------------]  14%
[===================>------------------------------------------------------------------------------------------------------------------------]  14%
[===================>------------------------------------------------------------------------------------------------------------------------]  15%
[====================>-----------------------------------------------------------------------------------------------------------------------]  15%
[=====================>----------------------------------------------------------------------------------------------------------------------]  15%
[=====================>----------------------------------------------------------------------------------------------------------------------]  16%
[======================>---------------------------------------------------------------------------------------------------------------------]  16%
[======================>---------------------------------------------------------------------------------------------------------------------]  17%
[=======================>--------------------------------------------------------------------------------------------------------------------]  17%
[========================>-------------------------------------------------------------------------------------------------------------------]  18%
[=========================>------------------------------------------------------------------------------------------------------------------]  18%
[=========================>------------------------------------------------------------------------------------------------------------------]  19%
[==========================>-----------------------------------------------------------------------------------------------------------------]  19%
[==========================>-----------------------------------------------------------------------------------------------------------------]  20%
[===========================>----------------------------------------------------------------------------------------------------------------]  20%
[============================>---------------------------------------------------------------------------------------------------------------]  20%
[============================>---------------------------------------------------------------------------------------------------------------]  21%
[=============================>--------------------------------------------------------------------------------------------------------------]  21%
[=============================>--------------------------------------------------------------------------------------------------------------]  22%
[==============================>-------------------------------------------------------------------------------------------------------------]  22%
[===============================>------------------------------------------------------------------------------------------------------------]  23%
[================================>-----------------------------------------------------------------------------------------------------------]  23%
[================================>-----------------------------------------------------------------------------------------------------------]  24%
[=================================>----------------------------------------------------------------------------------------------------------]  24%
[=================================>----------------------------------------------------------------------------------------------------------]  25%
[==================================>---------------------------------------------------------------------------------------------------------]  25%
[===================================>--------------------------------------------------------------------------------------------------------]  25%
[===================================>--------------------------------------------------------------------------------------------------------]  26%
[====================================>-------------------------------------------------------------------------------------------------------]  26%
[====================================>-------------------------------------------------------------------------------------------------------]  27%
[=====================================>------------------------------------------------------------------------------------------------------]  27%
[======================================>-----------------------------------------------------------------------------------------------------]  28%
[=======================================>----------------------------------------------------------------------------------------------------]  28%
[=======================================>----------------------------------------------------------------------------------------------------]  29%
[========================================>---------------------------------------------------------------------------------------------------]  29%
[========================================>---------------------------------------------------------------------------------------------------]  30%
[=========================================>--------------------------------------------------------------------------------------------------]  30%
[==========================================>-------------------------------------------------------------------------------------------------]  30%
[==========================================>-------------------------------------------------------------------------------------------------]  31%
[===========================================>------------------------------------------------------------------------------------------------]  31%
[===========================================>------------------------------------------------------------------------------------------------]  32%
[============================================>-----------------------------------------------------------------------------------------------]  32%
[=============================================>----------------------------------------------------------------------------------------------]  33%
[==============================================>---------------------------------------------------------------------------------------------]  33%
[==============================================>---------------------------------------------------------------------------------------------]  34%
[===============================================>--------------------------------------------------------------------------------------------]  34%
[===============================================>--------------------------------------------------------------------------------------------]  35%
[================================================>-------------------------------------------------------------------------------------------]  35%
[=================================================>------------------------------------------------------------------------------------------]  35%
[=================================================>------------------------------------------------------------------------------------------]  36%
[==================================================>-----------------------------------------------------------------------------------------]  36%
[==================================================>-----------------------------------------------------------------------------------------]  37%
[===================================================>----------------------------------------------------------------------------------------]  37%
[====================================================>---------------------------------------------------------------------------------------]  38%
[=====================================================>--------------------------------------------------------------------------------------]  38%
[=====================================================>--------------------------------------------------------------------------------------]  39%
[======================================================>-------------------------------------------------------------------------------------]  39%
[======================================================>-------------------------------------------------------------------------------------]  40%
[=======================================================>------------------------------------------------------------------------------------]  40%
[========================================================>-----------------------------------------------------------------------------------]  40%
[========================================================>-----------------------------------------------------------------------------------]  41%
[=========================================================>----------------------------------------------------------------------------------]  41%
[=========================================================>----------------------------------------------------------------------------------]  42%
[==========================================================>---------------------------------------------------------------------------------]  42%
[===========================================================>--------------------------------------------------------------------------------]  43%
[============================================================>-------------------------------------------------------------------------------]  43%
[============================================================>-------------------------------------------------------------------------------]  44%
[=============================================================>------------------------------------------------------------------------------]  44%
[=============================================================>------------------------------------------------------------------------------]  45%
[==============================================================>-----------------------------------------------------------------------------]  45%
[===============================================================>----------------------------------------------------------------------------]  45%
[===============================================================>----------------------------------------------------------------------------]  46%
[================================================================>---------------------------------------------------------------------------]  46%
[================================================================>---------------------------------------------------------------------------]  47%
[=================================================================>--------------------------------------------------------------------------]  47%
[==================================================================>-------------------------------------------------------------------------]  48%
[===================================================================>------------------------------------------------------------------------]  48%
[===================================================================>------------------------------------------------------------------------]  49%
[====================================================================>-----------------------------------------------------------------------]  49%
[====================================================================>-----------------------------------------------------------------------]  50%
[=====================================================================>----------------------------------------------------------------------]  50%
[======================================================================>---------------------------------------------------------------------]  50%
[======================================================================>---------------------------------------------------------------------]  51%
[=======================================================================>--------------------------------------------------------------------]  51%
[=======================================================================>--------------------------------------------------------------------]  52%
[========================================================================>-------------------------------------------------------------------]  52%
[=========================================================================>------------------------------------------------------------------]  53%
[==========================================================================>-----------------------------------------------------------------]  53%
[==========================================================================>-----------------------------------------------------------------]  54%
[===========================================================================>----------------------------------------------------------------]  54%
[===========================================================================>----------------------------------------------------------------]  55%
[============================================================================>---------------------------------------------------------------]  55%
[=============================================================================>--------------------------------------------------------------]  55%
[=============================================================================>--------------------------------------------------------------]  56%
[==============================================================================>-------------------------------------------------------------]  56%
[==============================================================================>-------------------------------------------------------------]  57%
[===============================================================================>------------------------------------------------------------]  57%
[================================================================================>-----------------------------------------------------------]  58%
[=================================================================================>----------------------------------------------------------]  58%
[=================================================================================>----------------------------------------------------------]  59%
[==================================================================================>---------------------------------------------------------]  59%
[==================================================================================>---------------------------------------------------------]  60%
[===================================================================================>--------------------------------------------------------]  60%
[====================================================================================>-------------------------------------------------------]  60%
[====================================================================================>-------------------------------------------------------]  61%
[=====================================================================================>------------------------------------------------------]  61%
[=====================================================================================>------------------------------------------------------]  62%
[======================================================================================>-----------------------------------------------------]  62%
[=======================================================================================>----------------------------------------------------]  63%
[========================================================================================>---------------------------------------------------]  63%
[========================================================================================>---------------------------------------------------]  64%
[=========================================================================================>--------------------------------------------------]  64%
[=========================================================================================>--------------------------------------------------]  65%
[==========================================================================================>-------------------------------------------------]  65%
[===========================================================================================>------------------------------------------------]  65%
[===========================================================================================>------------------------------------------------]  66%
[============================================================================================>-----------------------------------------------]  66%
[============================================================================================>-----------------------------------------------]  67%
[=============================================================================================>----------------------------------------------]  67%
[==============================================================================================>---------------------------------------------]  68%
[===============================================================================================>--------------------------------------------]  68%
[===============================================================================================>--------------------------------------------]  69%
[================================================================================================>-------------------------------------------]  69%
[================================================================================================>-------------------------------------------]  70%
[=================================================================================================>------------------------------------------]  70%
[==================================================================================================>-----------------------------------------]  70%
[==================================================================================================>-----------------------------------------]  71%
[===================================================================================================>----------------------------------------]  71%
[===================================================================================================>----------------------------------------]  72%
[====================================================================================================>---------------------------------------]  72%
[=====================================================================================================>--------------------------------------]  73%
[======================================================================================================>-------------------------------------]  73%
[======================================================================================================>-------------------------------------]  74%
[=======================================================================================================>------------------------------------]  74%
[=======================================================================================================>------------------------------------]  75%
[========================================================================================================>-----------------------------------]  75%
[=========================================================================================================>----------------------------------]  75%
[=========================================================================================================>----------------------------------]  76%
[==========================================================================================================>---------------------------------]  76%
[==========================================================================================================>---------------------------------]  77%
[===========================================================================================================>--------------------------------]  77%
[============================================================================================================>-------------------------------]  78%
[=============================================================================================================>------------------------------]  78%
[=============================================================================================================>------------------------------]  79%
[==============================================================================================================>-----------------------------]  79%
[==============================================================================================================>-----------------------------]  80%
[===============================================================================================================>----------------------------]  80%
[================================================================================================================>---------------------------]  80%
[================================================================================================================>---------------------------]  81%
[=================================================================================================================>--------------------------]  81%
[=================================================================================================================>--------------------------]  82%
[==================================================================================================================>-------------------------]  82%
[===================================================================================================================>------------------------]  83%
[====================================================================================================================>-----------------------]  83%
[====================================================================================================================>-----------------------]  84%
[=====================================================================================================================>----------------------]  84%
[=====================================================================================================================>----------------------]  85%
[======================================================================================================================>---------------------]  85%
[=======================================================================================================================>--------------------]  85%
[=======================================================================================================================>--------------------]  86%
[========================================================================================================================>-------------------]  86%
[========================================================================================================================>-------------------]  87%
[=========================================================================================================================>------------------]  87%
[==========================================================================================================================>-----------------]  88%
[===========================================================================================================================>----------------]  88%
[===========================================================================================================================>----------------]  89%
[============================================================================================================================>---------------]  89%
[============================================================================================================================>---------------]  90%
[=============================================================================================================================>--------------]  90%
[==============================================================================================================================>-------------]  90%
[==============================================================================================================================>-------------]  91%
[===============================================================================================================================>------------]  91%
[===============================================================================================================================>------------]  92%
[================================================================================================================================>-----------]  92%
[=================================================================================================================================>----------]  93%
[==================================================================================================================================>---------]  93%
[==================================================================================================================================>---------]  94%
[===================================================================================================================================>--------]  94%
[===================================================================================================================================>--------]  95%
[====================================================================================================================================>-------]  95%
[=====================================================================================================================================>------]  95%
[=====================================================================================================================================>------]  96%
[======================================================================================================================================>-----]  96%
[======================================================================================================================================>-----]  97%
[=======================================================================================================================================>----]  97%
[========================================================================================================================================>---]  98%
[=========================================================================================================================================>--]  98%
[=========================================================================================================================================>--]  99%
[==========================================================================================================================================>-]  99%
[==========================================================================================================================================>-] 100%
[===========================================================================================================================================>] 100%
[============================================================================================================================================] 100%
                                                                                                                                                   

Running diff_rhyth() on seedlist of size: 2300
NOTE: Using batches in regression.

[==>-----------------------------------------------------------------------------------------------------------------------------------------]   2%
[===>----------------------------------------------------------------------------------------------------------------------------------------]   3%
[====>---------------------------------------------------------------------------------------------------------------------------------------]   3%
[====>---------------------------------------------------------------------------------------------------------------------------------------]   4%
[=====>--------------------------------------------------------------------------------------------------------------------------------------]   4%
[=====>--------------------------------------------------------------------------------------------------------------------------------------]   5%
[======>-------------------------------------------------------------------------------------------------------------------------------------]   5%
[=======>------------------------------------------------------------------------------------------------------------------------------------]   5%
[=======>------------------------------------------------------------------------------------------------------------------------------------]   6%
[========>-----------------------------------------------------------------------------------------------------------------------------------]   6%
[========>-----------------------------------------------------------------------------------------------------------------------------------]   7%
[=========>----------------------------------------------------------------------------------------------------------------------------------]   7%
[==========>---------------------------------------------------------------------------------------------------------------------------------]   8%
[===========>--------------------------------------------------------------------------------------------------------------------------------]   8%
[===========>--------------------------------------------------------------------------------------------------------------------------------]   9%
[============>-------------------------------------------------------------------------------------------------------------------------------]   9%
[============>-------------------------------------------------------------------------------------------------------------------------------]  10%
[=============>------------------------------------------------------------------------------------------------------------------------------]  10%
[==============>-----------------------------------------------------------------------------------------------------------------------------]  10%
[==============>-----------------------------------------------------------------------------------------------------------------------------]  11%
[===============>----------------------------------------------------------------------------------------------------------------------------]  11%
[===============>----------------------------------------------------------------------------------------------------------------------------]  12%
[================>---------------------------------------------------------------------------------------------------------------------------]  12%
[=================>--------------------------------------------------------------------------------------------------------------------------]  13%
[==================>-------------------------------------------------------------------------------------------------------------------------]  13%
[==================>-------------------------------------------------------------------------------------------------------------------------]  14%
[===================>------------------------------------------------------------------------------------------------------------------------]  14%
[===================>------------------------------------------------------------------------------------------------------------------------]  15%
[====================>-----------------------------------------------------------------------------------------------------------------------]  15%
[=====================>----------------------------------------------------------------------------------------------------------------------]  15%
[=====================>----------------------------------------------------------------------------------------------------------------------]  16%
[======================>---------------------------------------------------------------------------------------------------------------------]  16%
[======================>---------------------------------------------------------------------------------------------------------------------]  17%
[=======================>--------------------------------------------------------------------------------------------------------------------]  17%
[========================>-------------------------------------------------------------------------------------------------------------------]  18%
[=========================>------------------------------------------------------------------------------------------------------------------]  18%
[=========================>------------------------------------------------------------------------------------------------------------------]  19%
[==========================>-----------------------------------------------------------------------------------------------------------------]  19%
[==========================>-----------------------------------------------------------------------------------------------------------------]  20%
[===========================>----------------------------------------------------------------------------------------------------------------]  20%
[============================>---------------------------------------------------------------------------------------------------------------]  20%
[============================>---------------------------------------------------------------------------------------------------------------]  21%
[=============================>--------------------------------------------------------------------------------------------------------------]  21%
[=============================>--------------------------------------------------------------------------------------------------------------]  22%
[==============================>-------------------------------------------------------------------------------------------------------------]  22%
[===============================>------------------------------------------------------------------------------------------------------------]  23%
[================================>-----------------------------------------------------------------------------------------------------------]  23%
[================================>-----------------------------------------------------------------------------------------------------------]  24%
[=================================>----------------------------------------------------------------------------------------------------------]  24%
[=================================>----------------------------------------------------------------------------------------------------------]  25%
[==================================>---------------------------------------------------------------------------------------------------------]  25%
[===================================>--------------------------------------------------------------------------------------------------------]  25%
[===================================>--------------------------------------------------------------------------------------------------------]  26%
[====================================>-------------------------------------------------------------------------------------------------------]  26%
[====================================>-------------------------------------------------------------------------------------------------------]  27%
[=====================================>------------------------------------------------------------------------------------------------------]  27%
[======================================>-----------------------------------------------------------------------------------------------------]  28%
[=======================================>----------------------------------------------------------------------------------------------------]  28%
[=======================================>----------------------------------------------------------------------------------------------------]  29%
[========================================>---------------------------------------------------------------------------------------------------]  29%
[========================================>---------------------------------------------------------------------------------------------------]  30%
[=========================================>--------------------------------------------------------------------------------------------------]  30%
[==========================================>-------------------------------------------------------------------------------------------------]  30%
[==========================================>-------------------------------------------------------------------------------------------------]  31%
[===========================================>------------------------------------------------------------------------------------------------]  31%
[===========================================>------------------------------------------------------------------------------------------------]  32%
[============================================>-----------------------------------------------------------------------------------------------]  32%
[=============================================>----------------------------------------------------------------------------------------------]  33%
[==============================================>---------------------------------------------------------------------------------------------]  33%
[==============================================>---------------------------------------------------------------------------------------------]  34%
[===============================================>--------------------------------------------------------------------------------------------]  34%
[===============================================>--------------------------------------------------------------------------------------------]  35%
[================================================>-------------------------------------------------------------------------------------------]  35%
[=================================================>------------------------------------------------------------------------------------------]  35%
[=================================================>------------------------------------------------------------------------------------------]  36%
[==================================================>-----------------------------------------------------------------------------------------]  36%
[==================================================>-----------------------------------------------------------------------------------------]  37%
[===================================================>----------------------------------------------------------------------------------------]  37%
[====================================================>---------------------------------------------------------------------------------------]  38%
[=====================================================>--------------------------------------------------------------------------------------]  38%
[=====================================================>--------------------------------------------------------------------------------------]  39%
[======================================================>-------------------------------------------------------------------------------------]  39%
[======================================================>-------------------------------------------------------------------------------------]  40%
[=======================================================>------------------------------------------------------------------------------------]  40%
[========================================================>-----------------------------------------------------------------------------------]  40%
[========================================================>-----------------------------------------------------------------------------------]  41%
[=========================================================>----------------------------------------------------------------------------------]  41%
[=========================================================>----------------------------------------------------------------------------------]  42%
[==========================================================>---------------------------------------------------------------------------------]  42%
[===========================================================>--------------------------------------------------------------------------------]  43%
[============================================================>-------------------------------------------------------------------------------]  43%
[============================================================>-------------------------------------------------------------------------------]  44%
[=============================================================>------------------------------------------------------------------------------]  44%
[=============================================================>------------------------------------------------------------------------------]  45%
[==============================================================>-----------------------------------------------------------------------------]  45%
[===============================================================>----------------------------------------------------------------------------]  45%
[===============================================================>----------------------------------------------------------------------------]  46%
[================================================================>---------------------------------------------------------------------------]  46%
[================================================================>---------------------------------------------------------------------------]  47%
[=================================================================>--------------------------------------------------------------------------]  47%
[==================================================================>-------------------------------------------------------------------------]  48%
[===================================================================>------------------------------------------------------------------------]  48%
[===================================================================>------------------------------------------------------------------------]  49%
[====================================================================>-----------------------------------------------------------------------]  49%
[====================================================================>-----------------------------------------------------------------------]  50%
[=====================================================================>----------------------------------------------------------------------]  50%
[======================================================================>---------------------------------------------------------------------]  50%
[======================================================================>---------------------------------------------------------------------]  51%
[=======================================================================>--------------------------------------------------------------------]  51%
[=======================================================================>--------------------------------------------------------------------]  52%
[========================================================================>-------------------------------------------------------------------]  52%
[=========================================================================>------------------------------------------------------------------]  53%
[==========================================================================>-----------------------------------------------------------------]  53%
[==========================================================================>-----------------------------------------------------------------]  54%
[===========================================================================>----------------------------------------------------------------]  54%
[===========================================================================>----------------------------------------------------------------]  55%
[============================================================================>---------------------------------------------------------------]  55%
[=============================================================================>--------------------------------------------------------------]  55%
[=============================================================================>--------------------------------------------------------------]  56%
[==============================================================================>-------------------------------------------------------------]  56%
[==============================================================================>-------------------------------------------------------------]  57%
[===============================================================================>------------------------------------------------------------]  57%
[================================================================================>-----------------------------------------------------------]  58%
[=================================================================================>----------------------------------------------------------]  58%
[=================================================================================>----------------------------------------------------------]  59%
[==================================================================================>---------------------------------------------------------]  59%
[==================================================================================>---------------------------------------------------------]  60%
[===================================================================================>--------------------------------------------------------]  60%
[====================================================================================>-------------------------------------------------------]  60%
[====================================================================================>-------------------------------------------------------]  61%
[=====================================================================================>------------------------------------------------------]  61%
[=====================================================================================>------------------------------------------------------]  62%
[======================================================================================>-----------------------------------------------------]  62%
[=======================================================================================>----------------------------------------------------]  63%
[========================================================================================>---------------------------------------------------]  63%
[========================================================================================>---------------------------------------------------]  64%
[=========================================================================================>--------------------------------------------------]  64%
[=========================================================================================>--------------------------------------------------]  65%
[==========================================================================================>-------------------------------------------------]  65%
[===========================================================================================>------------------------------------------------]  65%
[===========================================================================================>------------------------------------------------]  66%
[============================================================================================>-----------------------------------------------]  66%
[============================================================================================>-----------------------------------------------]  67%
[=============================================================================================>----------------------------------------------]  67%
[==============================================================================================>---------------------------------------------]  68%
[===============================================================================================>--------------------------------------------]  68%
[===============================================================================================>--------------------------------------------]  69%
[================================================================================================>-------------------------------------------]  69%
[================================================================================================>-------------------------------------------]  70%
[=================================================================================================>------------------------------------------]  70%
[==================================================================================================>-----------------------------------------]  70%
[==================================================================================================>-----------------------------------------]  71%
[===================================================================================================>----------------------------------------]  71%
[===================================================================================================>----------------------------------------]  72%
[====================================================================================================>---------------------------------------]  72%
[=====================================================================================================>--------------------------------------]  73%
[======================================================================================================>-------------------------------------]  73%
[======================================================================================================>-------------------------------------]  74%
[=======================================================================================================>------------------------------------]  74%
[=======================================================================================================>------------------------------------]  75%
[========================================================================================================>-----------------------------------]  75%
[=========================================================================================================>----------------------------------]  75%
[=========================================================================================================>----------------------------------]  76%
[==========================================================================================================>---------------------------------]  76%
[==========================================================================================================>---------------------------------]  77%
[===========================================================================================================>--------------------------------]  77%
[============================================================================================================>-------------------------------]  78%
[=============================================================================================================>------------------------------]  78%
[=============================================================================================================>------------------------------]  79%
[==============================================================================================================>-----------------------------]  79%
[==============================================================================================================>-----------------------------]  80%
[===============================================================================================================>----------------------------]  80%
[================================================================================================================>---------------------------]  80%
[================================================================================================================>---------------------------]  81%
[=================================================================================================================>--------------------------]  81%
[=================================================================================================================>--------------------------]  82%
[==================================================================================================================>-------------------------]  82%
[===================================================================================================================>------------------------]  83%
[====================================================================================================================>-----------------------]  83%
[====================================================================================================================>-----------------------]  84%
[=====================================================================================================================>----------------------]  84%
[=====================================================================================================================>----------------------]  85%
[======================================================================================================================>---------------------]  85%
[=======================================================================================================================>--------------------]  85%
[=======================================================================================================================>--------------------]  86%
[========================================================================================================================>-------------------]  86%
[========================================================================================================================>-------------------]  87%
[=========================================================================================================================>------------------]  87%
[==========================================================================================================================>-----------------]  88%
[===========================================================================================================================>----------------]  88%
[===========================================================================================================================>----------------]  89%
[============================================================================================================================>---------------]  89%
[============================================================================================================================>---------------]  90%
[=============================================================================================================================>--------------]  90%
[==============================================================================================================================>-------------]  90%
[==============================================================================================================================>-------------]  91%
[===============================================================================================================================>------------]  91%
[===============================================================================================================================>------------]  92%
[================================================================================================================================>-----------]  92%
[=================================================================================================================================>----------]  93%
[==================================================================================================================================>---------]  93%
[==================================================================================================================================>---------]  94%
[===================================================================================================================================>--------]  94%
[===================================================================================================================================>--------]  95%
[====================================================================================================================================>-------]  95%
[=====================================================================================================================================>------]  95%
[=====================================================================================================================================>------]  96%
[======================================================================================================================================>-----]  96%
[======================================================================================================================================>-----]  97%
[=======================================================================================================================================>----]  97%
[========================================================================================================================================>---]  98%
[=========================================================================================================================================>--]  98%
[=========================================================================================================================================>--]  99%
[==========================================================================================================================================>-]  99%
[==========================================================================================================================================>-] 100%
[===========================================================================================================================================>] 100%
[============================================================================================================================================] 100%
                                                                                                                                                   

Running diff_rhyth() on seedlist of size: 10783
NOTE: Using batches in regression.

[>-------------------------------------------------------------------------------------------------------------------------------------------]   1%
[=>------------------------------------------------------------------------------------------------------------------------------------------]   1%
[=>------------------------------------------------------------------------------------------------------------------------------------------]   2%
[==>-----------------------------------------------------------------------------------------------------------------------------------------]   2%
[===>----------------------------------------------------------------------------------------------------------------------------------------]   3%
[====>---------------------------------------------------------------------------------------------------------------------------------------]   3%
[====>---------------------------------------------------------------------------------------------------------------------------------------]   4%
[=====>--------------------------------------------------------------------------------------------------------------------------------------]   4%
[=====>--------------------------------------------------------------------------------------------------------------------------------------]   5%
[======>-------------------------------------------------------------------------------------------------------------------------------------]   5%
[=======>------------------------------------------------------------------------------------------------------------------------------------]   5%
[=======>------------------------------------------------------------------------------------------------------------------------------------]   6%
[========>-----------------------------------------------------------------------------------------------------------------------------------]   6%
[========>-----------------------------------------------------------------------------------------------------------------------------------]   7%
[=========>----------------------------------------------------------------------------------------------------------------------------------]   7%
[==========>---------------------------------------------------------------------------------------------------------------------------------]   8%
[===========>--------------------------------------------------------------------------------------------------------------------------------]   8%
[===========>--------------------------------------------------------------------------------------------------------------------------------]   9%
[============>-------------------------------------------------------------------------------------------------------------------------------]   9%
[============>-------------------------------------------------------------------------------------------------------------------------------]  10%
[=============>------------------------------------------------------------------------------------------------------------------------------]  10%
[==============>-----------------------------------------------------------------------------------------------------------------------------]  10%
[==============>-----------------------------------------------------------------------------------------------------------------------------]  11%
[===============>----------------------------------------------------------------------------------------------------------------------------]  11%
[===============>----------------------------------------------------------------------------------------------------------------------------]  12%
[================>---------------------------------------------------------------------------------------------------------------------------]  12%
[=================>--------------------------------------------------------------------------------------------------------------------------]  13%
[==================>-------------------------------------------------------------------------------------------------------------------------]  13%
[==================>-------------------------------------------------------------------------------------------------------------------------]  14%
[===================>------------------------------------------------------------------------------------------------------------------------]  14%
[===================>------------------------------------------------------------------------------------------------------------------------]  15%
[====================>-----------------------------------------------------------------------------------------------------------------------]  15%
[=====================>----------------------------------------------------------------------------------------------------------------------]  15%
[=====================>----------------------------------------------------------------------------------------------------------------------]  16%
[======================>---------------------------------------------------------------------------------------------------------------------]  16%
[======================>---------------------------------------------------------------------------------------------------------------------]  17%
[=======================>--------------------------------------------------------------------------------------------------------------------]  17%
[========================>-------------------------------------------------------------------------------------------------------------------]  18%
[=========================>------------------------------------------------------------------------------------------------------------------]  18%
[=========================>------------------------------------------------------------------------------------------------------------------]  19%
[==========================>-----------------------------------------------------------------------------------------------------------------]  19%
[==========================>-----------------------------------------------------------------------------------------------------------------]  20%
[===========================>----------------------------------------------------------------------------------------------------------------]  20%
[============================>---------------------------------------------------------------------------------------------------------------]  20%
[============================>---------------------------------------------------------------------------------------------------------------]  21%
[=============================>--------------------------------------------------------------------------------------------------------------]  21%
[=============================>--------------------------------------------------------------------------------------------------------------]  22%
[==============================>-------------------------------------------------------------------------------------------------------------]  22%
[===============================>------------------------------------------------------------------------------------------------------------]  23%
[================================>-----------------------------------------------------------------------------------------------------------]  23%
[================================>-----------------------------------------------------------------------------------------------------------]  24%
[=================================>----------------------------------------------------------------------------------------------------------]  24%
[=================================>----------------------------------------------------------------------------------------------------------]  25%
[==================================>---------------------------------------------------------------------------------------------------------]  25%
[===================================>--------------------------------------------------------------------------------------------------------]  25%
[===================================>--------------------------------------------------------------------------------------------------------]  26%
[====================================>-------------------------------------------------------------------------------------------------------]  26%
[====================================>-------------------------------------------------------------------------------------------------------]  27%
[=====================================>------------------------------------------------------------------------------------------------------]  27%
[======================================>-----------------------------------------------------------------------------------------------------]  28%
[=======================================>----------------------------------------------------------------------------------------------------]  28%
[=======================================>----------------------------------------------------------------------------------------------------]  29%
[========================================>---------------------------------------------------------------------------------------------------]  29%
[========================================>---------------------------------------------------------------------------------------------------]  30%
[=========================================>--------------------------------------------------------------------------------------------------]  30%
[==========================================>-------------------------------------------------------------------------------------------------]  30%
[==========================================>-------------------------------------------------------------------------------------------------]  31%
[===========================================>------------------------------------------------------------------------------------------------]  31%
[===========================================>------------------------------------------------------------------------------------------------]  32%
[============================================>-----------------------------------------------------------------------------------------------]  32%
[=============================================>----------------------------------------------------------------------------------------------]  33%
[==============================================>---------------------------------------------------------------------------------------------]  33%
[==============================================>---------------------------------------------------------------------------------------------]  34%
[===============================================>--------------------------------------------------------------------------------------------]  34%
[===============================================>--------------------------------------------------------------------------------------------]  35%
[================================================>-------------------------------------------------------------------------------------------]  35%
[=================================================>------------------------------------------------------------------------------------------]  35%
[=================================================>------------------------------------------------------------------------------------------]  36%
[==================================================>-----------------------------------------------------------------------------------------]  36%
[==================================================>-----------------------------------------------------------------------------------------]  37%
[===================================================>----------------------------------------------------------------------------------------]  37%
[====================================================>---------------------------------------------------------------------------------------]  38%
[=====================================================>--------------------------------------------------------------------------------------]  38%
[=====================================================>--------------------------------------------------------------------------------------]  39%
[======================================================>-------------------------------------------------------------------------------------]  39%
[======================================================>-------------------------------------------------------------------------------------]  40%
[=======================================================>------------------------------------------------------------------------------------]  40%
[========================================================>-----------------------------------------------------------------------------------]  40%
[========================================================>-----------------------------------------------------------------------------------]  41%
[=========================================================>----------------------------------------------------------------------------------]  41%
[=========================================================>----------------------------------------------------------------------------------]  42%
[==========================================================>---------------------------------------------------------------------------------]  42%
[===========================================================>--------------------------------------------------------------------------------]  43%
[============================================================>-------------------------------------------------------------------------------]  43%
[============================================================>-------------------------------------------------------------------------------]  44%
[=============================================================>------------------------------------------------------------------------------]  44%
[=============================================================>------------------------------------------------------------------------------]  45%
[==============================================================>-----------------------------------------------------------------------------]  45%
[===============================================================>----------------------------------------------------------------------------]  45%
[===============================================================>----------------------------------------------------------------------------]  46%
[================================================================>---------------------------------------------------------------------------]  46%
[================================================================>---------------------------------------------------------------------------]  47%
[=================================================================>--------------------------------------------------------------------------]  47%
[==================================================================>-------------------------------------------------------------------------]  48%
[===================================================================>------------------------------------------------------------------------]  48%
[===================================================================>------------------------------------------------------------------------]  49%
[====================================================================>-----------------------------------------------------------------------]  49%
[====================================================================>-----------------------------------------------------------------------]  50%
[=====================================================================>----------------------------------------------------------------------]  50%
[======================================================================>---------------------------------------------------------------------]  50%
[======================================================================>---------------------------------------------------------------------]  51%
[=======================================================================>--------------------------------------------------------------------]  51%
[=======================================================================>--------------------------------------------------------------------]  52%
[========================================================================>-------------------------------------------------------------------]  52%
[=========================================================================>------------------------------------------------------------------]  53%
[==========================================================================>-----------------------------------------------------------------]  53%
[==========================================================================>-----------------------------------------------------------------]  54%
[===========================================================================>----------------------------------------------------------------]  54%
[===========================================================================>----------------------------------------------------------------]  55%
[============================================================================>---------------------------------------------------------------]  55%
[=============================================================================>--------------------------------------------------------------]  55%
[=============================================================================>--------------------------------------------------------------]  56%
[==============================================================================>-------------------------------------------------------------]  56%
[==============================================================================>-------------------------------------------------------------]  57%
[===============================================================================>------------------------------------------------------------]  57%
[================================================================================>-----------------------------------------------------------]  58%
[=================================================================================>----------------------------------------------------------]  58%
[=================================================================================>----------------------------------------------------------]  59%
[==================================================================================>---------------------------------------------------------]  59%
[==================================================================================>---------------------------------------------------------]  60%
[===================================================================================>--------------------------------------------------------]  60%
[====================================================================================>-------------------------------------------------------]  60%
[====================================================================================>-------------------------------------------------------]  61%
[=====================================================================================>------------------------------------------------------]  61%
[=====================================================================================>------------------------------------------------------]  62%
[======================================================================================>-----------------------------------------------------]  62%
[=======================================================================================>----------------------------------------------------]  63%
[========================================================================================>---------------------------------------------------]  63%
[========================================================================================>---------------------------------------------------]  64%
[=========================================================================================>--------------------------------------------------]  64%
[=========================================================================================>--------------------------------------------------]  65%
[==========================================================================================>-------------------------------------------------]  65%
[===========================================================================================>------------------------------------------------]  65%
[===========================================================================================>------------------------------------------------]  66%
[============================================================================================>-----------------------------------------------]  66%
[============================================================================================>-----------------------------------------------]  67%
[=============================================================================================>----------------------------------------------]  67%
[==============================================================================================>---------------------------------------------]  68%
[===============================================================================================>--------------------------------------------]  68%
[===============================================================================================>--------------------------------------------]  69%
[================================================================================================>-------------------------------------------]  69%
[================================================================================================>-------------------------------------------]  70%
[=================================================================================================>------------------------------------------]  70%
[==================================================================================================>-----------------------------------------]  70%
[==================================================================================================>-----------------------------------------]  71%
[===================================================================================================>----------------------------------------]  71%
[===================================================================================================>----------------------------------------]  72%
[====================================================================================================>---------------------------------------]  72%
[=====================================================================================================>--------------------------------------]  73%
[======================================================================================================>-------------------------------------]  73%
[======================================================================================================>-------------------------------------]  74%
[=======================================================================================================>------------------------------------]  74%
[=======================================================================================================>------------------------------------]  75%
[========================================================================================================>-----------------------------------]  75%
[=========================================================================================================>----------------------------------]  75%
[=========================================================================================================>----------------------------------]  76%
[==========================================================================================================>---------------------------------]  76%
[==========================================================================================================>---------------------------------]  77%
[===========================================================================================================>--------------------------------]  77%
[============================================================================================================>-------------------------------]  78%
[=============================================================================================================>------------------------------]  78%
[=============================================================================================================>------------------------------]  79%
[==============================================================================================================>-----------------------------]  79%
[==============================================================================================================>-----------------------------]  80%
[===============================================================================================================>----------------------------]  80%
[================================================================================================================>---------------------------]  80%
[================================================================================================================>---------------------------]  81%
[=================================================================================================================>--------------------------]  81%
[=================================================================================================================>--------------------------]  82%
[==================================================================================================================>-------------------------]  82%
[===================================================================================================================>------------------------]  83%
[====================================================================================================================>-----------------------]  83%
[====================================================================================================================>-----------------------]  84%
[=====================================================================================================================>----------------------]  84%
[=====================================================================================================================>----------------------]  85%
[======================================================================================================================>---------------------]  85%
[=======================================================================================================================>--------------------]  85%
[=======================================================================================================================>--------------------]  86%
[========================================================================================================================>-------------------]  86%
[========================================================================================================================>-------------------]  87%
[=========================================================================================================================>------------------]  87%
[==========================================================================================================================>-----------------]  88%
[===========================================================================================================================>----------------]  88%
[===========================================================================================================================>----------------]  89%
[============================================================================================================================>---------------]  89%
[============================================================================================================================>---------------]  90%
[=============================================================================================================================>--------------]  90%
[==============================================================================================================================>-------------]  90%
[==============================================================================================================================>-------------]  91%
[===============================================================================================================================>------------]  91%
[===============================================================================================================================>------------]  92%
[================================================================================================================================>-----------]  92%
[=================================================================================================================================>----------]  93%
[==================================================================================================================================>---------]  93%
[==================================================================================================================================>---------]  94%
[===================================================================================================================================>--------]  94%
[===================================================================================================================================>--------]  95%
[====================================================================================================================================>-------]  95%
[=====================================================================================================================================>------]  95%
[=====================================================================================================================================>------]  96%
[======================================================================================================================================>-----]  96%
[======================================================================================================================================>-----]  97%
[=======================================================================================================================================>----]  97%
[========================================================================================================================================>---]  98%
[=========================================================================================================================================>--]  98%
[=========================================================================================================================================>--]  99%
[==========================================================================================================================================>-]  99%
[==========================================================================================================================================>-] 100%
[===========================================================================================================================================>] 100%
[============================================================================================================================================] 100%
                                                                                                                                                   

Plotting

With custom regression, I have to remake clock face plots, genes plots, etc.

setwd(path_to_cyclops_ordering)
Warning: The working directory was changed to /Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit inside a notebook chunk. The working directory will be reset when the chunk is finished running. Use the knitr root.dir option in the setup chunk to change the working directory for notebook chunks.
if (!(dir.exists("downstream_output/plots"))){
  dir.create("downstream_output/plots")
}

#for cycling_in_CTL.csv
plot_clock_face(paste0(path_to_cyclops_ordering, "/downstream_output/plots/CTL_clock_face_plot.png"),"downstream_output/cosinor_results_CTL.csv", mouse_data = mouse_data, BHQ_cutoff = isCyclingBHQCutoff, amp_ratio_cutoff = 0)
[1] "ARNTL set to 0 radians"
NULL
setwd(path_to_cyclops_ordering)

#for cycling_in_AD.csv
plot_clock_face(paste0(path_to_cyclops_ordering, "/downstream_output/plots/AD_clock_face_plot.png"),"downstream_output/cosinor_results_AD.csv",mouse_data = mouse_data, BHQ_cutoff = isCyclingBHQCutoff, amp_ratio_cutoff = 0)
[1] "ARNTL set to 0 radians"
NULL
#plot cyclops histogram
plot_subject_histogram(path_to_cyclops_ordering, "cond_0")
null device 
          1 
plot_subject_histogram(path_to_cyclops_ordering, "cond_1")
null device 
          1 
#Gene Tracings
plot_core_clock_genes(path_to_tmm_file, path_to_cyclops_ordering, useBatch = useBatch, percentile = percentile, split_cond_plots = T)
[1] "NOTE: Using batches in plotting"

Plot Other Genes

plot_me = c("RPL18","RPS7","RPL22","RPS18","RPL8","RPL10A","RPS13","RPL14","RPS15A","RPS12","RPL36AL","RPL34","RPL30","RPL9","RPL21","RPS28","RPSA","RPL26","RPL27A","FAU","RPS3A","RPL7A","RPL17","RPS27A","RPS16","RPL15","RPLP1","RPL7","RPL10","RPL23A","RPL18A","RPL3","RPL12","RPL39","RPL37A","RPL37","RPS10","RPL13","RPL35A","RPS23","RPL13A","RPL38","RPS29","RPL24","RPL41","RPL35","RPS26","RPS2","RPS9","RPL32","RPL6","RPS21","RPS24","RPS25","RPL36","RPL19","RPS27","RPL29","RPS15","RPS3","RPS5","RPS8","RPL27","RSL24D1","RPS6","RPL11")

tmm = read_csv(path_to_tmm_file, show_col_types = FALSE)
  cyc_pred_file = list.files(path = paste0(path_to_cyclops_ordering, "/Fits/"), pattern = '*Fit_Output_*')
cyc_pred = read_csv(paste(path_to_cyclops_ordering, "Fits", cyc_pred_file[1], sep = '/'), show_col_types = FALSE)
plot_gene_trace(cyc_pred, tmm, plot_me , percentile = percentile, useBatch = useBatch)
[1] "NOTE: Using batches in plotting"

EnrichR

Call EnrichR API with python, supplying various lists from cycling and differential rhythm analysis.

path = getwd()
setwd(paste0(path_to_cyclops_ordering, "downstream_output"))
Warning: The working directory was changed to /Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit/downstream_output inside a notebook chunk. The working directory will be reset when the chunk is finished running. Use the knitr root.dir option in the setup chunk to change the working directory for notebook chunks.
path_no_space = gsub(" ", "\\ ", path, fixed = TRUE)
#create list of filenames from the enrichR_files subfolder
gene_lists = list.files("./enrichR_files/", pattern = "^(CTL_|AD_|diff_mesor)")

#run system command to call python script that uses enrichR API
sapply(gene_lists, function(x){ 
  system(paste0("python3 ", path_no_space, "/Python_EnrichR_for_AnalysisPipeline.py -g \"enrichR_files/" ,x, "\" -b \"enrichR_files/EnrichR_background.csv\""))})
{'lib': ['KEGG_2021_Human', 'GO_Biological_Process_2023', 'MSigDB_Hallmark_2020', 'WikiPathway_2021_Human', 'TRRUST_Transcription_Factors_2019', 'GO_Cellular_Component_2023', 'TRANSFAC_and_JASPAR_PWMs', 'Transcription_Factor_PPIs', 'TargetScan_microRNA_2017', 'miRTarBase_2017'], 'genes_list_file': 'enrichR_files/AD_CTL_mthd2_cyclers_AR25BHQ1.csv', 'background_list_file': 'enrichR_files/EnrichR_background.csv'}
Please select file containing list of genes (Symbols)
Path to genes list: enrichR_files/AD_CTL_mthd2_cyclers_AR25BHQ1.csv

Please select file containing list of BACKGROUND genes (Symbols)
Path to Background list: enrichR_files/EnrichR_background.csv
Background DataFrame (nrow: 4337):
           Ensembl Gene_Symbols
0  ENSG00000244405         ETV5
1  ENSG00000132326         PER2
2  ENSG00000171314        PGAM1
3  ENSG00000188766       SPRED3
4  ENSG00000100335        MIEF1
Background DataFrame (nrow: 13310):
           Ensembl Gene_Symbols
0  ENSG00000188976        NOC2L
1  ENSG00000187961       KLHL17
2  ENSG00000188290         HES4
3  ENSG00000187608        ISG15
4  ENSG00000188157         AGRN
Directory "EnrichR_results/EnrichR_ADCTLmthd2cyclersAR25BHQ1_BCKGRND-EnrichRbackground" already exists. Skipping.
{'lib': ['KEGG_2021_Human', 'GO_Biological_Process_2023', 'MSigDB_Hallmark_2020', 'WikiPathway_2021_Human', 'TRRUST_Transcription_Factors_2019', 'GO_Cellular_Component_2023', 'TRANSFAC_and_JASPAR_PWMs', 'Transcription_Factor_PPIs', 'TargetScan_microRNA_2017', 'miRTarBase_2017'], 'genes_list_file': 'enrichR_files/AD_cyclers_AR25BHQ1.csv', 'background_list_file': 'enrichR_files/EnrichR_background.csv'}
Please select file containing list of genes (Symbols)
Path to genes list: enrichR_files/AD_cyclers_AR25BHQ1.csv

Please select file containing list of BACKGROUND genes (Symbols)
Path to Background list: enrichR_files/EnrichR_background.csv
Background DataFrame (nrow: 3793):
           Ensembl Gene_Symbols
0  ENSG00000244405         ETV5
1  ENSG00000132326         PER2
2  ENSG00000171314        PGAM1
3  ENSG00000130540      SULT4A1
4  ENSG00000115310         RTN4
Background DataFrame (nrow: 13310):
           Ensembl Gene_Symbols
0  ENSG00000188976        NOC2L
1  ENSG00000187961       KLHL17
2  ENSG00000188290         HES4
3  ENSG00000187608        ISG15
4  ENSG00000188157         AGRN
Directory "EnrichR_results/EnrichR_ADcyclersAR25BHQ1_BCKGRND-EnrichRbackground" already exists. Skipping.
{'lib': ['KEGG_2021_Human', 'GO_Biological_Process_2023', 'MSigDB_Hallmark_2020', 'WikiPathway_2021_Human', 'TRRUST_Transcription_Factors_2019', 'GO_Cellular_Component_2023', 'TRANSFAC_and_JASPAR_PWMs', 'Transcription_Factor_PPIs', 'TargetScan_microRNA_2017', 'miRTarBase_2017'], 'genes_list_file': 'enrichR_files/AD_cyclers_AR33BHQ1.csv', 'background_list_file': 'enrichR_files/EnrichR_background.csv'}
Please select file containing list of genes (Symbols)
Path to genes list: enrichR_files/AD_cyclers_AR33BHQ1.csv

Please select file containing list of BACKGROUND genes (Symbols)
Path to Background list: enrichR_files/EnrichR_background.csv
Background DataFrame (nrow: 2036):
           Ensembl Gene_Symbols
0  ENSG00000244405         ETV5
1  ENSG00000132326         PER2
2  ENSG00000171314        PGAM1
3  ENSG00000130540      SULT4A1
4  ENSG00000115310         RTN4
Background DataFrame (nrow: 13310):
           Ensembl Gene_Symbols
0  ENSG00000188976        NOC2L
1  ENSG00000187961       KLHL17
2  ENSG00000188290         HES4
3  ENSG00000187608        ISG15
4  ENSG00000188157         AGRN
Directory "EnrichR_results/EnrichR_ADcyclersAR33BHQ1_BCKGRND-EnrichRbackground" already exists. Skipping.
{'lib': ['KEGG_2021_Human', 'GO_Biological_Process_2023', 'MSigDB_Hallmark_2020', 'WikiPathway_2021_Human', 'TRRUST_Transcription_Factors_2019', 'GO_Cellular_Component_2023', 'TRANSFAC_and_JASPAR_PWMs', 'Transcription_Factor_PPIs', 'TargetScan_microRNA_2017', 'miRTarBase_2017'], 'genes_list_file': 'enrichR_files/CTL_cyclers_AR25BHQ1.csv', 'background_list_file': 'enrichR_files/EnrichR_background.csv'}
Please select file containing list of genes (Symbols)
Path to genes list: enrichR_files/CTL_cyclers_AR25BHQ1.csv

Please select file containing list of BACKGROUND genes (Symbols)
Path to Background list: enrichR_files/EnrichR_background.csv
Background DataFrame (nrow: 3501):
           Ensembl Gene_Symbols
0  ENSG00000244405         ETV5
1  ENSG00000132326         PER2
2  ENSG00000188766       SPRED3
3  ENSG00000171314        PGAM1
4  ENSG00000124785         NRN1
Background DataFrame (nrow: 13310):
           Ensembl Gene_Symbols
0  ENSG00000188976        NOC2L
1  ENSG00000187961       KLHL17
2  ENSG00000188290         HES4
3  ENSG00000187608        ISG15
4  ENSG00000188157         AGRN
Directory "EnrichR_results/EnrichR_CTLcyclersAR25BHQ1_BCKGRND-EnrichRbackground" already exists. Skipping.
{'lib': ['KEGG_2021_Human', 'GO_Biological_Process_2023', 'MSigDB_Hallmark_2020', 'WikiPathway_2021_Human', 'TRRUST_Transcription_Factors_2019', 'GO_Cellular_Component_2023', 'TRANSFAC_and_JASPAR_PWMs', 'Transcription_Factor_PPIs', 'TargetScan_microRNA_2017', 'miRTarBase_2017'], 'genes_list_file': 'enrichR_files/CTL_cyclers_AR33BHQ1.csv', 'background_list_file': 'enrichR_files/EnrichR_background.csv'}
Please select file containing list of genes (Symbols)
Path to genes list: enrichR_files/CTL_cyclers_AR33BHQ1.csv

Please select file containing list of BACKGROUND genes (Symbols)
Path to Background list: enrichR_files/EnrichR_background.csv
Background DataFrame (nrow: 1678):
           Ensembl Gene_Symbols
0  ENSG00000244405         ETV5
1  ENSG00000132326         PER2
2  ENSG00000188766       SPRED3
3  ENSG00000171314        PGAM1
4  ENSG00000124785         NRN1
Background DataFrame (nrow: 13310):
           Ensembl Gene_Symbols
0  ENSG00000188976        NOC2L
1  ENSG00000187961       KLHL17
2  ENSG00000188290         HES4
3  ENSG00000187608        ISG15
4  ENSG00000188157         AGRN
Directory "EnrichR_results/EnrichR_CTLcyclersAR33BHQ1_BCKGRND-EnrichRbackground" already exists. Skipping.
{'lib': ['KEGG_2021_Human', 'GO_Biological_Process_2023', 'MSigDB_Hallmark_2020', 'WikiPathway_2021_Human', 'TRRUST_Transcription_Factors_2019', 'GO_Cellular_Component_2023', 'TRANSFAC_and_JASPAR_PWMs', 'Transcription_Factor_PPIs', 'TargetScan_microRNA_2017', 'miRTarBase_2017'], 'genes_list_file': 'enrichR_files/diff_mesor_all_genes_BHQ05.csv', 'background_list_file': 'enrichR_files/EnrichR_background.csv'}
Please select file containing list of genes (Symbols)
Path to genes list: enrichR_files/diff_mesor_all_genes_BHQ05.csv

Please select file containing list of BACKGROUND genes (Symbols)
Path to Background list: enrichR_files/EnrichR_background.csv
Background DataFrame (nrow: 6464):
           Ensembl Gene_Symbols
0  ENSG00000188976        NOC2L
1  ENSG00000187961       KLHL17
2  ENSG00000188290         HES4
3  ENSG00000187608        ISG15
4  ENSG00000078808         SDF4
Background DataFrame (nrow: 13310):
           Ensembl Gene_Symbols
0  ENSG00000188976        NOC2L
1  ENSG00000187961       KLHL17
2  ENSG00000188290         HES4
3  ENSG00000187608        ISG15
4  ENSG00000188157         AGRN
Directory "EnrichR_results/EnrichR_diffmesorallgenesBHQ05_BCKGRND-EnrichRbackground" already exists. Skipping.
AD_CTL_mthd2_cyclers_AR25BHQ1.csv           AD_cyclers_AR25BHQ1.csv           AD_cyclers_AR33BHQ1.csv          CTL_cyclers_AR25BHQ1.csv 
                                0                                 0                                 0                                 0 
         CTL_cyclers_AR33BHQ1.csv    diff_mesor_all_genes_BHQ05.csv 
                                0                                 0 
isCyclingBHQCutoff_str = str_extract(as.character(isCyclingBHQCutoff), "(?<=\\.)\\d+")

setwd("enrichR_files/")
all_DR_gene_lists = list.files( pattern = "^(DR_)")
method2_files = grepl("mthd2", all_DR_gene_lists)
DR_gene_lists_method1 = all_DR_gene_lists[!method2_files]
DR_backgrounds_method1 = paste0("diff_rhythms_CyclingBHQ", isCyclingBHQCutoff_str, "AmpRatio", str_replace(str_extract(DR_gene_lists_method1, pattern = "AR\\d+"), "AR", ""), ".csv")
DR_gene_lists_method2 = all_DR_gene_lists[method2_files]
DR_backgrounds_method2 = paste0("diff_rhythms_method2_CyclingBHQ", isCyclingBHQCutoff_str, "AmpRatio", str_replace(str_extract(DR_gene_lists_method2, pattern = "AR\\d+"), "AR", ""), ".csv")
DR_gene_lists = c(DR_gene_lists_method1, DR_gene_lists_method2)
DR_backgrounds = c(DR_backgrounds_method1, DR_backgrounds_method2)

# repeat with differential gene lists
setwd(paste0(path_to_cyclops_ordering, "downstream_output"))
mapply(function(x, y){ 
  system(paste0("python3 ", path_no_space, "/Python_EnrichR_for_AnalysisPipeline.py -g \"enrichR_files/" ,x, "\" -b \"",y , "\""))},
  DR_gene_lists, 
  DR_backgrounds)
{'lib': ['KEGG_2021_Human', 'GO_Biological_Process_2023', 'MSigDB_Hallmark_2020', 'WikiPathway_2021_Human', 'TRRUST_Transcription_Factors_2019', 'GO_Cellular_Component_2023', 'TRANSFAC_and_JASPAR_PWMs', 'Transcription_Factor_PPIs', 'TargetScan_microRNA_2017', 'miRTarBase_2017'], 'genes_list_file': 'enrichR_files/DR_cyclers_CyclingBHQ1AR1_DRBHQ2.csv', 'background_list_file': 'diff_rhythms_CyclingBHQ1AmpRatio1.csv'}
Please select file containing list of genes (Symbols)
Path to genes list: enrichR_files/DR_cyclers_CyclingBHQ1AR1_DRBHQ2.csv

Please select file containing list of BACKGROUND genes (Symbols)
Path to Background list: diff_rhythms_CyclingBHQ1AmpRatio1.csv
Background DataFrame (nrow: 54):
           Ensembl Gene_Symbols  ...      Bonf  Log_AD_CTL_ampRatio
0  ENSG00000161281       COX7A1  ...  0.020616            -0.704488
1  ENSG00000211450      SELENOH  ...  0.068573            -0.454206
2  ENSG00000164405        UQCRQ  ...  1.000000            -0.565959
3  ENSG00000019582         CD74  ...  1.000000             1.234542
4  ENSG00000177556        ATOX1  ...  1.000000            -0.541007

[5 rows x 10 columns]
Background DataFrame (nrow: 10528):
           Ensembl Gene_Symbols         p  ...       BHQ  Bonf  Log_AD_CTL_ampRatio
0  ENSG00000188976        NOC2L  0.418157  ...  0.787534   1.0             0.333543
1  ENSG00000187961       KLHL17  0.021538  ...  0.402040   1.0             0.371275
2  ENSG00000188290         HES4  0.509391  ...  0.824935   1.0            -0.120408
3  ENSG00000187608        ISG15  0.508770  ...  0.824906   1.0            -0.294342
4  ENSG00000131591     C1orf159  0.782443  ...  0.930163   1.0             0.084900

[5 rows x 10 columns]
Directory "EnrichR_results/EnrichR_DRcyclersCyclingBHQ1AR1DRBHQ2_BCKGRND-diffrhythmsCyclingBHQ1AmpRatio1" already exists. Skipping.
{'lib': ['KEGG_2021_Human', 'GO_Biological_Process_2023', 'MSigDB_Hallmark_2020', 'WikiPathway_2021_Human', 'TRRUST_Transcription_Factors_2019', 'GO_Cellular_Component_2023', 'TRANSFAC_and_JASPAR_PWMs', 'Transcription_Factor_PPIs', 'TargetScan_microRNA_2017', 'miRTarBase_2017'], 'genes_list_file': 'enrichR_files/DR_cyclers_CyclingBHQ1AR25_DRBHQ2.csv', 'background_list_file': 'diff_rhythms_CyclingBHQ1AmpRatio25.csv'}
Please select file containing list of genes (Symbols)
Path to genes list: enrichR_files/DR_cyclers_CyclingBHQ1AR25_DRBHQ2.csv

Please select file containing list of BACKGROUND genes (Symbols)
Path to Background list: diff_rhythms_CyclingBHQ1AmpRatio25.csv
Background DataFrame (nrow: 230):
           Ensembl Gene_Symbols  ...      Bonf  Log_AD_CTL_ampRatio
0  ENSG00000161281       COX7A1  ...  0.008463            -0.704488
1  ENSG00000211450      SELENOH  ...  0.028151            -0.454206
2  ENSG00000164405        UQCRQ  ...  0.441316            -0.565959
3  ENSG00000019582         CD74  ...  0.580904             1.234542
4  ENSG00000177556        ATOX1  ...  0.633190            -0.541007

[5 rows x 10 columns]
Background DataFrame (nrow: 4322):
           Ensembl Gene_Symbols         p  ...       BHQ  Bonf  Log_AD_CTL_ampRatio
0  ENSG00000188290         HES4  0.509391  ...  0.772823   1.0            -0.120408
1  ENSG00000187608        ISG15  0.508770  ...  0.772823   1.0            -0.294342
2  ENSG00000131591     C1orf159  0.782443  ...  0.906011   1.0             0.084900
3  ENSG00000176022      B3GALT6  0.360561  ...  0.690359   1.0            -0.083917
4  ENSG00000169972        PUSL1  0.216335  ...  0.599048   1.0            -0.298418

[5 rows x 10 columns]
Directory "EnrichR_results/EnrichR_DRcyclersCyclingBHQ1AR25DRBHQ2_BCKGRND-diffrhythmsCyclingBHQ1AmpRatio25" already exists. Skipping.
{'lib': ['KEGG_2021_Human', 'GO_Biological_Process_2023', 'MSigDB_Hallmark_2020', 'WikiPathway_2021_Human', 'TRRUST_Transcription_Factors_2019', 'GO_Cellular_Component_2023', 'TRANSFAC_and_JASPAR_PWMs', 'Transcription_Factor_PPIs', 'TargetScan_microRNA_2017', 'miRTarBase_2017'], 'genes_list_file': 'enrichR_files/DR_cyclers_CyclingBHQ1AR33_DRBHQ2.csv', 'background_list_file': 'diff_rhythms_CyclingBHQ1AmpRatio33.csv'}
Please select file containing list of genes (Symbols)
Path to genes list: enrichR_files/DR_cyclers_CyclingBHQ1AR33_DRBHQ2.csv

Please select file containing list of BACKGROUND genes (Symbols)
Path to Background list: diff_rhythms_CyclingBHQ1AmpRatio33.csv
Background DataFrame (nrow: 316):
           Ensembl Gene_Symbols  ...      Bonf  Log_AD_CTL_ampRatio
0  ENSG00000161281       COX7A1  ...  0.004578            -0.704488
1  ENSG00000211450      SELENOH  ...  0.015228            -0.454206
2  ENSG00000164405        UQCRQ  ...  0.238731            -0.565959
3  ENSG00000124614        RPS10  ...  0.260744            -0.525353
4  ENSG00000137168        PPIL1  ...  0.216483            -0.628644

[5 rows x 10 columns]
Background DataFrame (nrow: 2338):
           Ensembl Gene_Symbols         p  ...       BHQ  Bonf  Log_AD_CTL_ampRatio
0  ENSG00000188290         HES4  0.509391  ...  0.732026   1.0            -0.120408
1  ENSG00000176022      B3GALT6  0.360561  ...  0.629568   1.0            -0.083917
2  ENSG00000175756     AURKAIP1  0.016927  ...  0.155197   1.0            -0.334959
3  ENSG00000235098      ANKRD65  0.025834  ...  0.195471   1.0            -0.320663
4  ENSG00000160075        SSU72  0.160727  ...  0.451659   1.0             0.010688

[5 rows x 10 columns]
Directory "EnrichR_results/EnrichR_DRcyclersCyclingBHQ1AR33DRBHQ2_BCKGRND-diffrhythmsCyclingBHQ1AmpRatio33" already exists. Skipping.
{'lib': ['KEGG_2021_Human', 'GO_Biological_Process_2023', 'MSigDB_Hallmark_2020', 'WikiPathway_2021_Human', 'TRRUST_Transcription_Factors_2019', 'GO_Cellular_Component_2023', 'TRANSFAC_and_JASPAR_PWMs', 'Transcription_Factor_PPIs', 'TargetScan_microRNA_2017', 'miRTarBase_2017'], 'genes_list_file': 'enrichR_files/DR_gainAmpAD_CyclingBHQ1AR1_DRBHQ2.csv', 'background_list_file': 'diff_rhythms_CyclingBHQ1AmpRatio1.csv'}
Please select file containing list of genes (Symbols)
Path to genes list: enrichR_files/DR_gainAmpAD_CyclingBHQ1AR1_DRBHQ2.csv

Please select file containing list of BACKGROUND genes (Symbols)
Path to Background list: diff_rhythms_CyclingBHQ1AmpRatio1.csv
Background DataFrame (nrow: 9):
           Ensembl Gene_Symbols
0  ENSG00000143303       RRNAD1
1  ENSG00000162843        WDR64
2  ENSG00000226479     TMEM185B
3  ENSG00000019582         CD74
4  ENSG00000196422      PPP1R26
Background DataFrame (nrow: 10528):
           Ensembl Gene_Symbols         p  ...       BHQ  Bonf  Log_AD_CTL_ampRatio
0  ENSG00000188976        NOC2L  0.418157  ...  0.787534   1.0             0.333543
1  ENSG00000187961       KLHL17  0.021538  ...  0.402040   1.0             0.371275
2  ENSG00000188290         HES4  0.509391  ...  0.824935   1.0            -0.120408
3  ENSG00000187608        ISG15  0.508770  ...  0.824906   1.0            -0.294342
4  ENSG00000131591     C1orf159  0.782443  ...  0.930163   1.0             0.084900

[5 rows x 10 columns]
Directory "EnrichR_results/EnrichR_DRgainAmpADCyclingBHQ1AR1DRBHQ2_BCKGRND-diffrhythmsCyclingBHQ1AmpRatio1" already exists. Skipping.
{'lib': ['KEGG_2021_Human', 'GO_Biological_Process_2023', 'MSigDB_Hallmark_2020', 'WikiPathway_2021_Human', 'TRRUST_Transcription_Factors_2019', 'GO_Cellular_Component_2023', 'TRANSFAC_and_JASPAR_PWMs', 'Transcription_Factor_PPIs', 'TargetScan_microRNA_2017', 'miRTarBase_2017'], 'genes_list_file': 'enrichR_files/DR_gainAmpAD_CyclingBHQ1AR25_DRBHQ2.csv', 'background_list_file': 'diff_rhythms_CyclingBHQ1AmpRatio25.csv'}
Please select file containing list of genes (Symbols)
Path to genes list: enrichR_files/DR_gainAmpAD_CyclingBHQ1AR25_DRBHQ2.csv

Please select file containing list of BACKGROUND genes (Symbols)
Path to Background list: diff_rhythms_CyclingBHQ1AmpRatio25.csv
Background DataFrame (nrow: 46):
           Ensembl Gene_Symbols
0  ENSG00000133063        CHIT1
1  ENSG00000162843        WDR64
2  ENSG00000226479     TMEM185B
3  ENSG00000144395      CCDC150
4  ENSG00000114779      ABHD14B
Background DataFrame (nrow: 4322):
           Ensembl Gene_Symbols         p  ...       BHQ  Bonf  Log_AD_CTL_ampRatio
0  ENSG00000188290         HES4  0.509391  ...  0.772823   1.0            -0.120408
1  ENSG00000187608        ISG15  0.508770  ...  0.772823   1.0            -0.294342
2  ENSG00000131591     C1orf159  0.782443  ...  0.906011   1.0             0.084900
3  ENSG00000176022      B3GALT6  0.360561  ...  0.690359   1.0            -0.083917
4  ENSG00000169972        PUSL1  0.216335  ...  0.599048   1.0            -0.298418

[5 rows x 10 columns]
Directory "EnrichR_results/EnrichR_DRgainAmpADCyclingBHQ1AR25DRBHQ2_BCKGRND-diffrhythmsCyclingBHQ1AmpRatio25" already exists. Skipping.
{'lib': ['KEGG_2021_Human', 'GO_Biological_Process_2023', 'MSigDB_Hallmark_2020', 'WikiPathway_2021_Human', 'TRRUST_Transcription_Factors_2019', 'GO_Cellular_Component_2023', 'TRANSFAC_and_JASPAR_PWMs', 'Transcription_Factor_PPIs', 'TargetScan_microRNA_2017', 'miRTarBase_2017'], 'genes_list_file': 'enrichR_files/DR_gainAmpAD_CyclingBHQ1AR33_DRBHQ2.csv', 'background_list_file': 'diff_rhythms_CyclingBHQ1AmpRatio33.csv'}
Please select file containing list of genes (Symbols)
Path to genes list: enrichR_files/DR_gainAmpAD_CyclingBHQ1AR33_DRBHQ2.csv

Please select file containing list of BACKGROUND genes (Symbols)
Path to Background list: diff_rhythms_CyclingBHQ1AmpRatio33.csv
Background DataFrame (nrow: 62):
           Ensembl Gene_Symbols
0  ENSG00000222046       DCDC2B
1  ENSG00000228063   LYPLAL1-DT
2  ENSG00000162843        WDR64
3  ENSG00000171848         RRM2
4  ENSG00000226479     TMEM185B
Background DataFrame (nrow: 2338):
           Ensembl Gene_Symbols         p  ...       BHQ  Bonf  Log_AD_CTL_ampRatio
0  ENSG00000188290         HES4  0.509391  ...  0.732026   1.0            -0.120408
1  ENSG00000176022      B3GALT6  0.360561  ...  0.629568   1.0            -0.083917
2  ENSG00000175756     AURKAIP1  0.016927  ...  0.155197   1.0            -0.334959
3  ENSG00000235098      ANKRD65  0.025834  ...  0.195471   1.0            -0.320663
4  ENSG00000160075        SSU72  0.160727  ...  0.451659   1.0             0.010688

[5 rows x 10 columns]
Directory "EnrichR_results/EnrichR_DRgainAmpADCyclingBHQ1AR33DRBHQ2_BCKGRND-diffrhythmsCyclingBHQ1AmpRatio33" already exists. Skipping.
{'lib': ['KEGG_2021_Human', 'GO_Biological_Process_2023', 'MSigDB_Hallmark_2020', 'WikiPathway_2021_Human', 'TRRUST_Transcription_Factors_2019', 'GO_Cellular_Component_2023', 'TRANSFAC_and_JASPAR_PWMs', 'Transcription_Factor_PPIs', 'TargetScan_microRNA_2017', 'miRTarBase_2017'], 'genes_list_file': 'enrichR_files/DR_lostAmpAD_CyclingBHQ1AR1_DRBHQ2.csv', 'background_list_file': 'diff_rhythms_CyclingBHQ1AmpRatio1.csv'}
Please select file containing list of genes (Symbols)
Path to genes list: enrichR_files/DR_lostAmpAD_CyclingBHQ1AR1_DRBHQ2.csv

Please select file containing list of BACKGROUND genes (Symbols)
Path to Background list: diff_rhythms_CyclingBHQ1AmpRatio1.csv
Background DataFrame (nrow: 45):
           Ensembl Gene_Symbols
0  ENSG00000178796       RIIAD1
1  ENSG00000177954        RPS27
2  ENSG00000188010        MORN2
3  ENSG00000143947       RPS27A
4  ENSG00000135617       PRADC1
Background DataFrame (nrow: 10528):
           Ensembl Gene_Symbols         p  ...       BHQ  Bonf  Log_AD_CTL_ampRatio
0  ENSG00000188976        NOC2L  0.418157  ...  0.787534   1.0             0.333543
1  ENSG00000187961       KLHL17  0.021538  ...  0.402040   1.0             0.371275
2  ENSG00000188290         HES4  0.509391  ...  0.824935   1.0            -0.120408
3  ENSG00000187608        ISG15  0.508770  ...  0.824906   1.0            -0.294342
4  ENSG00000131591     C1orf159  0.782443  ...  0.930163   1.0             0.084900

[5 rows x 10 columns]
Directory "EnrichR_results/EnrichR_DRlostAmpADCyclingBHQ1AR1DRBHQ2_BCKGRND-diffrhythmsCyclingBHQ1AmpRatio1" created.
Saved KEGG_2021_Human results to file
Saved GO_Biological_Process_2023 results to file
Saved MSigDB_Hallmark_2020 results to file
Saved WikiPathway_2021_Human results to file
{'lib': ['KEGG_2021_Human', 'GO_Biological_Process_2023', 'MSigDB_Hallmark_2020', 'WikiPathway_2021_Human', 'TRRUST_Transcription_Factors_2019', 'GO_Cellular_Component_2023', 'TRANSFAC_and_JASPAR_PWMs', 'Transcription_Factor_PPIs', 'TargetScan_microRNA_2017', 'miRTarBase_2017'], 'genes_list_file': 'enrichR_files/DR_lostAmpAD_CyclingBHQ1AR25_DRBHQ2.csv', 'background_list_file': 'diff_rhythms_CyclingBHQ1AmpRatio25.csv'}
Please select file containing list of genes (Symbols)
Path to genes list: enrichR_files/DR_lostAmpAD_CyclingBHQ1AR25_DRBHQ2.csv

Please select file containing list of BACKGROUND genes (Symbols)
Path to Background list: diff_rhythms_CyclingBHQ1AmpRatio25.csv
Background DataFrame (nrow: 184):
           Ensembl Gene_Symbols
0  ENSG00000142676        RPL11
1  ENSG00000130770      ATP5IF1
2  ENSG00000142937         RPS8
3  ENSG00000173846         PLK3
4  ENSG00000173660        UQCRH
Background DataFrame (nrow: 4322):
           Ensembl Gene_Symbols         p  ...       BHQ  Bonf  Log_AD_CTL_ampRatio
0  ENSG00000188290         HES4  0.509391  ...  0.772823   1.0            -0.120408
1  ENSG00000187608        ISG15  0.508770  ...  0.772823   1.0            -0.294342
2  ENSG00000131591     C1orf159  0.782443  ...  0.906011   1.0             0.084900
3  ENSG00000176022      B3GALT6  0.360561  ...  0.690359   1.0            -0.083917
4  ENSG00000169972        PUSL1  0.216335  ...  0.599048   1.0            -0.298418

[5 rows x 10 columns]
Directory "EnrichR_results/EnrichR_DRlostAmpADCyclingBHQ1AR25DRBHQ2_BCKGRND-diffrhythmsCyclingBHQ1AmpRatio25" already exists. Skipping.
{'lib': ['KEGG_2021_Human', 'GO_Biological_Process_2023', 'MSigDB_Hallmark_2020', 'WikiPathway_2021_Human', 'TRRUST_Transcription_Factors_2019', 'GO_Cellular_Component_2023', 'TRANSFAC_and_JASPAR_PWMs', 'Transcription_Factor_PPIs', 'TargetScan_microRNA_2017', 'miRTarBase_2017'], 'genes_list_file': 'enrichR_files/DR_lostAmpAD_CyclingBHQ1AR33_DRBHQ2.csv', 'background_list_file': 'diff_rhythms_CyclingBHQ1AmpRatio33.csv'}
Please select file containing list of genes (Symbols)
Path to genes list: enrichR_files/DR_lostAmpAD_CyclingBHQ1AR33_DRBHQ2.csv

Please select file containing list of BACKGROUND genes (Symbols)
Path to Background list: diff_rhythms_CyclingBHQ1AmpRatio33.csv
Background DataFrame (nrow: 254):
           Ensembl Gene_Symbols
0  ENSG00000175756     AURKAIP1
1  ENSG00000235098      ANKRD65
2  ENSG00000116251        RPL22
3  ENSG00000116288        PARK7
4  ENSG00000142676        RPL11
Background DataFrame (nrow: 2338):
           Ensembl Gene_Symbols         p  ...       BHQ  Bonf  Log_AD_CTL_ampRatio
0  ENSG00000188290         HES4  0.509391  ...  0.732026   1.0            -0.120408
1  ENSG00000176022      B3GALT6  0.360561  ...  0.629568   1.0            -0.083917
2  ENSG00000175756     AURKAIP1  0.016927  ...  0.155197   1.0            -0.334959
3  ENSG00000235098      ANKRD65  0.025834  ...  0.195471   1.0            -0.320663
4  ENSG00000160075        SSU72  0.160727  ...  0.451659   1.0             0.010688

[5 rows x 10 columns]
Directory "EnrichR_results/EnrichR_DRlostAmpADCyclingBHQ1AR33DRBHQ2_BCKGRND-diffrhythmsCyclingBHQ1AmpRatio33" already exists. Skipping.
{'lib': ['KEGG_2021_Human', 'GO_Biological_Process_2023', 'MSigDB_Hallmark_2020', 'WikiPathway_2021_Human', 'TRRUST_Transcription_Factors_2019', 'GO_Cellular_Component_2023', 'TRANSFAC_and_JASPAR_PWMs', 'Transcription_Factor_PPIs', 'TargetScan_microRNA_2017', 'miRTarBase_2017'], 'genes_list_file': 'enrichR_files/DR_cyclers_mthd2_CyclingBHQ1AR1_DRBHQ2.csv', 'background_list_file': 'diff_rhythms_method2_CyclingBHQ1AmpRatio1.csv'}
Please select file containing list of genes (Symbols)
Path to genes list: enrichR_files/DR_cyclers_mthd2_CyclingBHQ1AR1_DRBHQ2.csv

Please select file containing list of BACKGROUND genes (Symbols)
Path to Background list: diff_rhythms_method2_CyclingBHQ1AmpRatio1.csv
Background DataFrame (nrow: 54):
           Ensembl Gene_Symbols  ...      Bonf  Log_AD_CTL_ampRatio
0  ENSG00000161281       COX7A1  ...  0.020637            -0.704488
1  ENSG00000211450      SELENOH  ...  0.068645            -0.454206
2  ENSG00000164405        UQCRQ  ...  1.000000            -0.565959
3  ENSG00000019582         CD74  ...  1.000000             1.234542
4  ENSG00000177556        ATOX1  ...  1.000000            -0.541007

[5 rows x 10 columns]
Background DataFrame (nrow: 10539):
           Ensembl Gene_Symbols         p  ...       BHQ  Bonf  Log_AD_CTL_ampRatio
0  ENSG00000188976        NOC2L  0.418157  ...  0.787375   1.0             0.333543
1  ENSG00000187961       KLHL17  0.021538  ...  0.400331   1.0             0.371275
2  ENSG00000188290         HES4  0.509391  ...  0.824404   1.0            -0.120408
3  ENSG00000187608        ISG15  0.508770  ...  0.824372   1.0            -0.294342
4  ENSG00000131591     C1orf159  0.782443  ...  0.930401   1.0             0.084900

[5 rows x 10 columns]
Directory "EnrichR_results/EnrichR_DRcyclersmthd2CyclingBHQ1AR1DRBHQ2_BCKGRND-diffrhythmsmethod2CyclingBHQ1AmpRatio1" already exists. Skipping.
{'lib': ['KEGG_2021_Human', 'GO_Biological_Process_2023', 'MSigDB_Hallmark_2020', 'WikiPathway_2021_Human', 'TRRUST_Transcription_Factors_2019', 'GO_Cellular_Component_2023', 'TRANSFAC_and_JASPAR_PWMs', 'Transcription_Factor_PPIs', 'TargetScan_microRNA_2017', 'miRTarBase_2017'], 'genes_list_file': 'enrichR_files/DR_cyclers_mthd2_CyclingBHQ1AR25_DRBHQ2.csv', 'background_list_file': 'diff_rhythms_method2_CyclingBHQ1AmpRatio25.csv'}
Please select file containing list of genes (Symbols)
Path to genes list: enrichR_files/DR_cyclers_mthd2_CyclingBHQ1AR25_DRBHQ2.csv

Please select file containing list of BACKGROUND genes (Symbols)
Path to Background list: diff_rhythms_method2_CyclingBHQ1AmpRatio25.csv
Background DataFrame (nrow: 229):
           Ensembl Gene_Symbols  ...      Bonf  Log_AD_CTL_ampRatio
0  ENSG00000161281       COX7A1  ...  0.008493            -0.704488
1  ENSG00000211450      SELENOH  ...  0.028249            -0.454206
2  ENSG00000164405        UQCRQ  ...  0.442847            -0.565959
3  ENSG00000019582         CD74  ...  0.582920             1.234542
4  ENSG00000177556        ATOX1  ...  0.635387            -0.541007

[5 rows x 10 columns]
Background DataFrame (nrow: 4337):
           Ensembl Gene_Symbols         p  ...       BHQ  Bonf  Log_AD_CTL_ampRatio
0  ENSG00000188290         HES4  0.509391  ...  0.771174   1.0            -0.120408
1  ENSG00000187608        ISG15  0.508770  ...  0.771174   1.0            -0.294342
2  ENSG00000131591     C1orf159  0.782443  ...  0.904076   1.0             0.084900
3  ENSG00000176022      B3GALT6  0.360561  ...  0.689397   1.0            -0.083917
4  ENSG00000169972        PUSL1  0.216335  ...  0.596316   1.0            -0.298418

[5 rows x 10 columns]
Directory "EnrichR_results/EnrichR_DRcyclersmthd2CyclingBHQ1AR25DRBHQ2_BCKGRND-diffrhythmsmethod2CyclingBHQ1AmpRatio25" already exists. Skipping.
{'lib': ['KEGG_2021_Human', 'GO_Biological_Process_2023', 'MSigDB_Hallmark_2020', 'WikiPathway_2021_Human', 'TRRUST_Transcription_Factors_2019', 'GO_Cellular_Component_2023', 'TRANSFAC_and_JASPAR_PWMs', 'Transcription_Factor_PPIs', 'TargetScan_microRNA_2017', 'miRTarBase_2017'], 'genes_list_file': 'enrichR_files/DR_mthd2_gainAmpAD_CyclingBHQ1AR1_DRBHQ2.csv', 'background_list_file': 'diff_rhythms_method2_CyclingBHQ1AmpRatio1.csv'}
Please select file containing list of genes (Symbols)
Path to genes list: enrichR_files/DR_mthd2_gainAmpAD_CyclingBHQ1AR1_DRBHQ2.csv

Please select file containing list of BACKGROUND genes (Symbols)
Path to Background list: diff_rhythms_method2_CyclingBHQ1AmpRatio1.csv
Background DataFrame (nrow: 9):
           Ensembl Gene_Symbols
0  ENSG00000143303       RRNAD1
1  ENSG00000162843        WDR64
2  ENSG00000226479     TMEM185B
3  ENSG00000019582         CD74
4  ENSG00000196422      PPP1R26
Background DataFrame (nrow: 10539):
           Ensembl Gene_Symbols         p  ...       BHQ  Bonf  Log_AD_CTL_ampRatio
0  ENSG00000188976        NOC2L  0.418157  ...  0.787375   1.0             0.333543
1  ENSG00000187961       KLHL17  0.021538  ...  0.400331   1.0             0.371275
2  ENSG00000188290         HES4  0.509391  ...  0.824404   1.0            -0.120408
3  ENSG00000187608        ISG15  0.508770  ...  0.824372   1.0            -0.294342
4  ENSG00000131591     C1orf159  0.782443  ...  0.930401   1.0             0.084900

[5 rows x 10 columns]
Directory "EnrichR_results/EnrichR_DRmthd2gainAmpADCyclingBHQ1AR1DRBHQ2_BCKGRND-diffrhythmsmethod2CyclingBHQ1AmpRatio1" already exists. Skipping.
{'lib': ['KEGG_2021_Human', 'GO_Biological_Process_2023', 'MSigDB_Hallmark_2020', 'WikiPathway_2021_Human', 'TRRUST_Transcription_Factors_2019', 'GO_Cellular_Component_2023', 'TRANSFAC_and_JASPAR_PWMs', 'Transcription_Factor_PPIs', 'TargetScan_microRNA_2017', 'miRTarBase_2017'], 'genes_list_file': 'enrichR_files/DR_mthd2_gainAmpAD_CyclingBHQ1AR25_DRBHQ2.csv', 'background_list_file': 'diff_rhythms_method2_CyclingBHQ1AmpRatio25.csv'}
Please select file containing list of genes (Symbols)
Path to genes list: enrichR_files/DR_mthd2_gainAmpAD_CyclingBHQ1AR25_DRBHQ2.csv

Please select file containing list of BACKGROUND genes (Symbols)
Path to Background list: diff_rhythms_method2_CyclingBHQ1AmpRatio25.csv
Background DataFrame (nrow: 45):
           Ensembl Gene_Symbols
0  ENSG00000133063        CHIT1
1  ENSG00000162843        WDR64
2  ENSG00000226479     TMEM185B
3  ENSG00000144395      CCDC150
4  ENSG00000240891       PLCXD2
Background DataFrame (nrow: 4337):
           Ensembl Gene_Symbols         p  ...       BHQ  Bonf  Log_AD_CTL_ampRatio
0  ENSG00000188290         HES4  0.509391  ...  0.771174   1.0            -0.120408
1  ENSG00000187608        ISG15  0.508770  ...  0.771174   1.0            -0.294342
2  ENSG00000131591     C1orf159  0.782443  ...  0.904076   1.0             0.084900
3  ENSG00000176022      B3GALT6  0.360561  ...  0.689397   1.0            -0.083917
4  ENSG00000169972        PUSL1  0.216335  ...  0.596316   1.0            -0.298418

[5 rows x 10 columns]
Directory "EnrichR_results/EnrichR_DRmthd2gainAmpADCyclingBHQ1AR25DRBHQ2_BCKGRND-diffrhythmsmethod2CyclingBHQ1AmpRatio25" already exists. Skipping.
{'lib': ['KEGG_2021_Human', 'GO_Biological_Process_2023', 'MSigDB_Hallmark_2020', 'WikiPathway_2021_Human', 'TRRUST_Transcription_Factors_2019', 'GO_Cellular_Component_2023', 'TRANSFAC_and_JASPAR_PWMs', 'Transcription_Factor_PPIs', 'TargetScan_microRNA_2017', 'miRTarBase_2017'], 'genes_list_file': 'enrichR_files/DR_mthd2_lostAmpAD_CyclingBHQ1AR1_DRBHQ2.csv', 'background_list_file': 'diff_rhythms_method2_CyclingBHQ1AmpRatio1.csv'}
Please select file containing list of genes (Symbols)
Path to genes list: enrichR_files/DR_mthd2_lostAmpAD_CyclingBHQ1AR1_DRBHQ2.csv

Please select file containing list of BACKGROUND genes (Symbols)
Path to Background list: diff_rhythms_method2_CyclingBHQ1AmpRatio1.csv
Background DataFrame (nrow: 45):
           Ensembl Gene_Symbols
0  ENSG00000178796       RIIAD1
1  ENSG00000177954        RPS27
2  ENSG00000188010        MORN2
3  ENSG00000143947       RPS27A
4  ENSG00000135617       PRADC1
Background DataFrame (nrow: 10539):
           Ensembl Gene_Symbols         p  ...       BHQ  Bonf  Log_AD_CTL_ampRatio
0  ENSG00000188976        NOC2L  0.418157  ...  0.787375   1.0             0.333543
1  ENSG00000187961       KLHL17  0.021538  ...  0.400331   1.0             0.371275
2  ENSG00000188290         HES4  0.509391  ...  0.824404   1.0            -0.120408
3  ENSG00000187608        ISG15  0.508770  ...  0.824372   1.0            -0.294342
4  ENSG00000131591     C1orf159  0.782443  ...  0.930401   1.0             0.084900

[5 rows x 10 columns]
Directory "EnrichR_results/EnrichR_DRmthd2lostAmpADCyclingBHQ1AR1DRBHQ2_BCKGRND-diffrhythmsmethod2CyclingBHQ1AmpRatio1" already exists. Skipping.
{'lib': ['KEGG_2021_Human', 'GO_Biological_Process_2023', 'MSigDB_Hallmark_2020', 'WikiPathway_2021_Human', 'TRRUST_Transcription_Factors_2019', 'GO_Cellular_Component_2023', 'TRANSFAC_and_JASPAR_PWMs', 'Transcription_Factor_PPIs', 'TargetScan_microRNA_2017', 'miRTarBase_2017'], 'genes_list_file': 'enrichR_files/DR_mthd2_lostAmpAD_CyclingBHQ1AR25_DRBHQ2.csv', 'background_list_file': 'diff_rhythms_method2_CyclingBHQ1AmpRatio25.csv'}
Please select file containing list of genes (Symbols)
Path to genes list: enrichR_files/DR_mthd2_lostAmpAD_CyclingBHQ1AR25_DRBHQ2.csv

Please select file containing list of BACKGROUND genes (Symbols)
Path to Background list: diff_rhythms_method2_CyclingBHQ1AmpRatio25.csv
Background DataFrame (nrow: 184):
           Ensembl Gene_Symbols
0  ENSG00000142676        RPL11
1  ENSG00000130770      ATP5IF1
2  ENSG00000142937         RPS8
3  ENSG00000173846         PLK3
4  ENSG00000173660        UQCRH
Background DataFrame (nrow: 4337):
           Ensembl Gene_Symbols         p  ...       BHQ  Bonf  Log_AD_CTL_ampRatio
0  ENSG00000188290         HES4  0.509391  ...  0.771174   1.0            -0.120408
1  ENSG00000187608        ISG15  0.508770  ...  0.771174   1.0            -0.294342
2  ENSG00000131591     C1orf159  0.782443  ...  0.904076   1.0             0.084900
3  ENSG00000176022      B3GALT6  0.360561  ...  0.689397   1.0            -0.083917
4  ENSG00000169972        PUSL1  0.216335  ...  0.596316   1.0            -0.298418

[5 rows x 10 columns]
Directory "EnrichR_results/EnrichR_DRmthd2lostAmpADCyclingBHQ1AR25DRBHQ2_BCKGRND-diffrhythmsmethod2CyclingBHQ1AmpRatio25" already exists. Skipping.
         DR_cyclers_CyclingBHQ1AR1_DRBHQ2.csv         DR_cyclers_CyclingBHQ1AR25_DRBHQ2.csv 
                                            0                                             0 
        DR_cyclers_CyclingBHQ1AR33_DRBHQ2.csv        DR_gainAmpAD_CyclingBHQ1AR1_DRBHQ2.csv 
                                            0                                             0 
      DR_gainAmpAD_CyclingBHQ1AR25_DRBHQ2.csv       DR_gainAmpAD_CyclingBHQ1AR33_DRBHQ2.csv 
                                            0                                             0 
       DR_lostAmpAD_CyclingBHQ1AR1_DRBHQ2.csv       DR_lostAmpAD_CyclingBHQ1AR25_DRBHQ2.csv 
                                            0                                             0 
      DR_lostAmpAD_CyclingBHQ1AR33_DRBHQ2.csv    DR_cyclers_mthd2_CyclingBHQ1AR1_DRBHQ2.csv 
                                            0                                             0 
  DR_cyclers_mthd2_CyclingBHQ1AR25_DRBHQ2.csv  DR_mthd2_gainAmpAD_CyclingBHQ1AR1_DRBHQ2.csv 
                                            0                                             0 
DR_mthd2_gainAmpAD_CyclingBHQ1AR25_DRBHQ2.csv  DR_mthd2_lostAmpAD_CyclingBHQ1AR1_DRBHQ2.csv 
                                            0                                             0 
DR_mthd2_lostAmpAD_CyclingBHQ1AR25_DRBHQ2.csv 
                                            0 
#AD severity
system(paste0("python3 ", path_no_space, "/Python_EnrichR_for_AnalysisPipeline.py -g \"enrichR_files/strong_cogdx_diffs_AR25.csv\" -b \"enrichR_files/AD_cyclers_AR25BHQ05.csv\""))
{'lib': ['KEGG_2021_Human', 'GO_Biological_Process_2023', 'MSigDB_Hallmark_2020', 'WikiPathway_2021_Human', 'TRRUST_Transcription_Factors_2019', 'GO_Cellular_Component_2023', 'TRANSFAC_and_JASPAR_PWMs', 'Transcription_Factor_PPIs', 'TargetScan_microRNA_2017', 'miRTarBase_2017'], 'genes_list_file': 'enrichR_files/strong_cogdx_diffs_AR25.csv', 'background_list_file': 'enrichR_files/AD_cyclers_AR25BHQ05.csv'}
Please select file containing list of genes (Symbols)
Path to genes list: enrichR_files/strong_cogdx_diffs_AR25.csv

Please select file containing list of BACKGROUND genes (Symbols)
Path to Background list: enrichR_files/AD_cyclers_AR25BHQ05.csv
Background DataFrame (nrow: 64):
           Ensembl Gene_Symbols
0  ENSG00000160075        SSU72
1  ENSG00000084623        EIF3I
2  ENSG00000163875        MEAF6
3  ENSG00000162783         IER5
4  ENSG00000162817     C1orf115
File not found.
Traceback (most recent call last):
  File "/Users/henryhollis/Box Sync/Henry_stuff/AD_project/scROSMAP/Rscripts/automatic_downstream_analysis/Python_EnrichR_for_AnalysisPipeline.py", line 90, in <module>
    background = list(background_df['Gene_Symbols'])
TypeError: 'NoneType' object is not subscriptable
path = getwd()
setwd(paste0(path_to_cyclops_ordering, "downstream_output"))
path_no_space = gsub(" ", "\\ ", path, fixed = TRUE)
#create list of filenames from the enrichR_files subfolder
gene_lists = list.files("./enrichR_files/", pattern = "^(CTL_|AD_|diff_mesor)")

#run system command to call python script that uses enrichR API
sapply(gene_lists, function(x){ 
  system(paste0("python3 ", path_no_space, "/Python_EnrichR_for_AnalysisPipeline.py -g \"enrichR_files/" ,x, "\" -b \"enrichR_files/EnrichR_background.csv\""))})
python3: can't open file '/Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit/downstream_output/Python_EnrichR_for_AnalysisPipeline.py': [Errno 2] No such file or directory
python3: can't open file '/Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit/downstream_output/Python_EnrichR_for_AnalysisPipeline.py': [Errno 2] No such file or directory
python3: can't open file '/Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit/downstream_output/Python_EnrichR_for_AnalysisPipeline.py': [Errno 2] No such file or directory
python3: can't open file '/Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit/downstream_output/Python_EnrichR_for_AnalysisPipeline.py': [Errno 2] No such file or directory
python3: can't open file '/Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit/downstream_output/Python_EnrichR_for_AnalysisPipeline.py': [Errno 2] No such file or directory
python3: can't open file '/Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit/downstream_output/Python_EnrichR_for_AnalysisPipeline.py': [Errno 2] No such file or directory
AD_CTL_mthd2_cyclers_AR25BHQ1.csv           AD_cyclers_AR25BHQ1.csv           AD_cyclers_AR33BHQ1.csv          CTL_cyclers_AR25BHQ1.csv 
                                2                                 2                                 2                                 2 
         CTL_cyclers_AR33BHQ1.csv    diff_mesor_all_genes_BHQ05.csv 
                                2                                 2 
isCyclingBHQCutoff_str = str_extract(as.character(isCyclingBHQCutoff), "(?<=\\.)\\d+")

setwd("enrichR_files/")
all_DR_gene_lists = list.files( pattern = "^(DR_)")
method2_files = grepl("mthd2", all_DR_gene_lists)
DR_gene_lists_method1 = all_DR_gene_lists[!method2_files]
DR_backgrounds_method1 = paste0("diff_rhythms_CyclingBHQ", isCyclingBHQCutoff_str, "AmpRatio", str_replace(str_extract(DR_gene_lists_method1, pattern = "AR\\d+"), "AR", ""), ".csv")
DR_gene_lists_method2 = all_DR_gene_lists[method2_files]
DR_backgrounds_method2 = paste0("diff_rhythms_method2_CyclingBHQ", isCyclingBHQCutoff_str, "AmpRatio", str_replace(str_extract(DR_gene_lists_method2, pattern = "AR\\d+"), "AR", ""), ".csv")
DR_gene_lists = c(DR_gene_lists_method1, DR_gene_lists_method2)
DR_backgrounds = c(DR_backgrounds_method1, DR_backgrounds_method2)

# repeat with differential gene lists
setwd(paste0(path_to_cyclops_ordering, "downstream_output"))
mapply(function(x, y){ 
  system(paste0("python3 ", path_no_space, "/Python_EnrichR_for_AnalysisPipeline.py -g \"enrichR_files/" ,x, "\" -b \"",y , "\""))},
  DR_gene_lists, 
  DR_backgrounds)
python3: can't open file '/Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit/downstream_output/Python_EnrichR_for_AnalysisPipeline.py': [Errno 2] No such file or directory
python3: can't open file '/Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit/downstream_output/Python_EnrichR_for_AnalysisPipeline.py': [Errno 2] No such file or directory
python3: can't open file '/Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit/downstream_output/Python_EnrichR_for_AnalysisPipeline.py': [Errno 2] No such file or directory
python3: can't open file '/Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit/downstream_output/Python_EnrichR_for_AnalysisPipeline.py': [Errno 2] No such file or directory
python3: can't open file '/Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit/downstream_output/Python_EnrichR_for_AnalysisPipeline.py': [Errno 2] No such file or directory
python3: can't open file '/Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit/downstream_output/Python_EnrichR_for_AnalysisPipeline.py': [Errno 2] No such file or directory
python3: can't open file '/Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit/downstream_output/Python_EnrichR_for_AnalysisPipeline.py': [Errno 2] No such file or directory
python3: can't open file '/Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit/downstream_output/Python_EnrichR_for_AnalysisPipeline.py': [Errno 2] No such file or directory
python3: can't open file '/Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit/downstream_output/Python_EnrichR_for_AnalysisPipeline.py': [Errno 2] No such file or directory
python3: can't open file '/Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit/downstream_output/Python_EnrichR_for_AnalysisPipeline.py': [Errno 2] No such file or directory
python3: can't open file '/Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit/downstream_output/Python_EnrichR_for_AnalysisPipeline.py': [Errno 2] No such file or directory
python3: can't open file '/Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit/downstream_output/Python_EnrichR_for_AnalysisPipeline.py': [Errno 2] No such file or directory
python3: can't open file '/Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit/downstream_output/Python_EnrichR_for_AnalysisPipeline.py': [Errno 2] No such file or directory
python3: can't open file '/Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit/downstream_output/Python_EnrichR_for_AnalysisPipeline.py': [Errno 2] No such file or directory
python3: can't open file '/Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit/downstream_output/Python_EnrichR_for_AnalysisPipeline.py': [Errno 2] No such file or directory
         DR_cyclers_CyclingBHQ1AR1_DRBHQ2.csv         DR_cyclers_CyclingBHQ1AR25_DRBHQ2.csv 
                                            2                                             2 
        DR_cyclers_CyclingBHQ1AR33_DRBHQ2.csv        DR_gainAmpAD_CyclingBHQ1AR1_DRBHQ2.csv 
                                            2                                             2 
      DR_gainAmpAD_CyclingBHQ1AR25_DRBHQ2.csv       DR_gainAmpAD_CyclingBHQ1AR33_DRBHQ2.csv 
                                            2                                             2 
       DR_lostAmpAD_CyclingBHQ1AR1_DRBHQ2.csv       DR_lostAmpAD_CyclingBHQ1AR25_DRBHQ2.csv 
                                            2                                             2 
      DR_lostAmpAD_CyclingBHQ1AR33_DRBHQ2.csv    DR_cyclers_mthd2_CyclingBHQ1AR1_DRBHQ2.csv 
                                            2                                             2 
  DR_cyclers_mthd2_CyclingBHQ1AR25_DRBHQ2.csv  DR_mthd2_gainAmpAD_CyclingBHQ1AR1_DRBHQ2.csv 
                                            2                                             2 
DR_mthd2_gainAmpAD_CyclingBHQ1AR25_DRBHQ2.csv  DR_mthd2_lostAmpAD_CyclingBHQ1AR1_DRBHQ2.csv 
                                            2                                             2 
DR_mthd2_lostAmpAD_CyclingBHQ1AR25_DRBHQ2.csv 
                                            2 

FGSEA

path = getwd()
setwd(path_to_cyclops_ordering)
Warning: The working directory was changed to /Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit inside a notebook chunk. The working directory will be reset when the chunk is finished running. Use the knitr root.dir option in the setup chunk to change the working directory for notebook chunks.
if (!(dir.exists("downstream_output/fGSEA"))){
  dir.create("downstream_output/fGSEA")
  dir.create(paste("downstream_output", "fGSEA", "rnk_files", sep = "/"))
  dir.create(paste("downstream_output" , "fGSEA", "fGSEA_results", sep = "/"))
}

isCyclingBHQCutoff_str = str_extract(as.character(isCyclingBHQCutoff), "(?<=\\.)\\d+")

#create rnk files
write_rnks(path_to_cyclops_ordering, isCyclingBHQCutoff_str)
[1] "Creating rnk files for fGSEA."
#path of parent folder
setwd(path)

#Pathways, downloaded from MsigDB
pathways <- c(gmtPathways("./MsigDB_gmts_for_GSEA/c2.cp.kegg.v2023.1.Hs.symbols.gmt"), gmtPathways("./MsigDB_gmts_for_GSEA/h.all.v2023.1.Hs.symbols.gmt"))


gene_remapping_dict = read.delim("./MsigDB_gmts_for_GSEA/Human_Gene_Symbol_with_Remapping_MSigDB.v2023.1.Hs.chip", sep = '\t')

setwd(paste(path_to_cyclops_ordering,"downstream_output", "fGSEA", "rnk_files", sep = "/"))

rnk_files = list.files(pattern = ".rnk")


if(!(dir.exists(paste(path_to_cyclops_ordering, "downstream_output","fGSEA","fGSEA_results", sep = "/")))){
  dir.create( paste(path_to_cyclops_ordering,"downstream_output", "fGSEA","fGSEA_results", sep = "/") )
  dir.create( paste(path_to_cyclops_ordering,"downstream_output", "fGSEA","fGSEA_results", "plots",sep = "/"))
}

run_fgsea(rnk_files, gene_remapping_dict, pathways)
[1] "../fGSEA_results/AD_cyclers_minusLogPRanked.csv"
[1] "Used pos fGSEA score type"
Saving 7 x 7 in image
[1] "../fGSEA_results/CTL_cyclers_minusLogPRanked.csv"
[1] "Used pos fGSEA score type"
[1] "../fGSEA_results/differential_mesor_all_genes_minusLogPRanked.csv"
[1] "Used pos fGSEA score type"
[1] "../fGSEA_results/DRgenesAmpRatio1_Log(AD-CTL)ranked.csv"
[1] "Used std fGSEA score type"
[1] "../fGSEA_results/DRgenesAmpRatio1_minusLogPRanked.csv"
[1] "Used pos fGSEA score type"
[1] "../fGSEA_results/DRgenesAmpRatio25_Log(AD-CTL)ranked_method2.csv"
[1] "Used std fGSEA score type"
[1] "../fGSEA_results/DRgenesAmpRatio25_Log(AD-CTL)ranked.csv"
[1] "Used std fGSEA score type"
[1] "../fGSEA_results/DRgenesAmpRatio25_minusLogPRanked_method2.csv"
[1] "Used pos fGSEA score type"
[1] "../fGSEA_results/DRgenesAmpRatio25_minusLogPRanked.csv"
[1] "Used pos fGSEA score type"
                                           AD_cyclers_minusLogPRanked.rnk                                           CTL_cyclers_minusLogPRanked.rnk 
                  "../fGSEA_results/plots/AD_cyclers_minusLogPRanked.png"                  "../fGSEA_results/plots/CTL_cyclers_minusLogPRanked.png" 
                         differential_mesor_all_genes_minusLogPRanked.rnk                                    DRgenesAmpRatio1_Log(AD-CTL)ranked.rnk 
"../fGSEA_results/plots/differential_mesor_all_genes_minusLogPRanked.png"           "../fGSEA_results/plots/DRgenesAmpRatio1_Log(AD-CTL)ranked.png" 
                                     DRgenesAmpRatio1_minusLogPRanked.rnk                           DRgenesAmpRatio25_Log(AD-CTL)ranked_method2.rnk 
            "../fGSEA_results/plots/DRgenesAmpRatio1_minusLogPRanked.png"  "../fGSEA_results/plots/DRgenesAmpRatio25_Log(AD-CTL)ranked_method2.png" 
                                  DRgenesAmpRatio25_Log(AD-CTL)ranked.rnk                             DRgenesAmpRatio25_minusLogPRanked_method2.rnk 
         "../fGSEA_results/plots/DRgenesAmpRatio25_Log(AD-CTL)ranked.png"    "../fGSEA_results/plots/DRgenesAmpRatio25_minusLogPRanked_method2.png" 
                                    DRgenesAmpRatio25_minusLogPRanked.rnk 
           "../fGSEA_results/plots/DRgenesAmpRatio25_minusLogPRanked.png" 

Run fgsea on sleep related pathways from https://doi.org/10.1152/physiolgenomics.00275.2006

Phase Set Enrichment Analysis (PSEA)

path = getwd()
path_no_space = gsub(" ", "\\ ", path, fixed = TRUE)

# List of pathways, downloaded from MsigDB, absolute paths
psea_pathways <- c(paste0(getwd(), "/MsigDB_gmts_for_GSEA/c2.cp.kegg.v2023.1.Hs.symbols.gmt"), paste0(getwd(), "/MsigDB_gmts_for_GSEA/h.all.v2023.1.Hs.symbols.gmt"))

# make sure no spaces when calling cmd
psea_pathways = gsub(" ", "\\ ", psea_pathways, fixed = TRUE) 

# change dir into PSEA directory
setwd(paste(path_to_cyclops_ordering,"downstream_output", "PSEA_files",sep = "/"))
Warning: The working directory was changed to /Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit/downstream_output/PSEA_files inside a notebook chunk. The working directory will be reset when the chunk is finished running. Use the knitr root.dir option in the setup chunk to change the working directory for notebook chunks.
#list PSEA input files (generated during cosinor regression/DR step)
psea_files = list.files(pattern = "*.txt")

#Save path to ordering as string for cleaning up call to system
path_to_cyclops_ordering_no_space = gsub(" ", "\\ ", path_to_cyclops_ordering, fixed = TRUE)

#dataframe where rows are arguments to system, I will apply system call across rows
arg_df = expand.grid(psea_files, psea_pathways) %>% 
      mutate(out_dirs = paste(str_extract(string = as.character(Var1),     pattern= ".+(?=\\.txt)"), basename(as.character(Var2)) %>% str_extract(pattern= ".+(?=\\.gmt)")%>% 
      gsub(pattern = "\\.", replacement = ""), sep = "_") )

#apply system call, calling PSEA jar file, using rows of arg_df as arguments
apply(arg_df, 1, function(a){
    system(paste0("java -jar ", path_no_space, "/PSEA_cmd.jar ", 
    paste0(path_to_cyclops_ordering_no_space, 
    "/downstream_output/PSEA_files/", a[[1]]), " ", a[[2]]," ",
    paste0(path_to_cyclops_ordering_no_space,
    "/downstream_output/PSEA_files/", a[[3]]), " 5 10000 pdf"))
        })
Absolute Path to gene-acrophase file: /Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit//downstream_output/PSEA_files/PSEA_AD_cyclers_AR25BHQ1.txt
Absolute Path to gmt file: /Users/henryhollis/Box Sync/Henry_stuff/AD_project/scROSMAP/Rscripts/automatic_downstream_analysis/MsigDB_gmts_for_GSEA/c2.cp.kegg.v2023.1.Hs.symbols.gmt
Absolute Path to output folder: /Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit//downstream_output/PSEA_files/PSEA_AD_cyclers_AR25BHQ1_c2cpkeggv20231Hssymbols
Min geneset size: 5
Max permutations: 10000
Using output img format: pdf
Successfully loaded 1728 items
Successfully loaded 186 items
Absolute Path to gene-acrophase file: /Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit//downstream_output/PSEA_files/PSEA_CTL_cyclers_AR25BHQ1.txt
Absolute Path to gmt file: /Users/henryhollis/Box Sync/Henry_stuff/AD_project/scROSMAP/Rscripts/automatic_downstream_analysis/MsigDB_gmts_for_GSEA/c2.cp.kegg.v2023.1.Hs.symbols.gmt
Absolute Path to output folder: /Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit//downstream_output/PSEA_files/PSEA_CTL_cyclers_AR25BHQ1_c2cpkeggv20231Hssymbols
Min geneset size: 5
Max permutations: 10000
Using output img format: pdf
Successfully loaded 1673 items
Successfully loaded 186 items
Absolute Path to gene-acrophase file: /Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit//downstream_output/PSEA_files/PSEA_DR_AR25BHQ2_acrodiffs.txt
Absolute Path to gmt file: /Users/henryhollis/Box Sync/Henry_stuff/AD_project/scROSMAP/Rscripts/automatic_downstream_analysis/MsigDB_gmts_for_GSEA/c2.cp.kegg.v2023.1.Hs.symbols.gmt
Absolute Path to output folder: /Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit//downstream_output/PSEA_files/PSEA_DR_AR25BHQ2_acrodiffs_c2cpkeggv20231Hssymbols
Min geneset size: 5
Max permutations: 10000
Using output img format: pdf
Successfully loaded 516 items
Successfully loaded 186 items
Absolute Path to gene-acrophase file: /Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit//downstream_output/PSEA_files/PSEA_AD_cyclers_AR25BHQ1.txt
Absolute Path to gmt file: /Users/henryhollis/Box Sync/Henry_stuff/AD_project/scROSMAP/Rscripts/automatic_downstream_analysis/MsigDB_gmts_for_GSEA/h.all.v2023.1.Hs.symbols.gmt
Absolute Path to output folder: /Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit//downstream_output/PSEA_files/PSEA_AD_cyclers_AR25BHQ1_hallv20231Hssymbols
Min geneset size: 5
Max permutations: 10000
Using output img format: pdf
Successfully loaded 1728 items
Successfully loaded 50 items
Absolute Path to gene-acrophase file: /Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit//downstream_output/PSEA_files/PSEA_CTL_cyclers_AR25BHQ1.txt
Absolute Path to gmt file: /Users/henryhollis/Box Sync/Henry_stuff/AD_project/scROSMAP/Rscripts/automatic_downstream_analysis/MsigDB_gmts_for_GSEA/h.all.v2023.1.Hs.symbols.gmt
Absolute Path to output folder: /Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit//downstream_output/PSEA_files/PSEA_CTL_cyclers_AR25BHQ1_hallv20231Hssymbols
Min geneset size: 5
Max permutations: 10000
Using output img format: pdf
Successfully loaded 1673 items
Successfully loaded 50 items
Absolute Path to gene-acrophase file: /Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit//downstream_output/PSEA_files/PSEA_DR_AR25BHQ2_acrodiffs.txt
Absolute Path to gmt file: /Users/henryhollis/Box Sync/Henry_stuff/AD_project/scROSMAP/Rscripts/automatic_downstream_analysis/MsigDB_gmts_for_GSEA/h.all.v2023.1.Hs.symbols.gmt
Absolute Path to output folder: /Users/henryhollis/Box Sync/Henry_stuff/AD_project/human_data/Cyclops_folders/training_output/scROSMAP/cogdx_controls/wAD/ExcitatoryNeurons/TMMs_w_batch/Exc3and5_FiltByEdgeRDefault_fixedipBulkChenZhang_condAndBatchCovs_04ContrVar4EGdefault_noTransferFit//downstream_output/PSEA_files/PSEA_DR_AR25BHQ2_acrodiffs_hallv20231Hssymbols
Min geneset size: 5
Max permutations: 10000
Using output img format: pdf
Successfully loaded 516 items
Successfully loaded 50 items
[1] 0 0 0 0 0 0

Transcription Factor Analysis

Now I search for transcription factors associated with my various gene lists, i.e, TF’s associated with genes that lose amplitude in AD, etc.

Write KEGG Files

Write Files for Kegg Pathway Diagrams

Metabolon

Proteomics

Make sure you have downloaded syn32835854, syn21266449, and syn51389605. You can download programmatically with python or R synapse client. Make sure you specify the path to those downloads with proteiomics_file_path.

Plot Other Proteins

plot_list = c("BPNT1|O95861", "PPM1B|O75688", "BLVRA|P53004", "HDHD2|Q9H0R4")

proteiomics_file_path = "~/Box Sync/Henry_stuff/AD_project/scROSMAP/Proteomics/SynapseDownloads/"
tmt_filename = "2b.unlogged_relative_abundance-ROSMAP610.R1+R2-2xTAMPORcorrected.csv" #syn51389605
# tmt_filename = "3.Regressed(age.sex.PMI-noDxProtect)_cleanDat-7814x610-ROSMAP610_R1+R2(Uniprot2019multiconsensus).csv"

plot_prot(proteiomics_file_path, tmt_filename, path_to_rosmap_clin, path_to_cyclops_ordering, c("RPL10|P27635", "RPL32|P62910", "NDUFB1|O75438"), percentile = percentile)

fGSEA on proteomics

LS0tCnRpdGxlOiAiUGVyZm9ybSBEb3duc3RyZWFtIEFuYWx5c2lzIgpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sKZWRpdG9yX29wdGlvbnM6IAogIGNodW5rX291dHB1dF90eXBlOiBpbmxpbmUKLS0tCgpUaGlzIG5vdGVib29rIGlzIGZvciBwZXJmb3JtaW5nIGRvd25zdHJlYW0gYW5hbHlzaXMgb2YgQ1lDTE9QUy4gSXQgaW5jbHVkZXMgcnVubmluZyBjb3Npbm9yIHJlZ3Jlc3Npb24gYW5kIGFzc2Vzc2luZyBkaWZmZXJlbnRpYWwgcmh5dGhtaWNpdHkgKHR3byB3YXlzIGZvciBib3RoKSwgY2FsbGluZyB0aGUgZW5yaWNoUiBBUEkgdmlhIHB5dGhvbiwgdXNpbmcgZmdzZWEgaW4gUiwgcnVubmluZyBQU0VBLCBvdXRwdXRpbmcgZmlsZXMgZm9yIG1ha2luZyBLRUdHIG1hcCBkaWFncmFtcywgZXRjLgoKIyMjIExvYWQgUGFja2FnZXMKCmBgYHtyIExvYWQgTGlicmFyaWVzLCBpbmNsdWRlPUZBTFNFfQpsaWJyYXJ5KHJldGljdWxhdGUpICN1c2UgdG8gZmluZCBjb25kYSBlbnYgZm9yIHB5dGhvbgpsaWJyYXJ5KHRpZHl2ZXJzZSkKc291cmNlKCIuL2RpZmZlcmVudGlhbF9yaHl0aF9nZW5lcmFsaXplZC5SIikKc291cmNlKCIuL2NyZWF0ZV9ybmtfZmlsZXMuUiIpCnNvdXJjZSgiLi9mZ3NlYS5SIikKc291cmNlKCIuL01lc29yX1Zlbm5fRGlhZ3JhbS5SIikKc291cmNlKCIuL0ZpbGVzX2Zvcl9Qc2Nhbi5SIikKc291cmNlKCIuL01ha2VfS0VHR19tYXBfZmlsZXMuUiIpCnNvdXJjZSgiLi9DaGVja19URl9mb3JfY3ljbGluZ19vcl9EUi5SIikKc291cmNlKCIuL2Znc2VhX3Bsb3RfZm9yX2Rvd25zdHJlYW0uUiIpCnNvdXJjZSgiLi9wbG90X2Nsb2NrX2ZhY2VfcGxvdHMuUiIpCnNvdXJjZSgiLi9wbG90X2dlbmVfdHJhY2VzX2dlbmVyYWxpemVkLlIiKQpzb3VyY2UoIi4vb3JkZXJfbWV0YWJvbG9uLlIiKQpzb3VyY2UoIi4vb3JkZXJfcHJvdGVvbWljcy5SIikKYGBgCgpUaGlzIG5vdGVib29rIG1ha2VzIHVzZSBvZiBweXRob24gdmlhIGFuYWNvbmRhLiBQbGVhc2UgbWFrZSBzdXJlIHlvdSBpbnN0YWxsIGFuYWNvbmRhIHRvIHByb2NlZWQ6CgpgYGB7cn0KY29uZGFfbGlzdCgpCmBgYAoKU2VsZWN0IHlvdXIgY29uZGEgZW52aXJvbm1lbnQgb2YgY2hvaWNlLCBpbiBteSBjYXNlLCBpdGVtIDIgZnJvbSBsaXN0IGFib3ZlOgoKYGBge3J9CmNvbmRhX2xpc3QoKVtbMl1dWzJdICU+JSAKICB1c2VfY29uZGFlbnYocmVxdWlyZWQgPSBUUlVFKQpwcmludChwYXN0ZSgiVXNpbmcgY29uZCBlbnY6Iixjb25kYV9saXN0KClbWzFdXVsyXSApKQpgYGAKCiMjIyBTZXR1cAoKU2VsZWN0IHRoZSBDWUNMT1BTIG9yZGVyaW5nIHBhdGggYW5kIHBhdGggdG8gZGF0YSBmaWxlIGFuZCBwYXRoIHRvIERFc2VxMiBkaWZmZXJlbnRpYWwgZXhwcmVzc2lvbiByZXN1bHRzOgoKYGBge3J9CiNwYXRoIHRvIGN5Y2xvcHMgb3JkZXJpbmcKcGF0aF90b19jeWNsb3BzX29yZGVyaW5nID0gIi9Vc2Vycy9oZW5yeWhvbGxpcy9Cb3ggU3luYy9IZW5yeV9zdHVmZi9BRF9wcm9qZWN0L2h1bWFuX2RhdGEvQ3ljbG9wc19mb2xkZXJzL3RyYWluaW5nX291dHB1dC9zY1JPU01BUC9jb2dkeF9jb250cm9scy93QUQvRXhjaXRhdG9yeU5ldXJvbnMvVE1Nc193X2JhdGNoL0V4YzNhbmQ1X0ZpbHRCeUVkZ2VSRGVmYXVsdF9maXhlZGlwQnVsa0NoZW5aaGFuZ19jb25kQW5kQmF0Y2hDb3ZzXzA0Q29udHJWYXI0RUdkZWZhdWx0X25vVHJhbnNmZXJGaXQvIgoKI3BhdGggdG8gbm9ybWFsaXplZCBjb3VudHMKcGF0aF90b190bW1fZmlsZSA9ICIvVXNlcnMvaGVucnlob2xsaXMvQm94IFN5bmMvSGVucnlfc3R1ZmYvQURfcHJvamVjdC9odW1hbl9kYXRhL0N5Y2xvcHNfZm9sZGVycy90bW1zL3NjUk9TTUFQL2NvZ2R4X2NvbnRyb2xzL0VkZ2VSX2ZpbHRfbm9ybWVkL0V4Y1N1YnR5cGVzM2FuZDVfRmlsdEJ5RXhwckRlZmF1bHRfVE1NLmNzdiIKCiMgb3B0aW9uYWwgcGF0aCB0byBmaWxlIGNvbnRhaW5pbmcgZGlmZmVyZW50aWFsIGV4cHJlc3Npb24gcmVzdWx0cwpkZXNlcV9kZV9maWxlbmFtZSA9ICIvVXNlcnMvaGVucnlob2xsaXMvQm94IFN5bmMvSGVucnlfc3R1ZmYvQURfcHJvamVjdC9zY1JPU01BUC9zaW1wbGVfZGlmZmVyZW50aWFsX2V4cHIvRXhjM181L0V4YzNfNV9maWx0QnlFeHByX2NvbWJhdFNlcV9lZGdlUl9ERV9leGFjdFRlc3QuY3N2IgoKcGF0aF90b19yb3NtYXBfY2xpbiA9ICJ+L0JveCBTeW5jL0hlbnJ5X3N0dWZmL0FEX3Byb2plY3Qvc2NST1NNQVAvTWV0YV9kYXRhL2NsZWFuZWRfcm9zbWFwX21ldGFfY29nZHhDb25kcy5jc3YiCm1ldGFib2xvbl9maWxlbmFtZSA9ICJ+L0JveCBTeW5jL0hlbnJ5X3N0dWZmL0FEX3Byb2plY3Qvc2NST1NNQVAvTWV0YWJvbG9uL1JPU01BUF9NZXRhYm9sb25fSEQ0X0JyYWluNTE0X2Fzc2F5X2RhdGEuY3N2IgptZXRhYm9sb25fZGF0YWtleSA9ICJ+L0JveCBTeW5jL0hlbnJ5X3N0dWZmL0FEX3Byb2plY3Qvc2NST1NNQVAvTWV0YWJvbG9uL1JPU01BUF9NZXRhYm9sb25fSEQ0X0RhdGFfRGljdGlvbmFyeS5jc3YiCgojQkhRIGNvcnJlY3RlZCBwIGN1dG9mZiB0byBiZSBjb25zaWRlcmVkIGN5Y2xpbmcgaW4gYW5hbHlzaXM6CmlzQ3ljbGluZ0JIUUN1dG9mZiA9IDAuMQpwZXJjZW50aWxlID0gMC4wMjUKdXNlQmF0Y2ggPSBUCmBgYAoKIyMjIENvc2lub3IgUmVncmVzc2lvbiBhbmQgRGlmZmVyZW50aWFsIFJoeXRobWljaXR5CgpSdW4gY29zaW5vciByZWdyZXNzaW9uIGFuZCBkaWZmZXJlbnRpYWwgcmh5dGhtaWNpdHkgYW5hbHlzaXMgd2l0aCBjdXN0b20gc2NyaXB0IGFuZCB3aXRoIGNvbXBhcmVSaHl0aG1zIGxpYnJhcnkuIEkgbWFrZSB1c2Ugb2YgRW5zZW1ibCBhbmQgRW50cmV6IElEcyBzbyBJIG5lZWQgdGhlIGZvbGxvd2luZyB0cmFuc2xhdGlvbiBkaWN0aW9uYXJpZXM6CgpgYGB7cn0KIyByZWFkIGluIHRyYW5zbGF0aW9uIGRpY3QgYmV0d2VlbiBnZW5lIHN5bWJvbHMgYW5kIG90aGVyIGFsaWFzCkVuc2VtYmxfZGljdCA9IHJlYWR4bDo6cmVhZF94bHN4KCJzY1JPU01BUF9FTlNFTUJMX2RpY3QueGxzeCIpCgojcnVuIGN1c3RvbSBjeWNsaW5nIGFuZCBkaWZmZXJlbnRpYWwgcmh5dGhtaWNpdHkgYW5hbHlzaXMKcHJpbnQoIioqKlJ1bm5pbmcgY29zaW5vciBhbmQgZGlmZmVyZW50aWFsIHJoeXRobXMgYW5hbHlzaXMqKioiKQpydW5fY3ljbGluZ19hbmRfZHJfYW5hbHlzaXMocGF0aF90b19jeWNsb3BzX29yZGVyaW5nLCBwYXRoX3RvX3RtbV9maWxlLCBpc0N5Y2xpbmdTaWdDdXRvZmYgPSBpc0N5Y2xpbmdCSFFDdXRvZmYsIHBlcmNlbnRpbGUgPSBwZXJjZW50aWxlLCB1c2VCYXRjaCA9IHVzZUJhdGNoKQpgYGAKCiMjIyBQbG90dGluZwoKV2l0aCBjdXN0b20gcmVncmVzc2lvbiwgSSBoYXZlIHRvIHJlbWFrZSBjbG9jayBmYWNlIHBsb3RzLCBnZW5lcyBwbG90cywgZXRjLgoKYGBge3J9CnNldHdkKHBhdGhfdG9fY3ljbG9wc19vcmRlcmluZykKaWYgKCEoZGlyLmV4aXN0cygiZG93bnN0cmVhbV9vdXRwdXQvcGxvdHMiKSkpewogIGRpci5jcmVhdGUoImRvd25zdHJlYW1fb3V0cHV0L3Bsb3RzIikKfQoKI2ZvciBjeWNsaW5nX2luX0NUTC5jc3YKcGxvdF9jbG9ja19mYWNlKHBhc3RlMChwYXRoX3RvX2N5Y2xvcHNfb3JkZXJpbmcsICIvZG93bnN0cmVhbV9vdXRwdXQvcGxvdHMvQ1RMX2Nsb2NrX2ZhY2VfcGxvdC5wbmciKSwiZG93bnN0cmVhbV9vdXRwdXQvY29zaW5vcl9yZXN1bHRzX0NUTC5jc3YiLCBtb3VzZV9kYXRhID0gbW91c2VfZGF0YSwgQkhRX2N1dG9mZiA9IGlzQ3ljbGluZ0JIUUN1dG9mZiwgYW1wX3JhdGlvX2N1dG9mZiA9IDApCgpzZXR3ZChwYXRoX3RvX2N5Y2xvcHNfb3JkZXJpbmcpCgojZm9yIGN5Y2xpbmdfaW5fQUQuY3N2CnBsb3RfY2xvY2tfZmFjZShwYXN0ZTAocGF0aF90b19jeWNsb3BzX29yZGVyaW5nLCAiL2Rvd25zdHJlYW1fb3V0cHV0L3Bsb3RzL0FEX2Nsb2NrX2ZhY2VfcGxvdC5wbmciKSwiZG93bnN0cmVhbV9vdXRwdXQvY29zaW5vcl9yZXN1bHRzX0FELmNzdiIsbW91c2VfZGF0YSA9IG1vdXNlX2RhdGEsIEJIUV9jdXRvZmYgPSBpc0N5Y2xpbmdCSFFDdXRvZmYsIGFtcF9yYXRpb19jdXRvZmYgPSAwKQoKI3Bsb3QgY3ljbG9wcyBoaXN0b2dyYW0KcGxvdF9zdWJqZWN0X2hpc3RvZ3JhbShwYXRoX3RvX2N5Y2xvcHNfb3JkZXJpbmcsICJjb25kXzAiKQpwbG90X3N1YmplY3RfaGlzdG9ncmFtKHBhdGhfdG9fY3ljbG9wc19vcmRlcmluZywgImNvbmRfMSIpCgojR2VuZSBUcmFjaW5ncwpwbG90X2NvcmVfY2xvY2tfZ2VuZXMocGF0aF90b190bW1fZmlsZSwgcGF0aF90b19jeWNsb3BzX29yZGVyaW5nLCB1c2VCYXRjaCA9IHVzZUJhdGNoLCBwZXJjZW50aWxlID0gcGVyY2VudGlsZSwgc3BsaXRfY29uZF9wbG90cyA9IFQpCmBgYAoKUGxvdCBPdGhlciBHZW5lcwoKYGBge3J9CnBsb3RfbWUgPSBjKCJSUEwxOCIsIlJQUzciLCJSUEwyMiIsIlJQUzE4IiwiUlBMOCIsIlJQTDEwQSIsIlJQUzEzIiwiUlBMMTQiLCJSUFMxNUEiLCJSUFMxMiIsIlJQTDM2QUwiLCJSUEwzNCIsIlJQTDMwIiwiUlBMOSIsIlJQTDIxIiwiUlBTMjgiLCJSUFNBIiwiUlBMMjYiLCJSUEwyN0EiLCJGQVUiLCJSUFMzQSIsIlJQTDdBIiwiUlBMMTciLCJSUFMyN0EiLCJSUFMxNiIsIlJQTDE1IiwiUlBMUDEiLCJSUEw3IiwiUlBMMTAiLCJSUEwyM0EiLCJSUEwxOEEiLCJSUEwzIiwiUlBMMTIiLCJSUEwzOSIsIlJQTDM3QSIsIlJQTDM3IiwiUlBTMTAiLCJSUEwxMyIsIlJQTDM1QSIsIlJQUzIzIiwiUlBMMTNBIiwiUlBMMzgiLCJSUFMyOSIsIlJQTDI0IiwiUlBMNDEiLCJSUEwzNSIsIlJQUzI2IiwiUlBTMiIsIlJQUzkiLCJSUEwzMiIsIlJQTDYiLCJSUFMyMSIsIlJQUzI0IiwiUlBTMjUiLCJSUEwzNiIsIlJQTDE5IiwiUlBTMjciLCJSUEwyOSIsIlJQUzE1IiwiUlBTMyIsIlJQUzUiLCJSUFM4IiwiUlBMMjciLCJSU0wyNEQxIiwiUlBTNiIsIlJQTDExIikKCnRtbSA9IHJlYWRfY3N2KHBhdGhfdG9fdG1tX2ZpbGUsIHNob3dfY29sX3R5cGVzID0gRkFMU0UpCiAgY3ljX3ByZWRfZmlsZSA9IGxpc3QuZmlsZXMocGF0aCA9IHBhc3RlMChwYXRoX3RvX2N5Y2xvcHNfb3JkZXJpbmcsICIvRml0cy8iKSwgcGF0dGVybiA9ICcqRml0X091dHB1dF8qJykKY3ljX3ByZWQgPSByZWFkX2NzdihwYXN0ZShwYXRoX3RvX2N5Y2xvcHNfb3JkZXJpbmcsICJGaXRzIiwgY3ljX3ByZWRfZmlsZVsxXSwgc2VwID0gJy8nKSwgc2hvd19jb2xfdHlwZXMgPSBGQUxTRSkKcGxvdF9nZW5lX3RyYWNlKGN5Y19wcmVkLCB0bW0sIHBsb3RfbWUgLCBwZXJjZW50aWxlID0gcGVyY2VudGlsZSwgdXNlQmF0Y2ggPSB1c2VCYXRjaCkKYGBgCgojIyMgRW5yaWNoUgoKQ2FsbCBFbnJpY2hSIEFQSSB3aXRoIHB5dGhvbiwgc3VwcGx5aW5nIHZhcmlvdXMgbGlzdHMgZnJvbSBjeWNsaW5nIGFuZCBkaWZmZXJlbnRpYWwgcmh5dGhtIGFuYWx5c2lzLgoKYGBge3J9CnBhdGggPSBnZXR3ZCgpCnNldHdkKHBhc3RlMChwYXRoX3RvX2N5Y2xvcHNfb3JkZXJpbmcsICJkb3duc3RyZWFtX291dHB1dCIpKQpwYXRoX25vX3NwYWNlID0gZ3N1YigiICIsICJcXCAiLCBwYXRoLCBmaXhlZCA9IFRSVUUpCiNjcmVhdGUgbGlzdCBvZiBmaWxlbmFtZXMgZnJvbSB0aGUgZW5yaWNoUl9maWxlcyBzdWJmb2xkZXIKZ2VuZV9saXN0cyA9IGxpc3QuZmlsZXMoIi4vZW5yaWNoUl9maWxlcy8iLCBwYXR0ZXJuID0gIl4oQ1RMX3xBRF98ZGlmZl9tZXNvcikiKQoKI3J1biBzeXN0ZW0gY29tbWFuZCB0byBjYWxsIHB5dGhvbiBzY3JpcHQgdGhhdCB1c2VzIGVucmljaFIgQVBJCnNhcHBseShnZW5lX2xpc3RzLCBmdW5jdGlvbih4KXsgCiAgc3lzdGVtKHBhc3RlMCgicHl0aG9uMyAiLCBwYXRoX25vX3NwYWNlLCAiL1B5dGhvbl9FbnJpY2hSX2Zvcl9BbmFseXNpc1BpcGVsaW5lLnB5IC1nIFwiZW5yaWNoUl9maWxlcy8iICx4LCAiXCIgLWIgXCJlbnJpY2hSX2ZpbGVzL0VucmljaFJfYmFja2dyb3VuZC5jc3ZcIiIpKX0pCgppc0N5Y2xpbmdCSFFDdXRvZmZfc3RyID0gc3RyX2V4dHJhY3QoYXMuY2hhcmFjdGVyKGlzQ3ljbGluZ0JIUUN1dG9mZiksICIoPzw9XFwuKVxcZCsiKQoKc2V0d2QoImVucmljaFJfZmlsZXMvIikKYWxsX0RSX2dlbmVfbGlzdHMgPSBsaXN0LmZpbGVzKCBwYXR0ZXJuID0gIl4oRFJfKSIpCm1ldGhvZDJfZmlsZXMgPSBncmVwbCgibXRoZDIiLCBhbGxfRFJfZ2VuZV9saXN0cykKRFJfZ2VuZV9saXN0c19tZXRob2QxID0gYWxsX0RSX2dlbmVfbGlzdHNbIW1ldGhvZDJfZmlsZXNdCkRSX2JhY2tncm91bmRzX21ldGhvZDEgPSBwYXN0ZTAoImRpZmZfcmh5dGhtc19DeWNsaW5nQkhRIiwgaXNDeWNsaW5nQkhRQ3V0b2ZmX3N0ciwgIkFtcFJhdGlvIiwgc3RyX3JlcGxhY2Uoc3RyX2V4dHJhY3QoRFJfZ2VuZV9saXN0c19tZXRob2QxLCBwYXR0ZXJuID0gIkFSXFxkKyIpLCAiQVIiLCAiIiksICIuY3N2IikKRFJfZ2VuZV9saXN0c19tZXRob2QyID0gYWxsX0RSX2dlbmVfbGlzdHNbbWV0aG9kMl9maWxlc10KRFJfYmFja2dyb3VuZHNfbWV0aG9kMiA9IHBhc3RlMCgiZGlmZl9yaHl0aG1zX21ldGhvZDJfQ3ljbGluZ0JIUSIsIGlzQ3ljbGluZ0JIUUN1dG9mZl9zdHIsICJBbXBSYXRpbyIsIHN0cl9yZXBsYWNlKHN0cl9leHRyYWN0KERSX2dlbmVfbGlzdHNfbWV0aG9kMiwgcGF0dGVybiA9ICJBUlxcZCsiKSwgIkFSIiwgIiIpLCAiLmNzdiIpCkRSX2dlbmVfbGlzdHMgPSBjKERSX2dlbmVfbGlzdHNfbWV0aG9kMSwgRFJfZ2VuZV9saXN0c19tZXRob2QyKQpEUl9iYWNrZ3JvdW5kcyA9IGMoRFJfYmFja2dyb3VuZHNfbWV0aG9kMSwgRFJfYmFja2dyb3VuZHNfbWV0aG9kMikKCiMgcmVwZWF0IHdpdGggZGlmZmVyZW50aWFsIGdlbmUgbGlzdHMKc2V0d2QocGFzdGUwKHBhdGhfdG9fY3ljbG9wc19vcmRlcmluZywgImRvd25zdHJlYW1fb3V0cHV0IikpCm1hcHBseShmdW5jdGlvbih4LCB5KXsgCiAgc3lzdGVtKHBhc3RlMCgicHl0aG9uMyAiLCBwYXRoX25vX3NwYWNlLCAiL1B5dGhvbl9FbnJpY2hSX2Zvcl9BbmFseXNpc1BpcGVsaW5lLnB5IC1nIFwiZW5yaWNoUl9maWxlcy8iICx4LCAiXCIgLWIgXCIiLHkgLCAiXCIiKSl9LAogIERSX2dlbmVfbGlzdHMsIAogIERSX2JhY2tncm91bmRzKQoKI0FEIHNldmVyaXR5CnN5c3RlbShwYXN0ZTAoInB5dGhvbjMgIiwgcGF0aF9ub19zcGFjZSwgIi9QeXRob25fRW5yaWNoUl9mb3JfQW5hbHlzaXNQaXBlbGluZS5weSAtZyBcImVucmljaFJfZmlsZXMvc3Ryb25nX2NvZ2R4X2RpZmZzX0FSMjUuY3N2XCIgLWIgXCJlbnJpY2hSX2ZpbGVzL0FEX2N5Y2xlcnNfQVIyNUJIUTA1LmNzdlwiIikpCmBgYAoKIyMjIEZHU0VBCgpgYGB7cn0KcGF0aCA9IGdldHdkKCkKc2V0d2QocGF0aF90b19jeWNsb3BzX29yZGVyaW5nKQppZiAoIShkaXIuZXhpc3RzKCJkb3duc3RyZWFtX291dHB1dC9mR1NFQSIpKSl7CiAgZGlyLmNyZWF0ZSgiZG93bnN0cmVhbV9vdXRwdXQvZkdTRUEiKQogIGRpci5jcmVhdGUocGFzdGUoImRvd25zdHJlYW1fb3V0cHV0IiwgImZHU0VBIiwgInJua19maWxlcyIsIHNlcCA9ICIvIikpCiAgZGlyLmNyZWF0ZShwYXN0ZSgiZG93bnN0cmVhbV9vdXRwdXQiICwgImZHU0VBIiwgImZHU0VBX3Jlc3VsdHMiLCBzZXAgPSAiLyIpKQp9Cgppc0N5Y2xpbmdCSFFDdXRvZmZfc3RyID0gc3RyX2V4dHJhY3QoYXMuY2hhcmFjdGVyKGlzQ3ljbGluZ0JIUUN1dG9mZiksICIoPzw9XFwuKVxcZCsiKQoKI2NyZWF0ZSBybmsgZmlsZXMKd3JpdGVfcm5rcyhwYXRoX3RvX2N5Y2xvcHNfb3JkZXJpbmcsIGlzQ3ljbGluZ0JIUUN1dG9mZl9zdHIpCgojcGF0aCBvZiBwYXJlbnQgZm9sZGVyCnNldHdkKHBhdGgpCgojUGF0aHdheXMsIGRvd25sb2FkZWQgZnJvbSBNc2lnREIKcGF0aHdheXMgPC0gYyhnbXRQYXRod2F5cygiLi9Nc2lnREJfZ210c19mb3JfR1NFQS9jMi5jcC5rZWdnLnYyMDIzLjEuSHMuc3ltYm9scy5nbXQiKSwgZ210UGF0aHdheXMoIi4vTXNpZ0RCX2dtdHNfZm9yX0dTRUEvaC5hbGwudjIwMjMuMS5Icy5zeW1ib2xzLmdtdCIpKQoKCmdlbmVfcmVtYXBwaW5nX2RpY3QgPSByZWFkLmRlbGltKCIuL01zaWdEQl9nbXRzX2Zvcl9HU0VBL0h1bWFuX0dlbmVfU3ltYm9sX3dpdGhfUmVtYXBwaW5nX01TaWdEQi52MjAyMy4xLkhzLmNoaXAiLCBzZXAgPSAnXHQnKQoKc2V0d2QocGFzdGUocGF0aF90b19jeWNsb3BzX29yZGVyaW5nLCJkb3duc3RyZWFtX291dHB1dCIsICJmR1NFQSIsICJybmtfZmlsZXMiLCBzZXAgPSAiLyIpKQoKcm5rX2ZpbGVzID0gbGlzdC5maWxlcyhwYXR0ZXJuID0gIi5ybmsiKQoKCmlmKCEoZGlyLmV4aXN0cyhwYXN0ZShwYXRoX3RvX2N5Y2xvcHNfb3JkZXJpbmcsICJkb3duc3RyZWFtX291dHB1dCIsImZHU0VBIiwiZkdTRUFfcmVzdWx0cyIsIHNlcCA9ICIvIikpKSl7CiAgZGlyLmNyZWF0ZSggcGFzdGUocGF0aF90b19jeWNsb3BzX29yZGVyaW5nLCJkb3duc3RyZWFtX291dHB1dCIsICJmR1NFQSIsImZHU0VBX3Jlc3VsdHMiLCBzZXAgPSAiLyIpICkKICBkaXIuY3JlYXRlKCBwYXN0ZShwYXRoX3RvX2N5Y2xvcHNfb3JkZXJpbmcsImRvd25zdHJlYW1fb3V0cHV0IiwgImZHU0VBIiwiZkdTRUFfcmVzdWx0cyIsICJwbG90cyIsc2VwID0gIi8iKSkKfQoKcnVuX2Znc2VhKHJua19maWxlcywgZ2VuZV9yZW1hcHBpbmdfZGljdCwgcGF0aHdheXMpCmBgYAoKUnVuIGZnc2VhIG9uIHNsZWVwIHJlbGF0ZWQgcGF0aHdheXMgZnJvbSBodHRwczovL2RvaS5vcmcvMTAuMTE1Mi9waHlzaW9sZ2Vub21pY3MuMDAyNzUuMjAwNgoKYGBge3J9CnNsZWVwX3dha2VfcGF0aHdheXMgPSBnbXRQYXRod2F5cygiLi9Nc2lnREJfZ210c19mb3JfR1NFQS9zbGVlcF93YWtlX3BhdGh3YXlzX3dfYWxsX2h1bWFuX2FsaWFzZXNfZmlsdGVyZWRfdG9fYmFja2dyb3VuZF9nZW5lcy5nbXQiKQoKaWYoIShkaXIuZXhpc3RzKHBhc3RlKHBhdGhfdG9fY3ljbG9wc19vcmRlcmluZywgImRvd25zdHJlYW1fb3V0cHV0IiwiZkdTRUEiLCJmR1NFQV9yZXN1bHRzX3NsZWVwX3dha2UiLCBzZXAgPSAiLyIpKSkpewogICAgZGlyLmNyZWF0ZSggcGFzdGUocGF0aF90b19jeWNsb3BzX29yZGVyaW5nLCJkb3duc3RyZWFtX291dHB1dCIsICJmR1NFQSIsImZHU0VBX3Jlc3VsdHNfc2xlZXBfd2FrZSIsIHNlcCA9ICIvIikgKQogICAgZGlyLmNyZWF0ZSggcGFzdGUocGF0aF90b19jeWNsb3BzX29yZGVyaW5nLCJkb3duc3RyZWFtX291dHB1dCIsICJmR1NFQSIsImZHU0VBX3Jlc3VsdHNfc2xlZXBfd2FrZSIsICJwbG90cyIsIHNlcCA9ICIvIikgKQp9CnNldHdkKHBhc3RlKHBhdGhfdG9fY3ljbG9wc19vcmRlcmluZywiZG93bnN0cmVhbV9vdXRwdXQiLCAiZkdTRUEiLCAicm5rX2ZpbGVzIiwgc2VwID0gIi8iKSkKCnJ1bl9mZ3NlYShybmtfZmlsZXMsIGdlbmVfcmVtYXBwaW5nX2RpY3QsIHNsZWVwX3dha2VfcGF0aHdheXMsIGRpciA9ICJmR1NFQV9yZXN1bHRzX3NsZWVwX3dha2UiLCBtYXhfc2l6ZSA9IDEwMDAwKQpgYGAKCiMjIyBQaGFzZSBTZXQgRW5yaWNobWVudCBBbmFseXNpcyAoUFNFQSkKCmBgYHtyfQpwYXRoID0gZ2V0d2QoKQpwYXRoX25vX3NwYWNlID0gZ3N1YigiICIsICJcXCAiLCBwYXRoLCBmaXhlZCA9IFRSVUUpCgojIExpc3Qgb2YgcGF0aHdheXMsIGRvd25sb2FkZWQgZnJvbSBNc2lnREIsIGFic29sdXRlIHBhdGhzCnBzZWFfcGF0aHdheXMgPC0gYyhwYXN0ZTAoZ2V0d2QoKSwgIi9Nc2lnREJfZ210c19mb3JfR1NFQS9jMi5jcC5rZWdnLnYyMDIzLjEuSHMuc3ltYm9scy5nbXQiKSwgcGFzdGUwKGdldHdkKCksICIvTXNpZ0RCX2dtdHNfZm9yX0dTRUEvaC5hbGwudjIwMjMuMS5Icy5zeW1ib2xzLmdtdCIpKQoKIyBtYWtlIHN1cmUgbm8gc3BhY2VzIHdoZW4gY2FsbGluZyBjbWQKcHNlYV9wYXRod2F5cyA9IGdzdWIoIiAiLCAiXFwgIiwgcHNlYV9wYXRod2F5cywgZml4ZWQgPSBUUlVFKSAKCiMgY2hhbmdlIGRpciBpbnRvIFBTRUEgZGlyZWN0b3J5CnNldHdkKHBhc3RlKHBhdGhfdG9fY3ljbG9wc19vcmRlcmluZywiZG93bnN0cmVhbV9vdXRwdXQiLCAiUFNFQV9maWxlcyIsc2VwID0gIi8iKSkKCiNsaXN0IFBTRUEgaW5wdXQgZmlsZXMgKGdlbmVyYXRlZCBkdXJpbmcgY29zaW5vciByZWdyZXNzaW9uL0RSIHN0ZXApCnBzZWFfZmlsZXMgPSBsaXN0LmZpbGVzKHBhdHRlcm4gPSAiKi50eHQiKQoKI1NhdmUgcGF0aCB0byBvcmRlcmluZyBhcyBzdHJpbmcgZm9yIGNsZWFuaW5nIHVwIGNhbGwgdG8gc3lzdGVtCnBhdGhfdG9fY3ljbG9wc19vcmRlcmluZ19ub19zcGFjZSA9IGdzdWIoIiAiLCAiXFwgIiwgcGF0aF90b19jeWNsb3BzX29yZGVyaW5nLCBmaXhlZCA9IFRSVUUpCgojZGF0YWZyYW1lIHdoZXJlIHJvd3MgYXJlIGFyZ3VtZW50cyB0byBzeXN0ZW0sIEkgd2lsbCBhcHBseSBzeXN0ZW0gY2FsbCBhY3Jvc3Mgcm93cwphcmdfZGYgPSBleHBhbmQuZ3JpZChwc2VhX2ZpbGVzLCBwc2VhX3BhdGh3YXlzKSAlPiUgCiAgICAgIG11dGF0ZShvdXRfZGlycyA9IHBhc3RlKHN0cl9leHRyYWN0KHN0cmluZyA9IGFzLmNoYXJhY3RlcihWYXIxKSwgICAgIHBhdHRlcm49ICIuKyg/PVxcLnR4dCkiKSwgYmFzZW5hbWUoYXMuY2hhcmFjdGVyKFZhcjIpKSAlPiUgc3RyX2V4dHJhY3QocGF0dGVybj0gIi4rKD89XFwuZ210KSIpJT4lIAogICAgICBnc3ViKHBhdHRlcm4gPSAiXFwuIiwgcmVwbGFjZW1lbnQgPSAiIiksIHNlcCA9ICJfIikgKQoKI2FwcGx5IHN5c3RlbSBjYWxsLCBjYWxsaW5nIFBTRUEgamFyIGZpbGUsIHVzaW5nIHJvd3Mgb2YgYXJnX2RmIGFzIGFyZ3VtZW50cwphcHBseShhcmdfZGYsIDEsIGZ1bmN0aW9uKGEpewogICAgc3lzdGVtKHBhc3RlMCgiamF2YSAtamFyICIsIHBhdGhfbm9fc3BhY2UsICIvUFNFQV9jbWQuamFyICIsIAogICAgcGFzdGUwKHBhdGhfdG9fY3ljbG9wc19vcmRlcmluZ19ub19zcGFjZSwgCiAgICAiL2Rvd25zdHJlYW1fb3V0cHV0L1BTRUFfZmlsZXMvIiwgYVtbMV1dKSwgIiAiLCBhW1syXV0sIiAiLAogICAgcGFzdGUwKHBhdGhfdG9fY3ljbG9wc19vcmRlcmluZ19ub19zcGFjZSwKICAgICIvZG93bnN0cmVhbV9vdXRwdXQvUFNFQV9maWxlcy8iLCBhW1szXV0pLCAiIDUgMTAwMDAgcGRmIikpCiAgICAgICAgfSkKYGBgCgojIyMgCgojIyMgVHJhbnNjcmlwdGlvbiBGYWN0b3IgQW5hbHlzaXMKCk5vdyBJIHNlYXJjaCBmb3IgdHJhbnNjcmlwdGlvbiBmYWN0b3JzIGFzc29jaWF0ZWQgd2l0aCBteSB2YXJpb3VzIGdlbmUgbGlzdHMsIGkuZSwgVEYncyBhc3NvY2lhdGVkIHdpdGggZ2VuZXMgdGhhdCBsb3NlIGFtcGxpdHVkZSBpbiBBRCwgZXRjLgoKYGBge3J9CnBhdGggPSBnZXR3ZCgpCgojY3JlYXRlIHBzY2FuIGRpciBpZiBkb2VzIG5vdCBleGlzdDoKaWYoIShkaXIuZXhpc3RzKHBhc3RlMChwYXRoX3RvX2N5Y2xvcHNfb3JkZXJpbmcsICJkb3duc3RyZWFtX291dHB1dC9wc2NhbiIpKSkpewogIGRpci5jcmVhdGUocGFzdGUwKHBhdGhfdG9fY3ljbG9wc19vcmRlcmluZywgImRvd25zdHJlYW1fb3V0cHV0L3BzY2FuIikgKQp9CnNldHdkKHBhc3RlMChwYXRoX3RvX2N5Y2xvcHNfb3JkZXJpbmcsICJkb3duc3RyZWFtX291dHB1dCIpKQoKI2dyYWIgZ2VuZSBsaXN0cyBmcm9tIGVucmljaFIgZmlsZXMgZGlyZWN0b3J5CmdlbmVfbGlzdHMgPSBsaXN0LmZpbGVzKCIuL2VucmljaFJfZmlsZXMvIiwgcGF0dGVybiA9ICJeKENUTF98QURffERSXykiKQoKI3dyaXRlIHBzY2FuIGlucHV0CmxhcHBseShnZW5lX2xpc3RzLCB3cml0ZV9wc2Nhbl9pbnB1dCkKcGF0aF9ub19zcGFjZSA9IGdzdWIoIiAiLCAiXFwgIiwgcGF0aCwgZml4ZWQgPSBUUlVFKQoKc2V0d2QoInBzY2FuLyIpCnBzY2FuX2ZpbGVzID0gbGlzdC5maWxlcygiLi9wc2Nhbl9maWxlcy8iLCBwYXR0ZXJuID0gIi5jc3YkIikKCnNhcHBseShwc2Nhbl9maWxlcywgZnVuY3Rpb24oeCl7CiAgc3lzdGVtKHBhc3RlMCgicHl0aG9uMyAiLCBwYXRoX25vX3NwYWNlLCAiL3BzY2FuX3JlcXVlc3RzLnB5IC0tZmlsZSBwc2Nhbl9maWxlcy8iICx4KSl9KQoKaXNDeWNsaW5nQkhRQ3V0b2ZmX3N0ciA9IHN0cl9leHRyYWN0KGFzLmNoYXJhY3Rlcihpc0N5Y2xpbmdCSFFDdXRvZmYpLCAiKD88PVxcLilcXGQrIikKCmlmKGV4aXN0cygiZGVzZXFfZGVfZmlsZW5hbWUiKSl7CiAgc2V0d2QoInBzY2FuX3Jlc3VsdHMiKQogIHBzY2FuX3Jlc3VsdF9maWxlcyA9IGxpc3QuZmlsZXMocGF0dGVybiA9ICIuY3N2JCIpCiAgcHJpbnQoIlNlYXJjaGluZyBjeWNsaW5nIGFuZCBEUiByZXN1bHRzIGZvciBwc2NhbiBURiIpCiAgc2FwcGx5KHBzY2FuX3Jlc3VsdF9maWxlcywgYXVnbWVudF90Zl9maWxlLCBkZXNlcV9kZV9maWxlbmFtZSwgaXNDeWNsaW5nQkhRQ3V0b2ZmX3N0cikKICBzZXR3ZChwYXN0ZTAocGF0aF90b19jeWNsb3BzX29yZGVyaW5nLCAiZG93bnN0cmVhbV9vdXRwdXQvZW5yaWNoUl9yZXN1bHRzIikpCiAgZGlyZWN0b3JpZXMgPSBsaXN0LmZpbGVzKCkKICBlbnJpY2hSX1RGX3BhdGhzIDwtIGxhcHBseShkaXJlY3RvcmllcywgZnVuY3Rpb24oZGlyKSBsaXN0LmZpbGVzKGRpciwgZnVsbC5uYW1lcyA9IFRSVUUsIHJlY3Vyc2l2ZSA9IFRSVUUsIHBhdHRlcm4gPSAiLipUUkFOU0ZBQ2FuZEpBU1BBUlBXTXNfQkNLR1JORC4qXFwuY3N2fC4qVHJhbnNjcmlwdGlvbkZhY3RvclBQSXNfQkNLR1JORC4qXFwuY3N2IikpCiAgZW5yaWNoUl9URl9wYXRocyA9IGZsYXR0ZW4oZW5yaWNoUl9URl9wYXRocykgJT4lIHBhc3RlMCgiLi8iLCAuKSAgCiAgc2FwcGx5KGVucmljaFJfVEZfcGF0aHMsIGF1Z21lbnRfdGZfZmlsZSwgZGVzZXFfZGVfZmlsZW5hbWUsICAgICBpc0N5Y2xpbmdCSFFDdXRvZmZfc3RyKQogIAp9CmBgYAoKIyMjIFdyaXRlIEtFR0cgRmlsZXMKCldyaXRlIEZpbGVzIGZvciBLZWdnIFBhdGh3YXkgRGlhZ3JhbXMKCmBgYHtyfQplbnRyZXpfZGljdCA9IHJlYWR4bDo6cmVhZF9leGNlbCgic2NST1NNQVBfRU5TRU1CTF9FTlRSRVpfZGljdC54bHN4IikKCnNldHdkKHBhdGhfdG9fY3ljbG9wc19vcmRlcmluZykKaXNDeWNsaW5nQkhRQ3V0b2ZmX3N0ciA9IHN0cl9leHRyYWN0KGFzLmNoYXJhY3Rlcihpc0N5Y2xpbmdCSFFDdXRvZmYpLCAiKD88PVxcLilcXGQrIikKCmlmICghKGRpci5leGlzdHMoImRvd25zdHJlYW1fb3V0cHV0L0tFR0dfbWFwX2RpYWdyYW1zIikpKXsKICBkaXIuY3JlYXRlKCJkb3duc3RyZWFtX291dHB1dC9LRUdHX21hcF9kaWFncmFtcyIpCn0Kd3JpdGVfa2VnZ19tYXBfZmlsZXMocGFzdGUwKCJkb3duc3RyZWFtX291dHB1dC9kaWZmX3JoeXRobXNfQ3ljbGluZ0JIUSIsIGlzQ3ljbGluZ0JIUUN1dG9mZl9zdHIsICJBbXBSYXRpbzI1LmNzdiIpLCJMb2dfQURfQ1RMX2FtcFJhdGlvIiwgdHJhbnNfZGljdCA9IGVudHJlel9kaWN0LEJIUV9jdXRvZmYgPSAwLjIpCgpgYGAKCiMjIyBNZXRhYm9sb24KCmBgYHtyfQppZiAoIShkaXIuZXhpc3RzKHBhc3RlMChwYXRoX3RvX2N5Y2xvcHNfb3JkZXJpbmcsICIvZG93bnN0cmVhbV9vdXRwdXQvbWV0YWJvbG9uIikpKSl7CiAgZGlyLmNyZWF0ZShwYXN0ZTAocGF0aF90b19jeWNsb3BzX29yZGVyaW5nLCAiL2Rvd25zdHJlYW1fb3V0cHV0L21ldGFib2xvbiIpKQp9CgpvcmRlcl9tZXRhYm9sb24obWV0YWJvbG9uX2ZpbGVuYW1lLCBtZXRhYm9sb25fZGF0YWtleSwgcGF0aF90b19yb3NtYXBfY2xpbiwgcGF0aF90b19jeWNsb3BzX29yZGVyaW5nKQpgYGAKCiMjIyBQcm90ZW9taWNzCgpNYWtlIHN1cmUgeW91IGhhdmUgZG93bmxvYWRlZCBbKipzeW4zMjgzNTg1NCoqXXsudW5kZXJsaW5lfSwgWyoqc3luMjEyNjY0NDkqKl17LnVuZGVybGluZX0sIGFuZCBbKipzeW41MTM4OTYwNSoqXXsudW5kZXJsaW5lfS4gWW91IGNhbiBkb3dubG9hZCBwcm9ncmFtbWF0aWNhbGx5IHdpdGggcHl0aG9uIG9yIFIgc3luYXBzZSBjbGllbnQuIE1ha2Ugc3VyZSB5b3Ugc3BlY2lmeSB0aGUgcGF0aCB0byB0aG9zZSBkb3dubG9hZHMgd2l0aCAqcHJvdGVpb21pY3NfZmlsZV9wYXRoLioKCmBgYHtyfQpwcm90ZWlvbWljc19maWxlX3BhdGggPSAifi9Cb3ggU3luYy9IZW5yeV9zdHVmZi9BRF9wcm9qZWN0L3NjUk9TTUFQL1Byb3Rlb21pY3MvU3luYXBzZURvd25sb2Fkcy8iCnRtdF9maWxlbmFtZSA9ICIyYi51bmxvZ2dlZF9yZWxhdGl2ZV9hYnVuZGFuY2UtUk9TTUFQNjEwLlIxK1IyLTJ4VEFNUE9SY29ycmVjdGVkLmNzdiIgI3N5bjUxMzg5NjA1CiMgdG10X2ZpbGVuYW1lID0gIjMuUmVncmVzc2VkKGFnZS5zZXguUE1JLW5vRHhQcm90ZWN0KV9jbGVhbkRhdC03ODE0eDYxMC1ST1NNQVA2MTBfUjErUjIoVW5pcHJvdDIwMTltdWx0aWNvbnNlbnN1cykuY3N2IgoKI21ha2Ugc3VyZSB5b3UgaGF2ZSBkb3dubG9hZGVkIHN5bjMyODM1ODU0LCAgc3luMjEyNjY0NDksIGFuZCBzeW41MTM4OTYwNQppZiAoIShkaXIuZXhpc3RzKHBhc3RlMChwYXRoX3RvX2N5Y2xvcHNfb3JkZXJpbmcsICIvZG93bnN0cmVhbV9vdXRwdXQvcHJvdGVvbWljcyIpKSkpewogIGRpci5jcmVhdGUocGFzdGUwKHBhdGhfdG9fY3ljbG9wc19vcmRlcmluZywgIi9kb3duc3RyZWFtX291dHB1dC9wcm90ZW9taWNzIikpCn0KCm9yZGVyX3Byb3QocHJvdGVpb21pY3NfZmlsZV9wYXRoLCB0bXRfZmlsZW5hbWUsIHBhdGhfdG9fcm9zbWFwX2NsaW4sIHBhdGhfdG9fY3ljbG9wc19vcmRlcmluZykKYGBgCgpQbG90IE90aGVyIFByb3RlaW5zCgpgYGB7cn0KcGxvdF9saXN0ID0gYygiQlBOVDF8Tzk1ODYxIiwgIlBQTTFCfE83NTY4OCIsICJCTFZSQXxQNTMwMDQiLCAiSERIRDJ8UTlIMFI0IikKCnByb3RlaW9taWNzX2ZpbGVfcGF0aCA9ICJ+L0JveCBTeW5jL0hlbnJ5X3N0dWZmL0FEX3Byb2plY3Qvc2NST1NNQVAvUHJvdGVvbWljcy9TeW5hcHNlRG93bmxvYWRzLyIKdG10X2ZpbGVuYW1lID0gIjJiLnVubG9nZ2VkX3JlbGF0aXZlX2FidW5kYW5jZS1ST1NNQVA2MTAuUjErUjItMnhUQU1QT1Jjb3JyZWN0ZWQuY3N2IiAjc3luNTEzODk2MDUKIyB0bXRfZmlsZW5hbWUgPSAiMy5SZWdyZXNzZWQoYWdlLnNleC5QTUktbm9EeFByb3RlY3QpX2NsZWFuRGF0LTc4MTR4NjEwLVJPU01BUDYxMF9SMStSMihVbmlwcm90MjAxOW11bHRpY29uc2Vuc3VzKS5jc3YiCgpwbG90X3Byb3QocHJvdGVpb21pY3NfZmlsZV9wYXRoLCB0bXRfZmlsZW5hbWUsIHBhdGhfdG9fcm9zbWFwX2NsaW4sIHBhdGhfdG9fY3ljbG9wc19vcmRlcmluZywgYygiUlBMMTB8UDI3NjM1IiwgIlJQTDMyfFA2MjkxMCIsICJORFVGQjF8Tzc1NDM4IiksIHBlcmNlbnRpbGUgPSBwZXJjZW50aWxlKQpgYGAKCmZHU0VBIG9uIHByb3Rlb21pY3MKCmBgYHtyfQpwYXRoID0gZ2V0d2QoKQpzZXR3ZChwYXRoX3RvX2N5Y2xvcHNfb3JkZXJpbmcpCmlmICghKGRpci5leGlzdHMoImRvd25zdHJlYW1fb3V0cHV0L3Byb3Rlb21pY3MvZkdTRUEiKSkpewogIGRpci5jcmVhdGUoImRvd25zdHJlYW1fb3V0cHV0L3Byb3Rlb21pY3MvZkdTRUEiKQogIGRpci5jcmVhdGUoImRvd25zdHJlYW1fb3V0cHV0L3Byb3Rlb21pY3MvZkdTRUEvcm5rX2ZpbGVzIikKICBkaXIuY3JlYXRlKCJkb3duc3RyZWFtX291dHB1dC9wcm90ZW9taWNzL2ZHU0VBL2ZHU0VBX3Jlc3VsdHMiKQp9Cgppc0N5Y2xpbmdCSFFDdXRvZmZfc3RyID0gc3RyX2V4dHJhY3QoYXMuY2hhcmFjdGVyKGlzQ3ljbGluZ0JIUUN1dG9mZiksICIoPzw9XFwuKVxcZCsiKQoKI3BhdGggb2YgcGFyZW50IGZvbGRlcgpzZXR3ZChwYXRoKQoKd3JpdGVfcHJvdF9ybmtzKHBhdGhfdG9fY3ljbG9wc19vcmRlcmluZykKCiNQYXRod2F5cywgZG93bmxvYWRlZCBmcm9tIE1zaWdEQgpwYXRod2F5cyA8LSBjKGdtdFBhdGh3YXlzKCIuL01zaWdEQl9nbXRzX2Zvcl9HU0VBL2MyLmNwLmtlZ2cudjIwMjMuMS5Icy5zeW1ib2xzLmdtdCIpLCBnbXRQYXRod2F5cygiLi9Nc2lnREJfZ210c19mb3JfR1NFQS9oLmFsbC52MjAyMy4xLkhzLnN5bWJvbHMuZ210IikpCgoKZ2VuZV9yZW1hcHBpbmdfZGljdCA9IHJlYWQuZGVsaW0oIi4vTXNpZ0RCX2dtdHNfZm9yX0dTRUEvSHVtYW5fVW5pUHJvdF9JRHNfTVNpZ0RCLnYyMDIzLjIuSHMuY2hpcCIsIHNlcCA9ICdcdCcpCgpzZXR3ZChwYXN0ZTAocGF0aF90b19jeWNsb3BzX29yZGVyaW5nLCIvZG93bnN0cmVhbV9vdXRwdXQvcHJvdGVvbWljcy9mR1NFQS9ybmtfZmlsZXMiKSkKCnJua19maWxlcyA9IGxpc3QuZmlsZXMocGF0dGVybiA9ICIucm5rIikKCnJ1bl9mZ3NlYShybmtfZmlsZXMsIGdlbmVfcmVtYXBwaW5nX2RpY3QsIHBhdGh3YXlzKQpgYGAK